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sin x 


the same thing in the limiting case. 


§1. In 1937, the senior author (Tsuso1: 1937) 
introduced a direct and objective method for 
educing the underground mass distribution 
M(x, y) that will produce a gravity anomaly 
distribution 49(z, y) given along the earth’s 
surface. Of course, there can be no unique 
solution for the problem. In order that a 
definite answer may be obtained, a simplifying 
and yet plausible assumption had to be intro- 
duced in that method regarding the mass dis- 
tribution. The assumption made was that 
the mass M(x, y) is distributed on a single 
horizontal plane situated at a depth d below 
the earth’s surface. This is a justifiable as- 
sumption, because in many actual cases, the 
underground mass anomaly is caused by an 
undulation of the interface between the lighter 
and the underlying denser materials and from 
the gravimetrical point of view, this arrange- 
ment can be replaced with a good approxima- 
tion by a plane mass which is situated at the 
average depth of the interface and which is 
equal to the density difference of the two 
materials multiplied by the amplitude of undu- 
lation, so far as the undulation of interface 
is not very large compared with its average 
depth. With such an assumption as stated 
above, it was shown that if 49(x, y) is ex- 
pressed by a double Fourtmr series such as 


1) 


COs COS 
49a, Y= DVB mn sin Mx sin ny, 


the mass distribution that will produce this 
Ag is given by another double Fourier series 
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The relation between the FOURIER series method and the 


The Relation between the FOURIER Series Method and the 


a’ Method for Gravity Interpretations. 
By 
Chuji Tsusor and Yoshibumi Tomopa 


Geophysical Institute, Tokyo University, Tokyo. 


Abstract 


sin x 
= method for gravity 


interpretations has been investigated and it is shown that the two methods are representing 


such as 
M(x, Y= : 
27k? 
cos cos eo 
x 2 2Bmn sin!” gin “Y EXP Vm2+n? d, 


(2) 
where k? is the NewTonian constant of gravita- 
tion and d is the depth of the mass plane. 
In two dimensional case, (1) and (2) reduce 
themselves respectively into 


Aq(x)= Ba vy, 
Ke) = 2Bm sin Mx 


(3) 


M@)=-55 (4) 


In order to deduce M(x, y) from a given 
Ag(z, y) according to this method, one proceeds 
as follows: 

1) Ag(x, y) [49(z)] is analysed into a double 
[single] Fourier series, 

2) Each of the Fourtr coefficients in 
this series is multiplied by appropriate 

expY m?-+n?d [exp md], 

3) The Fourtmr series having these new 
coefficients is synthesized, 

4) The value of the series at x=a, y=yY, 
[v=] is divided by 2zk’. 

5) The final result is the required M(x, y) 
[M(a)]. 

Although this Fourtsr series method for 
gravity interpretation has proved itself to be 
very useful and competent in applications 
(Tsupot, 1938, 1939, 1940, 1941, 1942, 1948, 
1950), admittedly it has a kind of awkwardness 


~>Bm 28 maz exp md. 
ae SINE 
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coming from the very nature of the Fourier 
series. As in many other instances in which 
it is used, the Fourtmr series expression of a 
quantity implies periodic repetition of the 
same pattern of distribution of that quantity 
even outside the domain in which we are inte- 
rested. This circumstance has caused a kind 
of hindrance for unconditional usage of the 
method. In order to get over this difficulty, 
T. Rixitake (1952) suggested the usage of 
Hermite function while Y. Sarto (1954) proposed 
a modification of the Fourtmr series method 
for this purpose. 


§2. In 1955, the junior author (Tomopa) and 
and K. Axt(Tomopa and Axt, 1955) introduced 
the sina/a method for gravity interpretations 
with the hope of getting rid of the said awk- 
wardness of the Fourier series method. This 
function sinz/z is characteristic in that it 
takes a unit value at x=0 and vanishes at all 
other grid points which are z apart. It is 
therefore fitted for representing the distribution 
of such a quantity as vanishes at all equally 
spaced grid points except at one of them 
where it takes a finite value. If applied to 
gravity problems, 


sin x 


OT NCD a pai (Ga) 


represents such a distribution of 4g(z) as it is 
b at «=0 and is 0 at e=nz, where m can be 
any integer other than zero. Since 

sin v 

eer (6) 
it is readily seen that no harmonic components 
having higher frequecies than 1 are contained 
in this function. Making use of the formula 
(3) and (4), the mass distribution M(x) that 


will give rise to the gravity anomaly distribu- 
tion: 


if 
-| cos ma dm, 
0 


Agx)=b sin ca 
wv 
is given by 
b i! 
Ma)= ais j.c0s mx exp mddm. Hey 
This can be integrated into 
M(«)=— o.. Se 
Qk? a+ d2 
x{(dcosx+asinxz)expd—d}. (8) 


| 


: 


Since sinz is zero and cosa is either +1 0 
—l at x=nz according as m is even or od 
M(x) at these points becomes simply 


b d 
2k? (nz)?+da? 
In this formula, d is expressed in radian i 
case when the distance between two consecu. 
tive grid points is takenasz. If this spacing 
is a and the depth is 6 in the actual scale, we 
may write | 


M(nz)= (+ exp d—1). 


| 
(9 


pune (10 
a 
Inserting this in (9), we get | 
b 1 d/a Ox ) 
Bee eo —1). 
a. 2xk? x n?+(d/a)? (exp a 
(11 


The numerical values of 


On, 0/a)= : : ae “ (zeexp oe =i (12 
z n?+(0d/a)? a 
for 6/a=0.5 and 1 are given in Table I. 
Table I. 
rea! O/a On 
Values of O(n, d/a)=— ALE. exp, 1) 
+n 6/a=0.5 0/a=1.0 
0 2.42583 7.04670 
1 — 0.73979 —3.84223 
2 0.14270 1.40948 
3 —0.10000 — 0.76840 
4 0.03730 0.41447 
5 — 0.03661 —0.29548 
6 0.01673 0.19041 
7 — 0.01877 —0.15378 
8 0.00945 0.10849 
) —0.01139 — 0.09367 
10 0.00605 0.06908 
11 —0.00761 —0.06301 
12 0.00419 0.04871 
13 —0.00546 —0.04514 
14 0.00309 0.03587 
15 —0.00413 — 0.03428 
16 0.00236 0.02745 
17 —0.00320 — 0.02655 
18 0.00187 0.02170 
19 — 0.00256 — 0.02124 
20 0.00149 0.01749 
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Any gravity anomaly distribution can be 
regarded as the linear superposition of 
by sin 2/2, 
with the origin of 2x shifted appropriately 
for each /. In a single mathematical expres- 
sion, the whole distribution can be written as 


dot a\ex Sb sin (x =o) 


CTL 


(13) 


where 0, is the anomaly at e=2;. Consequently 
the underground mass which is needed for 
producing this anomaly distribution will take 
the value 


1 
Me«x,)= ~ 
=H) 2nk? 
O/a On 
b ea ae << Lo < a agi a 
x 2 (l—m)?+ (6/a)? (exp a 1) 
1 a 
rey Dbi- -nmD(l—m, 36/a) (14) 
beneath w=2v. 


§3. Thus there have been introduced two 
methods for the same purpose. The question 
naturally arises what the relationship between 
the two methods, the Fourrmr series method 
and the sinz/z method, is. In what follows, 
it will be shown that the two methods are 
after all representing the same thing in the 
limiting case. 

The finite Fourter series expression of a 
function which takes a unit value at x«=0 
(w=2z) and which vanishes at all the (2p—1) 
equally spaced intervenient points lying be- 
tween x=0 and x=2z is as follows: 


fa=5 5 +008 f2Pe0s 26-+-- >> 


+408 (b—Die-+ 5008 pi}. (15) 


If this function is regarded to be representing 
a gravity anomaly distribution, the mass dis- 
tribution which will give rise to it is given by 


Ma)= ns s{gteos 2expd 
+cos 2a exp 2d+-::- 


----+cos (p—1)a exp (p—l)d 
+5008 px exp od ’ 


the mass being assumed to be condensed on 
a single horizontal plane at a depth d as be- 


(16) 


fore. The series (16) can be written as 


Me)=—+ = LB | L+exp (d+ iz) 


+exp 2(d+ix)++++- 
vee Oxp {(p—1 d+ ta) } 
dab ' 
wae exp p(d + is) 
where . means the real part of the expres- 
sion within the brackets which follow _ it. 
Since the terms in the brackets but the last 
two form a geometrical series with a com- 
mon ratio exp(d+zz), (17) can be written as 
il 1 1— exp p(d+iz) 
Mette alt 
one pe | Token aes 
al 
gtr exp na-+ie) | 
If d and x are both small as compared with 
1, this expression is transformed as follows: 


mae 


(17) 


(18) 


1—(cos pe+isin px) exp pd 


eeacn (l+d+iz) 


1| —d+(dcos px +asin px)exppd 
+x? 


il 
ak cos pa exp pa | 


(19) 


At the grid points where z is integral multiple 
of 27/2p, say 
2t NT 
s=n——_=—_ , 
Lp? 
sin px is zero and cos px is either 1 or —1 
according as ” is even or odd. Consequently, 


we get 


ate eae 


5 exp pa | 


(20) 


(21) 


In this expression, d is expressed in radian. 
If a is the actual distance between two con- 
secutive grid points, 2z in radian corresponds 
to 2pa in the actual scale. So, if @ is the 
actual depth of the mass plane, there is a 
relation 


ces (22) 
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Putting (22) into (21) and making p very large, 
we get finally 


p Qnk? 
where 
To reg ory (23 
OG, Ciay= -+(d/a? (+ exp ) (23) 


This is exactly nothing but the formula (11), 
which was derived already. Thus it has been 
proved that the sin a/z method may be regarded 
to be the limiting case of the Fourier series 
method in case the harmonic iy ee in 


Mise Fey oxy Vie 
+7008 p5-y exp V0 +p dy 


the series is taken up to a very high one. 


§4. This conclussion is useful if we wish t 
extend the the sinz/z method so as to be a 
plied to three dimensional gravity interpretae 
tions. In three dimensional cases, we have 


to evaluate the integral ) 


HEL : 
O= J,cos ma cos ny expV m?-+-n? d dm dn, — 
0 


corresponding to (7) for two dimensional cases: 
This expression does not appear to be inte: 
grable analytically. From the conclusion obtain 
ed in the above, this integral can safely be 
replaced by a double Fourier series such as; 


i 3 T ti 27 Tee a 
=—_ a ii (ane: ——y exp— V 02.22 d-+----- 
e fee cee d+cos”, y EXD 423 d+ 


1 = : 
+0055 | exp SVE+0 d+cos—y exp5VL+P d+---- 


p 
il ‘1 7 — 
+-5c0s ba yexp5V1 +p a 
+.---- 


(25: 


pees V (p—1)? rte le tear —V(p—1P+Pd+:- 


1 ua 1 
++-ycos bp, Abs exp ap +03 d+cos— 


1 Tt — 
+7 C0s bY exp Vpi-+d a} 


This function, if divided by 2zk?, represents 
a mass distribution at a depth d which will 
produce a unit gravity at x=0, y=0 and zero 
gravity at all other two dimensional grid 
points. This double Fourtmr series can be 
evaluated numerically (Tsuso1, OLDHAM, WaITH- 
MAN: 1958). 
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Numerical Tables Facilitating Three Dimensional 


Gravity Interpretations. 
By 
Chuji Tsusor 
Geophysical Institute, Tokyo University, Tokyo 
C. H. G. OtpHam* and V. B. WarrHMan 
California Research Corporation, La Habra, California, U.S.A. 


Abstract 


This is an extension of the sin x/x method which was introduced by Y. TOMODA and 


K. AKI for two dimensional gravity interpretations. 


Numerical tables facilitating three 


dimensional gravity interpretations are presented. 

Gravity values given at two dimensional grid points on the earth’s surface are to be 
multiplied by corresponding tabular values and the products added up. The sum when 
divided by 2nG will represent the surface mass density at the assumed depth and im- 


mediately beneath the central grid point. 


$1. The sin z/z method which was intro- 
duced by Y. Tomopa and K. Aki (1955) 
proved to be so useful a method for the pur- 
Dose of two dimensional gravity interpretations 
Tsupor: 1956, 1957) that its extension to three 
dimensional cases appeared highly desirable. 
As was stated elsewhere (Tsusor and Tomopa: 
1958), the evaluation of the integral 


Af igi 
Dx, Y, a=| | COS ma Cos ny 


0J0 


-expV m?+n?ddmdn (1) 


for various values of z, y and d is needed 
for this purpose. This function when divid- 
ed by 2zG represents the surface mass 
Jensity at square grid points over a horizon- 
‘al underground plane which is situated at a 
jlepth d below the surface. This is the 
surface density which will produce unit 
Sravity at xz==0, y=0 and zero gravity at all 
ther conformal grid points on the earth’s 
surface. 
Unlike the corresponding integral 


Dx, d)=\'cos maz exp md dm (CW), 


0 


n two dimensional cases, the expression (1) 
loes not appear to be integrable analytically 


* Now with the California Exploration Com- 
yany, San Francisco, California, U.S.A. 


and numerical calculations have to be resort- 
ed to for its evaluation. The necessary 
numerical integration was carried out on a 
Datatron digital computor in Pasadena, 
California, U.S.A. 


§2. The values of the function (1) which 
have been computed by means of this machine 
from Oth grid point up to 15th for both x 
and y are given in Table I and II for d=0.5 
and d=1.0 respectively. Each of the tables 
represents only one quadrant of the complete 
rectangular table which has a symmetry as 
shown in Fig. 1. 


y=15 


Fig. 1. Shaded portion is given in Table I and II. 


Table I is for d=0.5, that is for the case 
in which the mass plane is assumed to be 
situated at a depth equal to half the horizon- 
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Table I. Tsusor-OLpHAM 


0 1 3 4 5 6 
0 +3.65742 —0.74169 +0.19159 —0.09485 +0.05113 — 0.03430 -+0.02309 
1 —0.74169 +0.02208 —0.01953 +0 .00340 —0.00441 +0.00153 — 0.00182 
Z +0.19159 —0.01953 -+-0.00040 —0.00335 +-0.00009 —0.00112 +0.00012 
3 —0.09485 +0.00340 —0.00335 —0.00025 — 0.00106 —0.00012 — 0.00045 
4 +0.05113 —0.00441 +0.00009 —0.00106 —0.00019 — 0.00046 —0.00010 
5 — 0.03430 +0.00153 —0.00112 —0.00012 — 0.00046 — 0.00012 — 0.00024 
6 +-0.02309 —0.00182 +0.00012 —0.00045 —0.00010 —0.00024 —0.00008 
7 —0.01750 +0 .00086 —0.00052 —0.00004 — 0.00023 — 0.00007 — 0.00014 
8 -+-0.01309 —0.00098 +0.00010 —0.00024 — 0.00005 — 0.00013 —0.00005 
9 —0.01059 +0.00055 —0.00029 —0.00001 — 0.00013 — 0.00004 — 0.00008 
10 +0.00843 —0.00061 +0.00008 — 0.00014 — 0.00002 —0.00008 — 0.00003 
iu —0.00712 +0.00038 —0.00019 0 — 0.00008 —0.00002 — 0.00005 
12 +0.00593 —0.00042 +-0.00006 —0.00009 —0.00001 —0.00005 — 0.00002 
13 —0.00519 +0.00028 —0.00013 +0.00001 — 0.00006 — 0.00001 — 0.00004 
14 +0.00450 —0.00032 -+0.00005 — 0.00007 — 0.00001 — 0.00004 —0.00001 
5) —0.00412 +0.00023 — 0.00010 +0.00001 — 0.00004 —0.00001 — 0.00003 


i 


Table II. Tsuso1-OLpHAm 


0 1 7 3 + 5 6 
— : ane 3 
0 +15.78495 —5.84713 +-2.18517 — 1.08956 +0.63200 —0.41449 +0.28900 
1 — 5.84713 +1.35318 — 0.48663 +0.22434 —0.13492 +0.08545 —0.06161 
2 1- 2.18517 — 0.48663 +0.15329 — 0.07843 +0.04144 — 0.02851 +0.01833 
3 — 1.08956 +0.22434 —0.07843 +0.03557 =) 02328 +0.01390 — 0.01107 
4 +. 0.63200 —0.13492 +0.04144 — 0.02328 +0.01042 — 0.00756 +0.00420 
5 ~ 0.41449 +0.08545 —0.02851 +-0.01390 — 0.00756 +0.00543 — 0.00465 
6 t 0.28900 —0.06161 +0.01833 —0.01107 +0.00420 — 0.00465 +0.00165 
if | — 0.21440 +0.04485 — 0.01408 +0.00771 — 0.00346 +0.00315 — 0.00145 
8 | }- 0.16400 — 0.03535 +0.00998 — 0.00668 +-0.00200 — 0.00297 -+-0.00067 
9 — 0.13029 -+0.02781 — 0.00809 +0.00512 — 0.00174 +0.00222 — 0.00066 
10 + 0.10541 — 0.02315 +-0.00606 — 0.00464 +0.00098 —0.00217 -+-0.00022 
11 — 0.08761 +-0.01917 — 0.00508 +0.00381 — 0.00087 +0.00175 — 0.00025 
12 + 0.07381 — 0.01660 +0.00395 — 0.00355 +0.00043 —0.00175 — 0.00001 
13 — 0.06367 +0.01432 — 0.00340 +0.00306 — 0.00038 +0.00149 —0.00002 
14 | + 0.05575 — 0.01287 +0.00273 — 0.00294 +0.00012 — 0.00150 —0.00015 
15 — 0.05025 +0.01160 — 0.00247 +0.00264 —0.00011 +-0.00134 +0.00011 


eee 


F 
+g 
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i 8 9 10 iat 12 13 14 15 


—0.01750 +0.01309  -0.01059 +0.00843 -—0.00712 -+0.00593 —0.00519 +0.00450 —0.00412 


+0.00086 —0.00098 +0.00055  -0.00061 +0.00038 —0.00042 -+0.00028 —0 .00032  -+-0.00023 
— —9.00052 +-0.00010 —0.00029 +0.00008 -—0.00019 +0.00006 —0.00013 +0.00005 — 0.00010 
—0.00004 —0.00024 -—0.00001  —0.00014 0 —0.00009 +0.00001 -—0.00007 +0.00001 
—0.00023  —0.00005 -—0.00013 -—0.00002 -—0.00008 -—0.00001 -—0.00006 -—0.00001 —0.00004 
—0.00007 -—0.00013 —0.00004 -—0.00008 —0.00002  -—0.00005 -—0.00001 -—0.00004 ~—0.00001 


—0.00014 —0.00005 —0.00008 -—0.00003 -—0.00005 -—0.00002 -—0.00004  -—0.00001 —0.00003 
—0.00005 —0.00009 —0.00004 -—0.00006 -—0.00002 -—0.00004 -—0.00002 -—0.00003 —0.00001 
—0.00009 —0.00004 —0.00006 —0.00003 -—0.00004 —0.00002 -—0.00003 -—0.00001 -—0.00002 
—0.00004 —0.00006 -—0.00003 -—0.00004 -—0.00002  -—0.00003 -—0.00001 -—0.00002 —0.00001 
—0.00006 —0.00003 —0.00004 -—0.00002 -—0.00003  -—0.00002 -—0.00002 -—0.00001 —0.00002 


—0.00002 -—0.00004 —0.00002 -—0.00003 -—0.00002 —0.00002 -—0.00001  -—0.00002 —0.00001 
—0.00004 —0.00002 —0.00003 -—0.00002 -—0.00002 —0.00001 -—0.00002 —0.00001 -—0.00001 
—0.00002  —0.00003  —0.00001 —0.00002 -—0.00001 -—0.00002  -—0.00001  -—0.00001 -—0.00001 
—0.00003  —0.00001 -—0.00002 -—0.00001 -—0.00002  -—0.00001 -—0.00001 -—0.00001 -—0.00001 
—0.00001 —0.00002 —0.00001 ~0.00002 —0.00001 —0.00001 —0.00001  —0.00001 —0.00001 


WaituMaNn Coefficients for d=1.0. 


oe Ee eS 


7 8 9 10 me 12 13 14 15 


—0.21440 +0.16400 —0.13029 +0.10541 —0.08761 +0.07381 —0.06367 +0.05575 -—0.05025 


+0.04485 -—0.03535 +0.02781 -—0.02315 +0.01917 -—0.01660 +0.01432 —0.01287 +0.01160 
—0.01408 +0.00998 -—0.00809 -+-0.00606  —0.00508 +0.00395 —0.00340 +0.00273 -—0.00247 
+0.00771 -—0.00668 +0.00512 -0.00464 +0.00381 -—0.00355 +0.00306 -—0.00294 +0.00264 
“— 0.00346 +0.00200 —0.00174 +0.00098 —0.00087 +0.00043 —0.00038 +0.00012 -—0.00011 
+0.00315 —0.00297 +0.00222 -—0.00217 +0.00175 -—0.00175 +0.00149 —0.00150 +0.00134 


—0.00145 +0.00067 —0.00066 +0.00022 -—0.00025 -—0.00001 -—0.00002 -—0.00015 +0.00011 
+0.00160 -—0.00158 +0.00113  -—0.00116 +0.00089 -—0.00093 +0.00076 —0.00080 +0.00069 
~—0.00158 +0.00049 —0.00052 +0.00024 -—0.00029 +0.00011  -~0.00015 +0.00003 —0.00008 
+0.00113 -—0.00052 +0.00043  -—0.00047 -+0.00029 -—0.00033 +0.00021 -—0.00025 +0.00017 
—0.00116 +0.00024 -—0.00047 +0.00048 -—0.00053 +0.00040  -—0.00044 +0.00035 —0.00039 


+-0.00089 —0.00029 +0.00029 —0.00053 -—0.00021 +0.00017 -—0.00027 +0.00023 —0.00029 
—~0.00093 +0.00011 -—0.00033 +0.00040 +0.00017 +0.00072  -0.00076 +0.00069 —0.00072 
+0.00076 —0.00015 +0.00021 -—0.00044 -—0.00027 -—0.00076  -—0.00066 +0.00063 —0.00068 
—0.00080 +0.00003- —0.00025 +0.00035 +0.00023 +0.00069 +0.00063 +0.00098 —0.00100 
+0.00069 —0.00008 -+0.00017 -—0.00039  -0.00029  -—0.00072 0.00068 -—0.00100 -—0.00096 


a ee ee a 


a 
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Table IV. Sum of 
Column 
‘ yf c 
i oga x=0 cl i ts c=4 t=—5 26 
jc (0) 3.65742 —0.74169 0.19159 —0.09485 0.05113 —0.03430 0.02309 
i=l 2.17404 — 69753 15253 — 8805 4231 —- 3124 1945 
Vad ZeaDiee — 73659 15333 - 9485 4249 ° — 3348 1969 
Vie: EOONOS — 72979 14663 - 9525 4037 ~ 3372 1879 
y= 4 2.46978 — 73861 14681 =. 9737 3999 _ 3464 1859 
v=o 2.40118 — 73555 14457 _ 9761 3907 ~ 3488 1811 
y= 6 2.44736 — 73919 14481 — 9851 3887 _ 3536 1795 
= Te 2.41236 — 73747 14377 _ 9859 3841 — 3550 1767 
Y= 8 2.43854 — 73943 14397 = 9907 3831 ~ 3576 1757 
y= 9 2.41736 — 73833 14339 - 9909 3805 ~ 3584 1741 
y= 10 2.43422 — 73955 14355 — 9937 3801 — 3600 1735 
ih | 2.41998 — 73879 14317 — 9937 3785 — 3604 1725 
Y=12 2.43184 — 73963 14329 — 9955 3783 _ 3614 LQ 
y=13 | 2.42146 — 73907 14303 - 9953 3771 _ 3616 Las 
y=14 | 2.43046 — 73971 14313 — 9967 3769 oa 3624 weal 
y= ls BeAL ee — 73925 14293 _ 9965 3761 _ 3626 1705 
Corresponding 2.42583 — 73979 14270 — 10000 3730 — 3661 1673 


Coefficients 


eS 


Table V. Sum of 
eee 


Column 
Y up to Ss SS SS 
«=0 r= T=2 Cat aneeye C5 c= 
= (0) | 15.78495 —5.84713 2.18517 =—1.08956 0.63200 -—0.41449 0.28900 
y= 1 4.09069 -—3.14077 1.21191 — 0.64088 63784 — 24359 16578 
sv | 8.16403  -—4.11403 1.51849 — 79774 44504 — 30061 20244 
y= 3 | 6.28191 —3.66535 1.36163 — 72660 39848 To 18030 
=a 7.54591 —3.93519 1.44451 = 77316 41932 — 28793 18870 
U=—5 6.71693 -—3.76429 1.38749 — 74536 40420 oe VAM AOYE 17940 
T= 6 7.29493 —3.88751 1.42415 — 76750 41260 — 28637 18270 
y= 7 6.86613 -—3.79781 1.39599 — 75208 40568 = ~ 28007 17980 
y= 8 |} 7.19413 —3.86851 1.41595 — 76544 40968 — 28601 18114 
v= 9 | 6.93355 —3.81289 1, 39977 — 75520 40620 — 28157 17982 
Y= 10 7.14437 =—3.85919 1.41189 — 76448 44816 — 28591 18026 
i 6.96915 -3.82085 1.40173 — 75686 40642 — 28241 17976 
y= 12 7.11677 =—3.85405 1.40963 — 76396 40728 — 28591 17974 
Bole 6.98943 -—3.82541 1.40283 =) 75784 40652 — 28293 17970 
y=14 7.10093 -—3.85115 1.40829 — 76372 40676 — 28593 17940 
y=15 7.00043 —3.82795 1.403885 — 75844 40654 — 28325 17962 
sebtaepue o 7.04670 -—3.84223 1.40948 — 76840 41447 — 29548 19041 
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Values in Table I. 
a an la gs et 


Column 
Bail, G=8 s=9 c=10 Sees lih Giles T=13 e=14 aN) 
—0.01750 0.01309 —0.01059 0.00843 —0.00712 0.00593  -—0.00519 0.00450  -—0.00412 
— 1578 alain _ 949 721 = 636 509, — 463 386 — 366 
— 1682 1113 — 1007 737 - 674 521 - 489 396° = 386 
— 1690 1085 — 1009 709 - 674 503. — 487 382 = 384 
— 1736 1075 <= 41.085 705 ~ 690 501 — 499 380 392 
— 1750 1049 — 1043 689 — 694 491 ~ 501 3/2 394 
— 1778 1039 — 1059 683 - 704 ASG) 509 3/0 Re 400 
— 1788 1021 — 1067 671 — 708 479 — 513 364 — 402 
— 1806 1031 — 1079 665 — 716 475 = 519 3620 = A06 
— 1814 1001 — 1085 657 _ 720 469 — 521 308hs = 408 
— 1826 995 — 1093 653 - 726 465 — 525 356 — 412 
— 1830 987 — 1097 647 — 730 461 = 527 352 — 414 
— 1838 983 — 1103 643, — 734 459 — 531 350 — 416 
— 1842 977 — 1105 639 - 736 A550 = 533 348 418 
— 1848 975 — 1109 637 — 740 A53h) 535 346 — 420 
— 1850 971 =, » Atel 633 - 742 451 — 537 344. — 422 
— 1877 945 = * 1139 605 - 761 419 — 546 309 = 413 


Values in Table II. 


Column 
Avil x=8 co x=10 = ue GAZ T= 13, w=14 c=15 
—0.21440 0.16400 —0.13029 0.10541 —0.08761 0.07381 — 0.06367 0.05575 —0.05025 
— 12470 9330 — 7467 5911 — 4927 4061 — 3503 3001 — 2705 
— 15286 11326 — 9085 TAUZB — 5943 4851 —. 4183 3547 — 3199 
— 13744 9990 — 8061 6195 — 5181 4141 — 3571 2959 — 2671 
— 14436 10390 — 8409 6391 — 5355 4227 — 3647 2983 — 2693 
— 13806 9796 — 7965 5957 — 5005 3877 | — saad 2683 — 2425 
— 14096 9930 — 8097 6001 — 5055 3875 = 3353 2653 — 2403 
— 13776 9614 — 7871 5769 — 4877 3689 — 8201 2493 Ee 200) 
— 14092 9712 — 7975 5817 — 4935 3711 —eewozal 2499 — 2281 
— 13866 9608 — 7889 5723 — 4877 3645 — $189 2449 = 2247 
— 14098 9656 — 7983 5819 — 4983 S725 oer el 2519 — 2325 
— 13920 9598 — 7925 5713 — 5025 3759 — 3331 2565 — 2383 
— 14106 9620 — 7991 5793 — 4991 3903 = e453 2703 = 2527 
— 13954 9590 — 7949 5705 — 5045 3751 — 3615 2829 — 2663 
— 14114 9596 — 7999 5775 =e 4999 3889 — 3489 3025 — 2863 
— 13976 9580 — 7965 5697 — 5057 8745) =)) 8625 2825 — Bn) 
— 15378 10849 — 9367 6908 — 6301 4871 — 4514 3587 — 3428 


oe I ae ee 
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tal spacing between grid points. Likewise, 
Table II is for d=1.0, that is for the case 
in which the plane is assumed to be situated 
at a depth equal to the grid spacing. 

The usage of the table is simple in princi- 
ple. Prepare transparent stencils, say A, on 
which values given in Table I and II respec- 
tively are written down in a _ complete 
rectangular mesh form. Write down the 
observed gravity values (in mgals) to be 
investigated on a sheet of paper, say B, with 
the same spacing as A. Take whichever ap- 
propriate one of A’s as the case may be and 
place it on B so that the mesh on both will 
coincide. Calculate all the products of values 
which appear at the same grid points and add 
them up. The sum, when divided by 2zx 
6.67x10-" will represent the surface mass 
density in gr/cm? at the assumed depth and 
immediately beneath the central point of A 
(c=0, y=0). Shift A by one step relative to 
B and make the same calculation. Repeat 
this until the area of interest will completely 
be covered. 


§3. This method of projecting the gravity 
field downwards, or alternatively of comput- 
ing a surface mass distribution, is ideally 
suited for solution on an electronic computor. 
If a computor is not available, then the 
arithmetic labor involved is prohibitive, and 
some simplifying modification of the above 
process is needed. Such a modification must 
involve a compromise between accuracy and 
the practicability of computing the results. 
One way to shorten the calculations is to use 
only a portion of the tables. For example, 
we have found that by using the portion of 
the tables from #=0, y=0 to e=9, y=9, we 
achieve a reasonable compromise. If. still 
less accuracy can be tolerated in the final 
results, the tables can be shortened even 
further. 


§4. From the potential theory, it is evident 
that Ox, y, d) must add up to 1 if the sum- 
mation is extended to infinity both for x and 
y. Actual sums of the values of 9 given in 
Table I and II when added up to various x 
and y are listed in Table III. 


From the values given in Table IU, it is 
noted that while a good result can be obtain- 
ed by including only a small number of terms 
in the summation for d=0.5, many more 
terms need be included in order that a 
similarly good result may be obtained for 
d=1.0; 


Table III. 


Number of terms included | 4 
_in the summation (in both | Sum of coef 


x and y directions) | ficients 


.0107 
.1714 
-0324 
. 1045 
20332 
.0739 
.0306 
.0568 
.0277 
-0461 
.0248 
.0387 
.0220 


ee ee 


3 
4 
5 
6 
if 
8 
9 
10 
1g 
12 
13 
14 
15 


2213 
8569 
7269 
4511 
8936 
2802 
9629 
1937 
9936 
1468 
-0048 
1219 
.0036 


| 
| 


| | 


PRE ROR OrROROFRS 


3 
as 
5) 
6 
i 
8 
9 
10 
11 
12 
13 
14 
1s) 


Also, in Table IV and V, the sums of 
values in each column of the Table I and II 
are compared with the corresponding coef- 
ficients appearing in the two-dimensional 
sin a/a method, to which they should converge 
if the summation is extended to _ infinity. 
The asymptotic approach of the sums to their 
respective final values is generally good, but 
it cannot be said to be very satisfactory nu- 
merically, especially for d=1.0. This diver- 
gence may be due to some approximations 
used in the process of computing tabular 
values. 


§5. In conclusion, the writers wish to thank 
Professor C. H. Dix of the California Institu- 
te of Technology for suggesting the coopera- 
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ion of our two groups in this work. We 
lso wish to thank the California Research 
Yorporation for making it possible for us to 
ise a digital computor to compute the tables, 
ind also for permission to publish them. 
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Fault Origin Model of Earthquakes, with Special Reference 
to the Tango Earthquake, 1927. 


By 


Keichi KASAHARA 


Earthquake Research Institute, Tokyo University. 


Abstract 


Geodetic and seismological characteristics due to a fault movement have been studied 


on the basis of a simple model. 


It is intended to examine whether or not such a model 
is acceptable as a possible representation of earthquake origins. 


The model is proposed 


from the standpoint similar to that of REID’s hypothesis and the sudden occurrence of a 
fracture fault in the upper part of the earth’s crust is taken as the immediate cause of 
earthquake shocks as well as of crustal deformations. 

First, we deal with statical deformation of the earth’s crust caused by stress change 
(shear) on the fault plane, which is assumed to be developed from the surface vertically 


down to a certain depth. 


It has been found that by this model, the general decrease of 


horizontal displacement of triangulation points with distance from the fault is explained 


reasonably well. 


probable conditions at the actual earthquake fault can be surmised. 


From the comparison of the theory with the observed geodetic data, 


In the case of the 


Gémura fault (Tango earthquake), for instance, the depth (H) is estimated at 15 km, while 
the stress change on the fault plane (Y,)9 and the total amount of strain energy (Haut), 
at 3x10%c.g.s. and 4x10” ergs, respectively (length of the fault being assumed to be 


30 km). 


Discussion is also made about the characteristics of seismic waves which are generated 


by the supposed fault movement. 


The analysis based on a two-dimensional model has 


shown that the push-pull distribution of the initial P-wave is of the quadrant type, one of 


the nodal lines coinciding with the fault. 


The maximum of the spectral intensity (P- 


wave) falls on the component, whose wave-length (referred to S-wave) is approximately 


1.2 times as large as the fault’s length. 


These results also agree satisfactorily well with the observed seismological data. The de- 
duced conditions at the fault are also acceptable judging from the mechanical strength of the 


crust which is already known. 


The model adopted here affording satisfactory explanations 


of various aspects of the Tango earthquake, we might say that the physical conditions at 


the seismic origin do not differ much from those of the model. 


Possibility of applying 


the similar model to other cases is also discussed. 


1. Introduction 


It is an almost established fact that earth- 
uakes are caused by a sudden release of 
he strain energy which has been accumulated 
rithin some limited portion of the earth.” 
\Ithough the mechanism of such energy 
elease has not been known in detail yet, a 
jajority of seismologists believe that some 


1) The present paper deals only with such earth- 
uakes of large magnitude as sometimes called 
sctonic earthquakes. 


fracture-like event that occurs in the extremely 
strained portion would probably be the im- 
mediate cause of it. 

The above theory is supported most posi- 
tively by those seismologists who pay at- 
tention to the systematic distribution of the 
push-pull sense of initial earthquake motion 
recorded at many stations. P. ByEerty (1955), 
J. H. Hopeson (1957) and others have given 
the “fault-plane solutions” for a number of 
earthquakes, in which they obtained the 


16 Keichi KASAHARA 


directions of faulting assumed as the im- 
mediate cause of them. 

In fact, the solutions obtained by them for 
some destructive earthquakes are in good 
harmony with observed features of the re- 
spective seismic faults especially in regard to 
their direction of motion. The theory has, 
however, almost nothing to do with the other 
source of observational data, such as the period 
and the form of seismic waves, distribution 
of crustal deformation around the fault, and 
so on, from which much can be learned about 
the processes at earthquake origins. It is in- 
disputable that these deta provide valuable 
informations about what occurs at the seismic 
origin (f.2. H. Honpa, 1957). If we could 
examine the proposed model taking such addi- 
tional aspects of earthquakes into consid- 
eration, we would be more confident in our 
conclusion. 

It is desirable, therefore, to improve the 
model if necessary so as to make it agree 
with all the geodetic and seismological evi- 
dences of earthquakes. While it would be 
very difficult to find appropriate models for 
all the earthquakes, it is worth conducting 
an analysis for some special cases. In that 
case, more concrete knowledge about the 
physical conditions at seismic origin will be 
obtained. 


§2. Crustal deformation associated with a 
lateral fault 


Statistics shows that shallow earthquakes 
in the Japan area whose magnitude exceeds 
a certain limit are usually accompanied by 
crustal deformations (K. KasaHArA, 1957). 
Although the deformation is complicated in 
general, in several cases, such as the Tango 
earthquake of 1927, it is of a fairly simple 
feature (Fig. 1). 

It is characterized by a long fault running 
across the epicentral area as well as by the 
systematic distribution of the earth’s de- 
formation. Triangulation points on both sides 
of the fault were displaced parallel to the 
fault in opposite directions. The amplitude 
of the displacement diminishes monotonously 
with distance from the fault. If we disregard 


its detailed feature, such a tendency is likely 
to be interpreted by taking a simple mode 
of fault production. We therefore assure 
the conditions as follows (Fig. 2). 
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Fig. 1. Horizontal displacements of triangulation 
points in the Tango district (after C. TsuBo1); 
(thick lines: seismic faults). 


Fig. 2. Model of a lateral fault. 

A large part of the crust (half-space) is 
initially subject to uniform shear stress, which 
exceeds, at a certain moment, the limit of 


2 


strength of the medium. Accordingly, a 
single fracture plane takes place, which is 
assumed to be in the direction bisecting the 
angle between the two principal stress axes. We 
also assume that the fracture plane is infinitely 
long but is finite in depth. This plane being 
liberated from the stress (shear), deformation 
of the medium takes place around it due to 
the stress change. Weare able to investigate 
such a deformation analytically from the view- 
point of theory of elasticity (K. KASAHARA, 
1957). 

The solution, which satisfies the fundamental 
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equation of equilibrium as well as the boundary 
sonditions at the free surface, is given as 
‘ollows, 


u=0, 


1a: . 
v=>) (7. )isinh na—cosh nx) cos nz , 


ae 
where, uw, v, and w denote the displacements 
Nw, y, and z directions, while Fn(m=1, 2,3,---) 
ure the constants which are to be determined 
rom the following relation, 


(yo 1). Fin COS NZ . K Z ) 


Accordingly, if we know the stress distri- 
jution along the fault plane, F, and v are 
alculated. Fig. 3 illustrates the distribution 


—— xX 


ie} 05 


| (H) 


Fig. 3. Distribution of v in the vertical cross- 
section. 


f » in the vertical cross-section perpendicular 
o the fault plane, which is needed to make 
he plane free from the initial stress. We 
ee that v is the largest at =z=0 
ind decreases with x as well as with ae 
K. Sacisaka (1954), who investi- 
‘ated the crustal deformation ac- 
ompanying the North-Izu_ earth- 
uake (1930), assumed the constant 


one of the three is the most reasonable. 
In order to compare the theory with the 
observational results, we illustrate the relation 


20 30 


Fig. 4. Diminution of vp with distance (a). 


between vp and 2 in Fig. 4, wherein the depth 
of the fault, H, is taken asa parameter. The 
curves are in the tendency quite similar to 
the distribution of v) as actually observed (see 
Fig. 5). The fact that the displacement of 
triangulation points shows a characteristic 
tendency has already been pointed out by H. 
F. Reip (1911), who investigated the crustal 
deformation accompanying the California 
earthquake of 1906. However, no quantitative 
discussion of such a tendency has not been 
given yet. The present model provides a 
reasonable interpretation for the above- 
mentioned tendency and makes it possible to 
estimate the probable value of H, by applying 
the curves in Fig. 4 to the observed data. 
Fig. 5 shows that the diminution of v) can be 
well explained provided we take H=15km. 
In the foregoing discussion we dealt with the 


é (A) 


mplitude of horizontal displacement 
vith depth, whereas S. Homma (1952) 
rrived at the conclusion that v in- 
reases with z. The writer’s result 
; different from both of them. Since 
ve have no data about the features 
f deformation below the earth’s 
urface, we can not conclude which 
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Diminution of horizontal displacement (v9) with 


distance (the Tango district). 
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relative value of vp) taking v at e=z=0 asa 
unit. The absolute value of v(*#=z=0) is 
given, in the present theory, as the function 
of the rigidity of the crust, 4, and of the 
magnitude of stress change on the fracture 
plane, (Ye)ez5. Laking 2 as) >< l0M@elg-s., 
for instance, we get, (Y2)ro9=—3x10' c.g:s.. 
This is of the order of strength of 
ordinary rocks, so that the agree- 
ment of both the quantities may be 
regarded as another proof for the 
present model. Since the displace- 
ment is known from point to point, 
strain energy associated with the 
deformation can be estimated, too. 


- facts in regard to those points. 


origin. It is necessary to examine how well 
the proposed model can explain the observed 
Although it 


is difficult to deal with the dynamical charac- 
teristics of faulting strictly, an approximate 
examination for a two-dimensional case wil! 
be effected in the following. 


Taking the above-mentioned values 
into consideration, we estimate the 
energy at 4x10” ergs, where the 
integration of energy is carried out 
for the volume which is 30 km long 
along. the fault, eThisPis.can sap 
proximate value of the fault’s length, 
which is presumed on the basis of 
some aspects of the earthquake. 
Seismometrical studies of the Tango 
earthquake have concluded that the energy 
of seismic waves was as much as 1023 ergs, 
approximately. Such amount of wave energy 
can be supplied from the fault movement of 
the above-mentioned condition, provided it 
takes place in a very short time (see Section 4), 
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§ 3. Seismic waves caused by fault movement 


The wave-form and direction of the initial 
motion are regarded as another kind of in- 
formation about what occurs at the seismic 


errs 


Fig. 6. Type of the initial stress and the fracture 
(AB). 


Azimuthal distribution of displacements wu, and 


vz at a distant point caused by periodic origin force 
(As/2L=1.12). 


Let us consider that the medium is initially 
subject to uniform shear stress and that the 


fracture AB (length: 2Z) takes place in it 
abruptly (Fig. 6). Such a case can be treated 
by studying the wave equation in elliptic 
coordinates (&, 7). K. Sezawa (1928) had given 
the general solution for the equation, on the 
basis of which the writer worked out a study 
of wave generation from a fault (K. KASAHARA, 
1958). The condition that a fault appears 
abruptly may be represented by the apparent 
stress change working so as to cancel the 
initial stress on the fault plane. Since it is 
convenient to deal with the stress change in 
time with respect to its spectrum, we take 
the following equation as the boundary con- 
dition at the origin (&=0), 


EF=0,  Ey=Set, (3) 

In practice, we take several different values 
of p and carry out numerical calculation for 
the purpose of obtaining displacement due to 


P- and S-waves. Fig. 7 illustrates an example 
of azimuthal distribution of the real part of 
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u, and v; at a great distance, where mu, and 
v2 denote displacements due to P-(&-component) 
and S-wave (y-component), respectively. As 
to P-wave, it is evident that the distribution 
is of the quadrant type whose nodal lines ap- 
pear along and perpendicular to the fault. 
Such a tendency being noticed for all the 
cases, we may conclude that the model is 
consistent with the observed pattern of push- 
pull distribution as shown in Fig. 8. 


Fig. 8. Azimuthal effect of the push-pull distri- 
bution (initial P-wave) observed in the case of 
the Tango earthquake (after S. KUNITOMI). 


Amplitude 
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Fig. 9. Amplitude of wu plotted against A;/2D. 


The result of numerical calculation for vari- 
ous values of p enables us to deduce the 
spectrum of the wave which is to be gene- 
rated from the fault. In Fig. 9 are shown 
the amplitudes of wm (relative values) for dif- 
ferent periods of the origin force, wherein 
the abscissa is given with respect to the ratio 
of As (wave-length referred to the velocity of 
S-wave) to 2L. Taking the spectrum of an 
impulse into consideration, it is readily seen 
that Fig. 9 also indicates the spectrum of the 
wave caused by an impulsive excitation at the 
origin. Accordingly we know that the maxi- 
mum intensity of the spectrum of P-wave 


falls on the component whose wave-length, 
A;, is about 1.2 times as large as 2L. That 
is to say, the period of the maximum intensity 
(we called it as the characteristic period, 7) 
becomes as long as 11 sec, provided we take 
2L and v; as 30 km and 3.2 km/sec, respectively. 
It is more natural to consider that the stress 
change associated with faulting is step-like 
rather than impulsive. If we take this supposi- 
tion, the characteristic period will be estimated 
at a larger value, 20 sec or more, say. Seismo- 
grams of the Tango earthquake recorded at 
Hongo, Tokyo, indicates 23 sec as the charac- 
teristic period, so that the present model does 
not contradict with this fact. 


Synthetical examination of the model 
and discussion of the related problems 


§ 4. 


The foregoing two sections have dealt with 
geodetic and seismological characteristics of 
the fault model, separately. In the present 
section, therefore, we shall make synthetical 
examination taking all the sorts of evidence 
of the Tango earthquake into consideration. 
Discussion that has been made previously 
indicates that the probable condition of the 
faulting is as shown in the central part of 
Fig. 10 (model A), where, 2L, H, 6, and 9 
denote, respectively, length, depth, strike, 
and dip of the fault plane, whereas (Yz)po is 
the stress change (shear) to which the plane 
was subject. Although it is desirable to make 
all the discussion referring to model A, we 
shall take more simplified models, B and C, 
to avoid difficulties in mathematical treatment 
of the original one. Model B is derived from A 
by taking ZL as infinity while all the other 
conditions are kept unchanged. Model C is 
given in the similar way by taking H as in- 
finity. 

Items which are arranged on the left of 
the proposed models in Fig. 10 are geodetic 
aspects of the earthquake. The previous dis- 
cussion (see Section 2) has proved that the 
diminution of vy) can be well interpreted from 
the standpoint of model B. Geodetic obser- 
vations revealed, also, more complicated 
features of deformation distribution at the 
tails of the faults. Characteristics of these 
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Fig. 10. Relation between the proposed model and various sorts of aspects of the Tango earth- 


quake. 


secondary effects can be seen more clearly 
with respect to distribution of various sorts 
of strain components, such as dilatation, shear, 
and so on (C. Tsusor, 1933). As it was hardly 
possible to deal with such effects mathematical- 
ly, we carried out model experiments previ- 
ously, by which it was shown that some part 
of the above-mentioned features could be ex- 
plained from the viewpoint of model A (K. 
KASAHARA, 1957). 

It is a well-known fact that the mechanical 
strength of the crust is not infinitely large 
but there is a certain limit of strain energy 
to be accumulated there. This information 
has been got from the distribution of the 
maximum shear around seismic faults as well 
as from the upper limit of the magnitude of 
great earthquakes (C. Tsusor, 1956). It is 
notable that the value of (Y.)o as given in 
model A does not contradict with this fact. 

We shall now turn the point of discussion 
to seismological aspects of the earthquake, 
which are arranged in the right half of Fig. 


Items surrounded by broken lines are the assumptions taken in the discussion. 


10. The energy of seismic waves of the 
Tango earthquake has been estimated at about 
10° ergs. We know, after rough estimation 
based on model B, that the strain energy as- 
sociated with the crustal deformation amounts 
to 4x10" ergs. That is to say, the energy 
which is as much as 4x10 ergs is to be 
radiated as elastic wave, provided such a 
deformation arises instantly. The energy as- 
sociated with the vertical component of defor- 
mation and with other strain components is 
excluded from the present estimation, so that 
we shall find better agreement of both kinds 
of energy if they are taken into account. 
Seismometrical studies of the Tango earth- 
quake have found out that the epicentre was 
located on the fault. However, lack of 
knowledge about the depth of the fault plane 
has left uncertainty in our discussion of 
spatial relation of the seismic origin with the 
faults. The present study provides a concrete 
answer for the question that the hypocentre 
(focal depth: 14km) located at the lower part 
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of the fault plane. We can not explain what 
this fact means, but we may suppose, as one 
of the possible explanations, that the first 
breaking of the fault arose at that point. 
Spectrum and azimuthal effect of the initial 
P-wave have been investigated in the former 
Section. Seismograms recorded at Tokyo 
(Hongo) indicates the characteristic period, 7), 
to be 23 sec. Mathematical study of model 
C proved that 7) is about 11 sec in the case 
of impulsive excitation at the origin while 20 
sec or more, for the step-like one. We have no 
concrete basis of presuming the mode of 
stress change, but we might imagine, ap- 
proximately, that the faulting causes step-like 
stress change. The azimuthal effect which is 
derived from the model is also consistent with 
the observed result (see Fig. 8). Thus we may 
conclude that the proposed model is acceptable 
with respect to spectrum and azimuthal effect, 
provided it is permitted to apply the two- 
dimensional model to the actual case. 
Results of the foregoing discussion appear, 
therefore, to justify the proposed conditions 
for the Tango earthquake. We have toremark, 
however, that no discussions have been made 
about some other factors notwithstanding their 
importance. These are shown in Fig. 10 
being surrounded with broken lines. First, we 
did not discuss the physical meaning of the 
vertical component of deformation. It is 
needless to say that the factor brings, also, 
valuable information about the nature of the 
rigin, so that further improvement of the 
fault model should be made so as to take the 
mentioned factor into consideration. Second, 
we assumed that the faulting took place in a 
very short time. Rapidity of fault production 
yas much influence on the efficiency of energy 
-adiation as well as on the feature of spectrum. 
Although we have no observational data ac- 
urate enough for discussing the rapidity, some 
macroseismic studies of other earthquakes 
ire likely to suggest that main part of the 
‘rustal deformation is completed within a 
short time, some ten seconds, say. If it is so, 
he above-mentioned conclusions would not 
uffer any drastic alteration. The third as- 
umption was taken on the possibility of strain 


accumulation. Geophysicists have discussed 
the problem in relation to the tectonophysical 
conditions of the earth (f.7. H. Jerrreys, 1952; 
T. Martuzawa, 1953), who found out several 
possible causes of strain accumulation. Al-: 
though the problem has not been solved com- 
pletely yet, we may say that we have good 
reason to justify the above-mentioned as- 
sumption. The writer worked out a study of 
Martuzawa’s model for the purpose of examin- 
ing whether or not the strain accumulation 
is much influenced by viscous flow of rocks 
(K. Kasanara, 1956). According to the study, 
energy dissipation due to the viscous flow is 
not so large as to deny the possibility of strain 
accumulation. 

Thus it has been shown that the proposed 
conditions of faulting is acceptable as a pos- 
sible model of the origin of the Tango earth- 
quake. The question that is to be raised on 
the next step is if the similar model is ap- 
plicable to the other cases. Unfortunately, 
we are far away from getting a definite 
answer for the question. However, we know 
some special cases, such as the North-Izu 
earthquake (1930), whose aspects (geodetic and 
seismological) were similar to those of the 
Tango earthquake. It is very likely that the 
present model is also applicable to them 
without great alteration of its fundamental 
conditions. 

On the contrary, it would be hardly possible 
to find out concrete conditions at origins for 
the other earthquakes which are not ac- 
companied by crustal deformation of simple 
features. The writer would like here only to 
mention the result of his previous analysis, 
which treated the density of energy accumu- 
lation around seismic origins and proved that 
there was no significant difference, from such 
a standpoint, between earthquakes occurring 
at different depth in the earth (K. Kasanara, 
1957). This fact might be the evidence for 
their mechanism of occurrence being similar 
with each other. 


§5. Conclusions and acknowledgement 


The writer worked out a study of geodetic 
and seismological characteristics of a modei 
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of faulting with the intention of finding out 
the probable physical conditions at seismic 
origins. It is assumed that a vertical fracture 
plane of finite dimension takes place in the 
extremely strained crust (half-space), which 
results in radiation of seismic waves as well 
as in crustal deformation around the fault. 

We applied the results of the writer’s 
previous studies to the case of the Tango 
earthquake, and determined, numerically, 
various conditions at the seismic origin as 
shown in Fig. 10. Synthetical examination of 
the proposed conditions with the observational 
data has proved that the model adopted here 
is consistent not only with the geodetic aspects 
but also with seismological ones. Taking this 
result into consideration, we might conclude 
that the physical conditions at the actual origin 
do not differ much from those of the model. 

We can not say much about further ap- 
plication of the model to other cases. 
Some of the earthquakes, such as the North- 
Izu earthquake, are likely to be of the similar 
features with those of the Tango earthquake. 
In such cases, there is a possibility of inter- 
preting the mechanism at their origins from 
the present standpoint. 

The writer expresses his sincere thanks to 
Prof. Chuji Tsusor and Prof. Takahiro Hacr- 
wara of Tokyo University for their encour- 
agement and support throughout the prepa- 
ration of this paper. Also the writer owes 
many thanks to Dr. Tsuneji RIKITAKE who 
gave him valuable advice. 

Notes added in proof. Crustal deformation 
around a fault has also been investigated by 
L. Knoporr (1958, Geophys. Journ., R.A.S., 1, 
44-52) and P. Bygeriy and J. DENoyeEr (1958, 
Contributions in Geophysics in Honor of Beno 
GUTENBERG, 17-35), recently. They dealt with 
che diminution of vp with distance and arrived 
at a conclusion which does not differ much 
from the writer’s. 


References 


lennon, 1B % 
1955 Nature of faulting as deduced from seismo- 


grams. Geol. Soc. Amer., Special Paper 
62, 75-85. 
Hopaeson, J. H. 

1957 Nature of faulting in large earthquakes. 

Bull. Geol. Soc. Amer., 68, 611-644. 
HOMMA, S.: 

1952 Energy of the crustal deformation. Quart. 
Journ. Seismology, C.M.O., 16, 57-66 (in 
Japanese). 

HONDA, H.: 
1957 The mechanism of the earthquakes. Sci. 
Rep. Tohoku Univ., Ser. 5, 9, Supplement. 

JEFFREYS, H.: 

1952 “The Earth” 3rd ed. (Cambridge). 
KASAHARA, K.: 

1956 Strain energy in the visco-elastic crust. 
Bull. Earthq. Res. Inst., 34, 157-165. 
The nature of seismic origins as inferred 
from seismological and geodetic observations 
(1). Bull. Earthq. Res. Inst., 35, 473-532. 
The nature of seismic origins as inferred 
from seismological and geodetic observations 
(2). Bull. Earthq. Res. Inst., 36, 21-53. 
KUNITOMI, S.: 

1930 Note on the Tango earthquake of March 
7, 1927. Geophys. Mag., 2, 65-89. 

MATUZAWA, T.: 

1953 Feldtheorie der Erdbeben. 

Res. Inst., 31, 179-201. 


Loaf 


1958 


Bull. Earthq. 


REID, H. F.: 
1911 The  elastic-rebound theory of  earth- 
quakes. Bull. Dept. Geol. Sci., Univ. 


California, 6, 413-444. 
SAGISAKA, K.: 

1954 On the energy of earthquakes. Geophys. 

Mag., 26, 53-82. 
SEZAWA, K.: 

1928 Propagation of elastic waves from an 
elliptic or a_ spheroidal Bull. 
Earthq. Res. Inst., 2, 29-48. 

TSUBOI, C.: 

1933 Investigation on the deformation of the 
earth’s crust found by precise geodetic 
means. Jap. Journ. Astro. Geophys., 10, 
93-248. 

1956 Earthquake energy, earthquake volume, 
aftershock area, and strength of the earth’s 
crust. Journ. Phys. Earth, 4, 63-66. 


origin. 


JOURNAL OF PHYSICS OF THE EARTH, VOL. 6, No. 1, 1958 23 


Observation of Near-by Microearthquakes with Ultra 


Sensitive Seismometers at Matsushiro, Japan. 


By 
Toshi Asapa 
Geophysical Institute, Tokyo University 


Shigeji SuyEHIRO 


Japan Meteorological Agency 


Kei AKAMATU 


Geophysical Institute, Tokyo University 


Abstract 


A temporary observation of near-by microearthquakes was made with Ultra Sensitive 


Seismometers at Matsushiro. 


The purpose of this experiment was to see the relation between 


the frequency of near-by microearthquakes and the annual number of near-by shocks felt 


at Matsushiro. 


The results and conclusions are as follows. 
The average annual number of earthquakes felt at Matsushiro is about 
As M of the microearthquakes are from less than —1 to 0, and M of near-by felt 


occurs every hour. 
ee 


An average of one microearthquake 


earthquakes are about 4~5, we have an equation, 
log N=3.6—0.77M , 
where WN is the annual number of earthquakes whose M are between M-0.5 and M+0.5. 
This equation holds good in the relation between N and M of earthquakes with M from 
—1 to 5 occurring in the vicinity of Matsushiro (area=60 x 60 km?). 
As the equation is valid in such a wide range of M, it will be possible to determine 
the seismicity of a certain region within a rather short period of time through observing 


microearthquakes. 


The # of the minimal recorded earthquake was about 109 ergs in this experiment. Its 


M was approximately —2. 


Introduction 


$1. 

One of the present writers made a tempo- 
rary observation of near-by microearthquakes 
at Mt. Tsukuba with Ultra Sensitive Seismo- 
meters in 1955 (Asapa, 1957a). The principal 
purpose of this experiment was to estimate 
the energy of the smallest detectable shock 
and to see if the equation between the number 
and magnitude, log N=a-bM, holds good in 
a wider range of M than has ever been in- 
vestigated. 

The conclusions of the experiment at Mt. 
Tsukuba are as follows. 
1) The energy of the minimal recorded shock 
was of the order of 10 ergs. This value 
was calculated according to the formula of 


the kinetic energy of spherical waves, assum- 
ing a proper constant of absorption. Also the 
M of this shock was determined to be —1.1. 
2) In the vicinity of Mt. Tsukuba an average 
of 200 shocks of energy (10!°~10'* ergs) occur 
every 24 hours. Shocks of M 4~5 occur 40 
to 50 times every year in the same area. 
Accordingly, we can conclude safely that the 
equation, log N=a-bM, (b=+-0.8), is valid in the 
range of M from less than 0 up to 5. 

3) The space distribution of the shocks of 
M Z0 agrees with that of the shocks of M 
4~5. 

4) The predominant frequencies in the mo- 
tions due to near-by microearthquakes are 
about several tens of a cycle per second, Ac- 
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cordingly, these shocks can be recorded only 
by a seismometer sensitive enough in the 
frequency band from 10 cps to 100 cps. 

In order to check the above conclusions the 
present writers installed the Ultra Sensitive 
Seismometers at the Matsushiro Observatory 
of Japan Meteorological Agency and made a 
temporary observation in August 1956. The 
location of Matsushiro is shown in the map in 


Fig. 1. 
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Fig. 2. Frequency response of ten-cps velocity 
seismometer. 
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Two ten-cps velocity seismometers, vertica 
and horizontal, and one one-second horizonta 
velocity seismometers were set up in a vault 
(a), and a one-second horizontal velocity seis- 
mometer was set up at a_position-(b) about. 
200 meters apart from the position-(a) in order 
to discriminate small earthquakes from dis- 
turbing noises:(AsApDA, 1957a). | 

Their electrical outputs were amplified and 
recorded side by side on a photograhic paper 
with an oscillograph. The out put of a ten-cps 
vertical seismometer was amplified and differ- 
enciated in another channel simultaneously and 
recorded on the same photographic paper. A 
visible recorder was also used. Its galvano- 
meter was fed by a power amplifier connected 
to the preamplifier of one-second seismometer 
(No. 3, in Table 1). The visible recorder was 
operated continuously for 13 days. The oscillo- 
graph was operated only when the ground was 
very quiet. 

Frequency responses of the ten-cps and one- 
second velocity seismometers are given in Fig. 
2 and Fig. 3. The sensitivity of the ten-cps 


percent response 


cm/sec mm 


— 
100- 


Fig. 3. Frequency response of one-second velo- 
city seismometer. 


Table 1. Natural periods, components and positions of seismometers. 
re ne 
No. i" T Component Position Trace in photographic paper 

No. 1 10 cps Vertical a i A 
Be] 
No. 2 10 cps Horizontal a . Cc 
No. 3 liesec Horizontal a : D ’ 
. - 4 (Power amplifier—Smoked paper Merit, 
No. 4 1 sec Horizontal b Gi - : ;: } 


E 


* The output of the No. 1 ten-cps seismometer was amplified and differentiated in another 
channel and recorded in Trace B on the same paper. 
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seismometer is 8.6 10-?cm/sec/mm or 1.2x 
10° mm/cm/sec in the frequency band from 10 


Table 2. Frequency distribution of trace 
amplitude of seismograms recorded with 
one-second velocity seismometer with visi- 
ble recorder. 


cps to 200 cps, The sensitivity of the one- 
second seismometer is about the same order 
as that of the ten-cps seismometer in the fre- 
quency band from 2 cps to 20 cps. However, 
the former is not appropriate for recording the 
shocks of M less than 0. The frequency 


Price Amplitude Piegieney responses of Galitzin and other seismometers 
P for routine observation compared with those 
seek Sie a of Ultra Sensitive Seismometers are given in 
S11 4.9 22 | 86 Fig. 4. 
pine “ §2. Minimal recorded earthquake and its 
Glen.) 6 ea 
7T1— 8.0 11 One hundred and twenty one shocks of hypo- 
Siar a 29 central distances smaller than 70~80 km were 
10,1—11.0 2 recorded by the visible recorder during the 
1151 —12.0 4 course of 13 days. Forty-three shocks occur- 
. ere : : a ring within the distance of 70~80 km from the 
14.1-15.9 1 observatory were recorded on the oscillograph 
15.1—16.0 1 during 45 hours in total, or an average of 
ton i7. 0 1 about 24 shocks occurred every 24 hours. 
eer we : 6 The frequency distribution of the maximum 
19.1—20.0 iA trace amplitude of the seismograms on the 
20. 1—21.0 » visible recorder of the one-second 
velocity seismometer is shown in 
| Table 2 and Fig. 5. This distribution 
: Ultra Sensitive can satisfy the equation (SuzukKI, 
Magnification SeiSmometer Pe 10, cps 1953), 

A N(A)-dA=kA-"dA, (2.1) 
ite) where the value of m is approxima- 
10° Frequency 

To=10S 
10° 
Te 1.08 
10° Te= 3-05 
Galitzin 
10° 


yeoKerey Keys Is 


Period 


1Ocps !lOOcps 


iB 10. 
Trace Amplitude ( mm) 


100. 


Fig. 5. Frequency distribution dia- 
gram of trace amplitudes of seismo- 


Fig. 4. Frequency responses of Ultra Sensitive seismo- 
meters, SUYEHIRO vertical seismometer (T';=1.0s, Te—o.0s) 
and Galitzin seismometer. 


grams recorded by one-second velocity 
seismometer with visible recorder. 
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swe ergs (Asapa, 1957a). 


E one-sec H 


D one-sec H 


Bten-cpsV_ § 
(Accel.) 


A ten-cps V 


Fig. 6. Seismograms of minimal recorded earth- 
quake recorded at Matsushiro. 


tely 1.8. The E of these shocks are from 10” 
ergs to 10” ergs as will be described below. 

In Fig. 6 the seismograms of one of the mini- 
mal recorded earthquakes are given. These 
seismograms give us the information on the 
predominant periods because the seismometers 
of different responses are used. Comparing 
the seismograms A, C, and B with D and E, 
one can easily see that their predominant fre- 
quency is higher than several tens of a cycle 
per second. 

The P-S interval of this shocks is 0.96 sec, 
and the maximum velocity amplitude is about 
5x8.6x10-7cm/sec, or 4x10-®cm/sec. The 
minimal detectable earthquake recorded at Mt. 
Tsukuba in 1955 had a P-S interval of 4.0 sec 
and a maximum velocity amplitude of 3.5x 
10-® cm/sec. Its E was estimated at about 101° 
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6 if the wave attenuation in the crust is as- 
sumed. Now assuming that the attenuation 
is proportional to 2nd or 1.5th power of the 


4 distance (Den, 1953) and E is proportional to 
~ the square of the amplitude, E of this minimal 


recorded shock is 

10" x (4/3.5)? x (0.96/4)2=6.2 x 108 ergs, or 

101° x (4/3.5)? x (0.96/4)?-®=1.1x 10° ergs. 
In determining E of a shock, the accuracy 
being very low, we may say that EF of this 
shock is about of the order of 10° ergs. The 
approximate value of its M can be obtained 
also, if we assume the decay of waves is 
proportional to 1.5th power of the hypocentral 
distance. The maximum trace amplitude of 
this shock should be about 0.012 microns if 
it were recorded by a “ standard ” seismometer 
at the hypocentral distance of 100km as its 
maximum displacement at the hypocentral dis- 
tance of about 7km is 2.3x10-8cm at 30 cps 
or 4x10-®cm/sec. Accordingly, the values of 
its M is about —2. 

As the sensitivity of the seismometer with 
power amplifier and visible recorder decreases 
in the frequency range higher than about 20 
cps, E of the shocks (P-S< 10 sec) recorded 
on the smoked paper are from 10” ergs to 
10”? ergs. 


$3. Space distribution of microearthquakes 
and small shallow shocks felt at Ma- 
tsushiro 


The frequency distribution of P-S intervals 
of microearthquakes occurring in the vicinity 
of Matsushiro is given in Table 3 and Fig. 7. 
The frequency distribution of P-S intervals 
of the microearthquakes recorded at Mt. 
Tsukuba is given in Fig. 8. The diagram in 
Fig. 8 shows that most of the shocks have 
P-S intervals larger than 5 sec, although the 
greater part of the shocks recorded at Matsu- 
shiro have P-S intervals smaller than 5 sec. 
From this we interprete that the earthquakes 
occur not in the crust but in the upper part 
of the mantle in the vicinity of Mt. Tsukuba, 
while in the vicinity of Matsushiro the shocks 
occur in the crust only (Asapa and Den, 1954). 


Accordingly, we may © 
- determine the E of the shock shown in Fig. 
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The list of the shocks felt at Matsushiro is 
given in Table 4. Comparing the frequency 
diagram in Fig. 7 with the data in Table 4, 
it may be seen that the microearthquakes 
occur almost in the same part of the crust as 
the shocks of M 4~5 listed in the table. 

The list in Table 4 shows that the number 
of near-by earthquakes felt at Matsushiro is 
an average of 1.2 every one year. As the 
magnitudes of microearthquakes occurring in 
the vicinity of Matsushiro (area~60 x 60km?) 
are approximately from —1 to 0 and the an- 
nual frequency is about 8700 (24 eqks/day), we 
can calculate the values of a and 0 in the equa- 
tion, log N=a-bM (GuTEnBERG and RIcHTER, 
1942, 1956), for the range of M from —1 to 
5. Putting into the equation the number of 
shocks of M 4~5, we get, 


G—o.0, and) 0—OMir. 
Of course the second digits of these values 


are of little significance. 
The equation 


log N=3.6—0.77M (3.1) 


holds good in the range of M from —0.5 to 
4.5 where W is the annual number of earth- 
quakes whose M is between M-0.5 and M+ 
0.5. As the equation is valid in such a wide 
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Fig. 7. Frequency distribution of P-S intervals 
of near-by microearthquakes recorded at Matsu- 
shiro. 
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Fig. 8. Frequency distribution of P-S intervals of 
near-by microearthquakes recorded at Tsukuba. 


Table 3. Frequency distribution of P-S intervals of shocks recorded at Matsushiro 
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Table 4. Minor shocks felt at Matsushiro and its visinity 
a nn ne en] DRENDEN 


d. h. m. (J.8.T.)| C.M.O. Intensity | £8 
1952 July 17 07 32 | 0 
Dec 27 14 59 ll 1.4 | [36.5N, 138.2D] 
Feb. 1 15 47 0 2.3 |felt at Yashiro 
1953 May 7 01 52 0 3.0 | felt at Nagano 
July 6 03 35 0 Hill 
July 6 03 38 1 les 
Feb. ib 14 44 1 4.1 |[36.7N, 137.8E] 
ane Sept 10 04 56 0 felt ot Nagano 
Aug. I 12 50 0 felt at Nagano. 
1955 Oct. 6 07 21 0 2.1 | felt at Nagano 
Oct. 24 03 40 0 ies 
Jan 27 15 04 i ? 
1956 
May 23 13 7A 0 4.6 


range of M, it is possible to determine the 
seismicity of a certain region within a rather 
short period of time through observing micro- 
earthquakes. 

This was the case in the “shallow ” earth- 
quakes occurring in the vicinity of Mt. Tsu- 
kuba in the depth from 30 to 70km, although 
the value of coefficient 6 is slightly different 
from the present one. 

The frequencies of bigger earthquakes oc- 


Table 5. Number of shocks calculated ac- 
cording to the equation, (3.1). 


M Frequency 
4~5 1.2/1 year 
5~6 2/10 year 
6~7 3.4/100 years 
7~8 5. 7/1000 year 


Table 6. Earthquakes (MV=>6.0) occurred 


curring in this area calculated according to this 
equation (3.1), are shown in Table 5. The ma- 
jor earthquakes which have actually occurred 
in this region are listed in Table 6. The 
calculated values seem to agree with the actual 
ones. The difference is of course statistically 
insignificant. The size of the data is, however, 
too small to conclude decisively that the equa- 
tion (3.1) holds good in the magnitude range 
from —0.5 to 7.5. 


§4. Seismograms 


The seismograms in Fig. 6 show very clearly 
that their predominant frequencies are more 
than several tens of acps. It can be seen that 
the one-second velocity seismometer cannot 
register the minimal recorded earthquake well. 
Most of electro-magnetic seismometers coupled 
directly with galvanometers have their maxi- 


in the vicinity of Nagano and Matsushiro. 


_—_— SSS 


1953 October 13 


Nojiriko Wse20k 36.8°N M=6.0 
1941 January 15 Nagano 138 .2°E 36.7°N M=6.4 
1897 January 17 Chikuma-gawa M=6.3 
1858 April 23 Shinano-Matsushiro M=5.9 
1853 January 26 Shinano-Hanishinagun M=5.9 
1847 May 8 Zenkoji M=7.4 
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Fig. 9. Seismograms of shocks whose hypocentral distances are smaller than 180 km. 
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Records of distant earthquakes. 
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mum magnification at 0.2 sec or longer periods 
and their magnification mostly 10° or smaller. 
Accordingly, the earthqukes with E (10" ergs 
or smaller) cannot be registered well by these 
seismometers even if the hypocentral distances 
are 10 to 20 km. 

The present writers have obtained many 
seismograms of near-by microearthquakes in 
this experiment. These have enabled us to 
see qualitatively the relation between the pre- 
dominant frequency and the magnitude of micro- 
earthquakes (Asapa, 1953) (Ax1 1955). Several 
examples of near-by earthquakes are given 
in Fig. 9. In comparing these seismograms 
recorded by the seismometers with different 
responses, one can see that the waves of 
several tens of acps are predominant. Only the 
earthquakes of FE (10" ergs or bigger) (e.g., 
(d), (e), and (f) in Fig. 9) have the wave trains 
with periods of about 0.5 sec. 

The seismograms in Fig. 10 are those of 
shocks with hypocentral distances from about 
200km to 400km. The waves of these shocks 
are not recorded by ten-cps velocity seismo- 
meters. 

In Fig. 11 the examples of seismograms of 
distant earthquakes are given. Fig. 12 shows 
the seismograms of the same earthquakes re- 
corded by a seismometer of 7)=1 sec and 
T,=3 sec whose response is shown in Fig. 4. 

The seismogram A in Fig. 12 is due to the 
same earthquake as No. 3 in Fig. 10. In Fig. 
12, B is the same as No. 2 in Fig. 10, and C 
in Fig. 12 is the same as No. 1 in Fig. 11. The 
epicenter of the earthquake, No. 1 in Fig. 10, 
is in Bolivia, 15°S 68°W, h=100km M=63~ 
64). The earthquakes shown in Fig. 10 occur 
red in Molucca Passage (M5). 

The earthquakes shown in Fig. 9 were not 
recorded by the seismometer of 7)>=1 sec and 
T,=3 sec because of its low sensitivity. 

It should be stated also that the seismo- 
grams of the earthquakes with hypocentral 
distances from several hundred to one thou- 
sand km recorded at Mt. Tsukuba contain the 
waves with frequencies shorter than 10 to 
20cps. This stands in a striking contrast 


with the case of; Matsushiro and may be 
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because the attenuation of wave in the crust 
or in the upper part of the mantle under 
Mt. Tsukuba is smaller than that under Ma- 
tsushiro. 


§5. Conclusion 


The results of the present experiment made 
at Matsushiro lead to the following conclu- 
sions. 

(1) Microearthquakes of E from 10” to 10” 
ergs and bigger shocks of magnitudes up to 
4 or 5 occur in the shallow part of the crust 


ial. 
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Seismograms recorded by a SUYEHIRO vertical seismometer. 
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No. 2, C—Fig. 11, No. 1. 
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in the vicinity of Matsushiro. No earthquake 
occurs under the Mohorovici¢ discontinuity 
here. This is the case in the greater part of 
other areas in the world although almost all 
earthquakes occur in the upper part of the 
mantle in the vicinity of Mt. Tsukuba. 

(2) The energy of the minimal detectable 
earthquakes recorded at Matsushiro is about 
10° ergs. 

(3) In the vicinity of Matsushiro an average 
of 24 microearthquakes occur every 24 hours. 
Six earthquakes of M 4~5 occurred during 
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these 5 years in the same area (area is ap- 
proximately 60x60 kmz?). 
Accordingly, the equation, 


log N=3.6—0.77M, 


holds good in the range of M from —0.5 to 
4.5, where N is the annual number of earth- 
quakes whose magnitude is between M—0.5 
and M-+-0.5 occurring in this area. 

The seismicity of this area is about 1/5~ 
1/10 of that in the vicinity of Mt. Tsukuba. 
(4) The problem still is to see whether or 
not that the extrapolated values of N in the 
equation written above for larger values of M 
coincide with the actual values of frequency 
of earthquakes of M from 4.5 to 7.5. No 
statistically significant difference between the 
calculated and actual values is found but the 
size of the present data is small and insufficient 
to conclude that the range of M in which the 
equation holds good can cover from —0.5 up 
LOulno% 


Aknowledgment 


The authors wish to express their thanks 
to Drs. WapatiI and SacisaKA who permitted 
them to carry their experiment at The Ma- 
tsushiro Observatory of Japan Meteorological 
Agency. The authors also thank the members 
of the Observatory for their kind help and 
Messrs. Kurimoro and TerasHima for their 
yenerous assistance. 


Bibliography 
AKI, K., 

1955 Correlogram Analysis of Seismograms. Zi- 
sin If (Journal of Seism. Soc. Japan), 8 
pp 99-107. 

ASADA, T., 

1953 On the Relation between the Predominant 
Period and the Maximum Amplitude of 
Earthquake Motions. Zisin II (Journal of 
Seism. Soc. Japan), 6, pp 69-73. 

1957 a. Observation of Near-by Earthquakes 

with Ultra Sensitive Seismometers. 
Journal of Physics of the Earth, 5, pp 
83-113. 

b. Frequency Distribution of Earthquake 
Magnitude, Seismicity and Related Pro- 
blems. Zisin II (Journal of Seism. Soc. 
Japan), 10, pp 24-34. 

ASADA, T. and DEN, N., 

1954 On the Frequency Distribution of the P-S 
Intervals of Earthquakes Recorded at a 
Certain Station. Zisin II (Journal of Seism. 
Soc. Japan), 7, pp 37-44. 

DEN, N., 

1953 On Very Low Frequency Amplifiers to Ob- 
serve the Explosion Seismic Waves. Zisin 
II (Journal of Seism. Soc. Japan), 6, pp 
101-108. 

GUTENBERG, B. and RICHTER, C. F., 

1951 Earthquake Magnitude, Intensity, Energy 
and Acceleration B.S.S.A., 32, pp 163-191. 

1956 Earthquake Magnitude, Intensity, Energy 
and Acceleration B.S.S.A., 46, pp 105-146. 

SUVA GIZA, JA 

1953 A Statistical Study on the Occurrence of 
Small Earthquakes, I, Sci. Rep. Tohoku 
Univ., Ser. 5, Geophys., 6, pp 177-182. 


? 


me) ae “ 


a 


a 


JOURNAL OF PHYSICS OF THE EARTH, VOL. 6, No. 1, 1958 


Earthquake Province—Domain of Sympathetic Seismic Activities 


By 
Chuji Tsusor 
Geophysical Institute, Tokyo University, Tokyo. 


Abstract 


The notion of earthquake province has been introduced. It is a domain within which 
various portions show sympathetic seismic activities. It may be interpreted to be repre- 
senting a stationary site over which one seismic field continues to exist, although its inten- 
sity may vary according to time. By means of correlation studies of earthquake numbers 
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in various compartments, six such earthquake provinces have been established in Japan. 


§1. The term “earthquake zone” or “earth- 
quake belt” has long been used in seismo- 
logical literatures for describing the charac- 
teristic pattern of distribution of earthquake 
epicentres on the earth’s surface. As is well- 
known, epicentres are not evenly distributed 
all over the earth’s surface, but are signifi- 
cantly crowded in several well-defined long 
strips of area which really deserve the name 
of zone or belt. Either of the terms is un- 
doubtedly appropriate and useful for repre- 
senting the global picture of epicentre distri- 
bution. The Circum-Pacific Zone is one of the 
well-known examples. But going into a little 
further detail, one will easily find that the 
belt or zone is not one continuous figure at all, 
but is actually an array of clusters of epicentres 
arranged at irregular longitudinal intervals. 
As an illustration, the epicentres of 3,147 
important earthquakes (felt beyond 200 km) 
which took place from 1900 through 1950 in 
and near Japan are plotted in Fig. 1. Even 
in this Japanese area, which is noted for its 
extremely high seismic activity and which, as 
a whole, is regarded to be a representative 
portion of the Circum-Pacific Earthquake Zone, 
the number of earthquake epicentres is seen 
to vary considerably along the island arc. At 
several places, the epicentres are clustered to 
form more or less definite groups, while there 
are several other places from which earth- 
quakes occur rarely. In a case like this, the 
favorite term earthquake zone or belt will 
easily cease to have its original significance. 


$2. In order to make this uneveness in the 
distribution of epicentres in and near Japan ap- 
pear more clearly, the following simple 
method of representation has been attempted. 

Japanese islands as a whole are of a form 
of approximately circular arc, the centre being 
located near Vladivostok. From the approxi- 
mate centre of this circle, equally spaced 128 
radiating lines were drawn by means of which 
the whole plane was divided into as many 
sectors. The sectors are called (1), (2),---, 
(128), the numbering starting from the one 
that lies to the due east of the centre and 
proceeding in clockwise direction (Fig. 2). 
These sectors into which Japanese islands are 
divided will be called “ compartments” in this 
paper. The compartments are all approxi- 
mately 40km long as measured longitudinally 
along the island arc. The compartments that 
cover the area under study are 44 in number, 
namely, (122), (123), 2°; (28), 1). (2), es, 
(37). 

In Table I are given the number of earth- 
quake epicentres in each of the compartments. 
The earthquakes taken in this statistics are 
those which took place from 1912 through 1955 
and which were reported as felt beyond 200 
km from the respective epicentres. Here- 
after, these earthquakes will simply be re- 
ferred to as earthquakes. In column A of 
the table are given the total numbers of 
earthquakes for the whole length of the said 
period of time, while in column B and C; the 
numbers for its first half (1912-1933) and for 
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its latter half (1934-1955) separately. The 
values given in Table I are shown graphically 
by the curves in Fig. 3. 

Two things of importance are noticed in 
Fig. 3. First, the curve A is not a flat and 


Epicentres of 3,147 earthquakes in and near Japan. (Felt beyond 200 km. 1900-1950). 


smooth one at all and its several maxima cor- 
respond to epicentre groups. Secondly, the 
curve B and C which represent the numbers 
of earthquakes during two different time in- 
tervals are similar in shape to each other in- 
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Fig. 2. Division of Japanese islands into compart- 
ments. 


dicating that the biased spatial distribution of 
epicentres along the island arc is a time- 
persistent feature. So, if we investigate this 
distribution, it is unlikely that we are dealing 
with phenomena of accidental or sporadic na- 
ture. On the contrary, from this very distri- 
bution, something significant about the nature 
of earthquake occurrence is likely to be 
learned 


§3. In the foregoing sections, the word “ epi- 
centre group” has vaguely been used. Many 
authors have also used this term. But the 
question in what way a great many of epi- 
centres can significantly divided into such 
groups will not be easy to answer. No clear 
criterion for this purpose appears to have been 
given. 

It may be possible to establish groups of 


Table I. Number of earthquake epicentres 
(perceptible beyond 200 km 


from epicentres) 


A 1912-1955 
B_ 1912-1933 
C 1934-1955. 
ee eee eee 

Compartments | A B | (e 
22 12 4 8 
3 31 8 23 
4 40 15 25 
5 91 48 | 43 
6 124 60 | 64 
if 125 69 | 56 
128 100 69 31 
1 62 28 34 
2 56 | 37 19 
3 74 | 44 30 
4 94 54 40 
5 61 49 12 
6 58 38 20 
7 53 26 27 
8 61 27 34 
9 64 28 36 
10 158 82 76 
i 91 61 30 
Z 149 94 55 
3 96 69 27 
4 167 123 44 
5 102 67 35 
6 112 85 2 
7 67 54 13 
8 29 10 | 19 
9 22 11 11 
20 41 ils; 28 
i 41 12 29 
P, 52 34 18 
3 28 df Zi 
4 22 11 iit 
) 19 4 15 
6 24 7 NG 
7 18 0 18 
8 6 2 4 
9 3 3 0 
30 9 5 4 
1 15 9 6 
2 63 47 16 
3 23 15 8 
4 20 13 df 
5 7 ily 0 
6 10 10 0 
it 6 2 4 
Total 2516 1471 1045 


epicentres merely by seeing where their areal 
densities are high and simply state there is a 
group here and there is another there. As 
will be illustrated later, there are several in- 
stances, however, that although two adjoining 
areas show equally high seismic activities, 
earthquake frequencies in them vary in no 
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sympathetic way with time indicating that 
any direct causative relationship is not likely 
to be existing between the seismic activities. 
In a case like this, if earthquake epicentres 
are plotted on a map on the basis of data of 
ten odd years, two such areas with numerous 
epicentres will appear to form one apparently 
single group of epicentres. But this process 
of grouping cannot be meaningful if we are 
interested in causative relationship which 
might be existing among the occurrences of 
earthquakes. 


Chuji TSUBOI 


What will be described below is the result 
of an attempt to establish groups of earth- 
quake epicentres in and near Japan on a phy- 
sical basis. By “ona physical basis ” is meant 
here that particular effort is going to be made 
to find groups of epicentres which are not 
only located in a more or less confined area 
but also there are reasons to believe are cau- 
satively related. Each of such areas, if es- 
tablished, will be called an “ Earthquake Pro- 
vince ”. 

§4. As was stated before, Japanese islands 


Table II. Yearly number of 
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have been divided into a number of compart- 
ments as shown in Fig. 2. The yearly num- 
bers of earthquake epicentres in each of the 
compartments are listed in Table II]. This 
will be the datum for the subsequent studies. 
First of all, the year-to-year variations in 
earthquake frequencies in various compart- 
ments will be compared. 

Let the number of earthquakes in the /-th 
compartment in the ¢-th year be denoted by 
Ni, t). Generally speaking, MN(/, t) and 
N(i+s, t) vary with ¢ similarly if s is small 


epicentres in various compartments. 
Bel eae Pee ioe: 
18 | 22AP2Sue2A 25: \e2Og) 


Peat tae) ey 


16 | 17 27 
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(that is if the two compartments are close to 
each other) and do less so if s becomes larger 
(that is if the two compartments are farther 
apart). This phenomenon is well illustrated 
by an example shown in Fig. 4 in which the 
yearly earthquake numbers in the compart- 
ments (4), (5), (6), (7), (8) and (9) are compa- 
ratively shown. While the graph for M65, 2) 
is very similar in shape to that for MN(4, f), 
that for N(6,¢f) is a little less similar to it, 
and that for N(7, t) is still less similar and so 
on until that for N(9, t) shows no similarity 
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O—O A 1912-1915 
KK B 1912-1933 
e-* C 1934-1955 
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Fig. 3. Number of epicentres in various compartments. 
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Year 
Fig. 4. Yearly earthquake numbers in compartments (4), (5), (6), (7), (8) and (9). 


at all. 


where x and 7 represent respectively the 
A useful numerical measure for expressing 


nerica average values of « and y taken over the 
the degree of similarity of two variable quan- whole range of variables. If z(t) in the above 


tities 2 and y is the correlation coefficient expression is representing the number of 
r(x, y) which is defined as follows: earthquakes in the /-th compartment in the 
oe D(e—z)y—9) t-th year which is N(/, t), and y(t) that in the 

V Maz? Sy—zyP’ m-th compartment which is N(m, ft), the cor- 
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Fig. 5. Symbolic representation of correlation coefficients. 


relation coefficient between Ni/, t) and Nim, t), 
which will be denoted by ,r(/, m), 


=I, t)- ND) }{N(m, t) - N(m)} 


a VEIN, )-NOPE(N(n, )— Nom)? 


will represent the similarity of year-to-year 
variation in earthquake numbers in the two 
compartments / and m. In calculating the cor- 
relation coefficients 7(4, 5), 7(4, 6), 7(4, 7), r(4, 8) 
and r(4, 9), m in the above formula is 44 in 
our case, since we are dealing with yearly 
numbers of earthquakes during 44 years from 
1912 through 1955. The correlation coeffici- 


Table III. Correlation coefficients 7(4, ™). 
r(4, 4) 1.00 
r(4, 5) | 0.74 
r(4, 6) | 0.62 
r(4, 7) | 0.04 
r(4, 8) —0.05 
r(4, 9) Says 


I 


Pe 


ents which have thus been deduced are given 
in Table III. 

The correlation coefficient is seen to de- 
crease monotonously with the distance be- 
tween the two compared compartments. 

In order to investigate the problem from a 
more general standpoint, that is to find out 
the rule of decrease of the correlation coeffi- 
cient with distance, rv(/, m) for all possible 
combinations of / and m should be known. 
For this purpose, Table IV has been prepared. 
The value of r(J, m) is to be found at the in- 
tersection of the /-th column and the m-th 
row of the table. The values given in the 
table are naturally symmetrical with respect 
to the diagonal from upper left to lower right. 
Fig. 5 shows the same thing symbolically by 
means of different marks as indicated in it. 

First, of all, we wish to investigate how the 
correlation coefficient depends on the distance 
between the two compared compartments. 
For instance, corresponding to the unit dis- 


tance (40 km), 
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Correlation Coefficient 


Fig. 6. Average correlation coefficient 7 as a function of distance s. 


7(122, 122+1)=r7(122, 123) Table V. Average correlation coefficients 
r(123, 123+1)=r(123, 124) for various s. 
Ae Mie nah BS e 

7(36, 36-+1)=7(36, 37) “pA » be eeh- ep Ree ee 
should be averaged. Generally the value of oo 1.00 22 —0.05 
y(l, 14s) averaged over / will give the cor- ae 0.34 3 —0.03 
relation coefficient corresponding to the dis- 2 0.19 4 0.02 
tance s and this should be calculated for all 3 0.13 5 0.04 
of s=1, 2,---, 43. The results of calculations 4 0.06 6 —0.03 
are given in Table V. The values given in 5 | 0.04 7 —0.04 
Table V are the weighted means of r(/, /+s) 6 | =On02 8 —0.06 
averaged over /, the weight being given propor- 7 —0.01 9 0.00 
tional to the sum of earthquake numbers in 8 —0.02 30 0.00 
the compartments / and (/+s). Fig. 6 in which 9 —0.06 1 0.04 
these average values are plotted clearly shows 10 ~0.02 2 0.08 
that the correlation becomes lower and lower 1 ~0.03 3 0.00 
with distance (Tsupor: 1949). The coefficient 2 0.04 4 0.05 
is naturally 1.0 for s=0 and decreases smoothly 3 _0.02 5 oon 
with distance until it is almost 0 at s=4 or 4 0.00 6 ’ 

; ’ . —0.06 

5. Since the compartment is each 40 km long, 5 Gage 
the conclusion is that earthquakes taking 6 Aras 8 ee 
place within the range of 40x4=160km or : } cS 
40x5=200km are positively correlated in | ae nee 
their occurrence and that they may be regarded a, bibs oe ae 
to be common offspring of some specific seis- : | oe | Lue ee 
mic condition prevailing within and not be- - ae - | —0.05 
yond this extent. The site of this “some : —0.01 3 a 0.40 


. . : Oat ” ‘ a 
specific seismic condition ” may appropriately 
be called a seismic field. It may be added have a small maximum around s=20. This 
that the curve shown in Fig. 6 appears to point will be discussed later. 
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§5. What has been stated above is just an 
apparent conclusion which could be drawn at 
once from the materials given in Table IV. 
But going into a little further detail, Table 
IV involves several other facts of interest. 
First of all, the correlation coefficients be- 
tween earthquake numbers in two adjoining 
compartments (s=1) are extracted in Table 
VI and shown in Fig. 7, together with earth- 
quake numbers involved. The average of the 


Table VI. Correlation coefficient » between 
earthquake numbers in two adjoining 
compartments. 

Comp. | A Total earthquake 
number 

122 12 
3 ae 29 
4 | (23 36 
5 ee | 89 
6 tee 120 
,, aa 125 
128 er | 99 
1 wer 61 
: 0.44 a 
3 0.59 S 
4 0.74 | _ 
5 os | 61 
6 0.23 | a 
7 0.04 2 
8 0.12 60 
9 0.21 62 
a | 0.35 a 
l 0.48 - 
2 | 034 146 
3 | 0.45 - 
4 0.47 164 
5 poe 101 
6 0.51 112 
7 0.16 fe 
8 0.03 oe 
9 0.22 & 
20 0.69 he 
I 0.22 | Ea 
2 0.22 os 
3 | 0.26 | 2 
4 0.08 ae 
5 0.32 oe 
6 0.68 oe 
i 0.23 te 
: —0.08 : 
0.01 9 

30 20012 
1 0.30 1S 
2 0.16 os 
3 0.44 a 
. 0.24 
2 ~0.15 a 
S =e} 6 

average 0.34 


a 


coefficients is 0.34 but each individual value 
varies considerably along the island arc. The 
largest coefficient is 0.80, while the lowest is 
0.01. 

It is noteworthy that the coefficient does 
vary in no parallel way with earthquake num- 
bers. There are instances that although two 
adjoining compartments both show high seis- 
mic activities, the correlation coefficient be- 
tween them is low, like 7(127, 128)=0.09 with 
224 earthquakes, while there are other in- 
stances that although two adjoining compart- 
ments both show low activities, the correla- 
tion coefficient is high, like 7(26, 27)=0.68 with 
39 earthquakes. This fact indicates above 
all that the correlation coefficient is not a 
function of the distance between the two com- 
pared compartments alone. 

The above argument is just regarding two 
adjoining compartments. In order to investi- 
gate this problem from a more general stand- 
point, let us look into Table IV again. In 
Table IV, the correlation coefficients that are 
larger than 0.3 in absolute values are printed 
in Gothic type. The mode in which these 
relatively high correlation coefficients are ar- 
ranged in the table and in Fig. 5 is remarkable 
in that they are clustered to form several 
squares along the diagonal (Fig. 8a), instead 
of dispersing generally away from the diagonal 
with statis-tical irrezularities (Fig. 8). 

If the correlation coefficient is a function of 
the distance between the two compared com- 
partment alone, the latter distribution (Fig. 
8b) is likely to occur but actually this is not 
the case. This peculiar distribution of rela- 
tively higher correlation coefficients would in- 
dicate that the several compartments belong- 
ing to each one of the square groups are close- 
ly related in seismic activities showing sym- 
pathetic variations. In this sense, these com- 
partments may be said to form an “earth- 
quake province ”. 

Referring to Table IV from this standpoint, 
six earthquake provinces may be established, 
one of which, however, is doubtful. These 
provinces are listed in Table WHT 


$6. The overall average of the correlation 
coefficients is 0.48, if the compartments {are 
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ae 
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Correlation coefficient between earthquake numbers in adjoining compartments. 
m: Earthquake numbers. 


both inside an earthquake province, while it 
is 0.13 if they are both outside it. 

The overall decrease of the correlation co- 
efficient with distance was shown in Fig. 6. 
In that calculation, no distinction was made 
as to whether or not the two compartments 
compared are within one earthquake province. 
Now that six earthquake provinces have been 
established, the decreases of the correlation 
coefficient with distance inward and outward 
them can be calculated separately. They are 


(a) 


Table VII. Earthquake provinces 


. | . . 

Province | Compartments compared in Table VIII and in Fig. 9. The 
i. a 125 126 127 ; decrease is notably different in the two cases. 
B 9 3 gies Se Inward the earthquake province, the correla- 
C 11 2 13 1 15 16 tion coefficient does not decrease greatly with 
D (doubtful) {9° 3g San distance while outward it falls sharply down 
: 25 96 27 to a very small value. 

Bi oe ey ny, All these findings equally point to the phys- 
i ical significance of the earthquake province. 
Table VIM. Correlation coefficients inside and outside earthquake province 
Inside earthquake | | Outside earthquake 
| province | | province 
Distance (unit H=40km) 1) 5 | 4°) 3] 2 | 2) OG.) 1 Weegee ee 
Lie ancy a, ‘a SS ie| — = = = | | = - 
Correlation coefficient 0.28 | 0.34 | 0.42 | 0.42 | 0.46 | 1.00 | 0.10 | 0.03 | 0.04 | 0.04 | 0.03 


a ee ees ee See 
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s=5 
200 km 


inside 


s=5 
200 km 


outside 


Fig. 9. Correlation coefficient inside and outside 
earthquake province. 


The earthquake province may be said to re- 
present a stationary site over which one seis- 
mic field continues to exist, although its in- 
tensity may vary according to time. Earth- 
quake numbers in the compartments which 
belong to the same one earthquake province 
are commonly affected by this field intensity 
and this will explain why they vary sympathe- 
tically. The earthquake province is sharply 
bounded on its both sides and the boundaries 
appear to act as barriers which are hard for 
a seismic field to trepass. 

The earthquake numbers in each of the six 
earthquake provinces are given in Table IX. 
The correlation coefficients among earthquake 
numbers in the six provinces are given in 
Table X. From the low values of coefficients 
given in the table, it is obvious that the seis- 
-mic activities in the six earthquake provinces 
are almost independent events. 
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Table IX. Earthquake numbers in various 
earthquake provinces 


a enh a eatin tee enn 
[arn B re D E F 
[125-127 | 2-6 -| 11-16 | 19-21 | 25-27 | 31-34 

1912 | 10 q. | 286 2 3 6 
3 9 5 | 30 2 0 | 13 
4 5 9 | 30 0 0 2 
5 3 22 | 40 2 1 2 
6 5 16 | 27 0 0 6 
7 Fi 13 | 24 2 0 0 
8 5 Aer Oi 2 0 6 
9 0 alts 3 2 3 
20 14 Se 6] 1 0 3 
1 12 al at) 0 0 3 
2 9 4 | 26 0 1 2 
3 8 ae 6 0 9 
4 10 3 | 54 4 2 3 
5 6 1 i 24 0 0 2 
6 11 ae oats 3 0 3 
7 9 Slee 2 : 1 
8 3 1a 2 1 | 2 
9 6 4 9 2 0 4 
30 14 2 8 0 0 4 
1 14 10 | 15 0 0 5 
2 10 7 3 0 0 4 
3 7 82 2 4 0 1 
4 y) yi 2 0 0 3 
5 9 13 4 1 0 1 
6 2 0 3 0 0 0 
is 2 6 Z 1 0 6 
8 1 Ai 19 0 3 0 
9 6 | 6 8 0 0 2 
40 1 5 HT 0 1 0 
1 3 7 9 4 1 3 
2 6 0 9 5 1 4 
3 3 tae 4 ig Oe 1 
4 5 2 2 1 0 1 
5 6 2 ff | Scab 0 0 
Giulio 6 3 7 2 7 1 
Pele 2 5 als 4 
8 2 a 9 | 13 2 4 
9 10 2a ae 1 2 2 
50 iG 2; 5 2 2 2 
i oleae 2 8 4 0 1 
ee RS) 27 13 3 1 0 
sa 17 Ra eALS eS 1 1 
Nae cy alka 0 1 1 

1955 ss 6 113) 1 10 | 20 3 0 

eT 


Table X. Correlation coefficients between earthquake numbers 
in various earthquake provinces. 


| A B 
A 1.00 0.16 
B 0.16 | 1.00 
c Gries 0.07 
D =6 11 ~0.02 
E Oni —0.14 
F ~0.04 EO Ad 


c D E F 
0.13 zi. 11 -0.11 | 0.04 
—0.07 =.02 =O, 14 —0.14 
1.00 —0.05 | -—0.08 0.38 
—0.05 1.00 0.05 —0.03 
—0.08 0.05 1.00 Aig 
0:38 —0.03 =Oeit 1.00 
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§7. Referring back to Fig. 5 again, it will be 
noted that the correlation coefficient becomes 
a little higher at the distance of s=20 or 
800km. Fig. 10 is simply a diagrammatic 
rearrangement of Fig. 5 and the correlation 
coefficients larger than 0.3 in absolute value 
are specifically indicated. The diagram shows 
at what distances the relatively high positive 
correlation coefficients appear in relation to 
each of the compartments. The lower curve 
in the same figure shows the number of such 
coefficients divided by the total number of 
coefficients corresponding to that distance. 


The percentage is very low between 6th to 
14th compartment, and beyond that a blunt 
but well recognisable maximum is seen at 
about 18th—-20th compartment. The question 
what this blunt maximum really means is as 
yet hard to answer. 
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A New Formula connecting Magnitude 


and Number of Earthquakes. 


By 


Chuji Tsusot 


Geophysical Institute, Tokyo University 


Abstract 


A new formula 


log n=(1— M/M,,) log 7 


connecting the magnitude M and the number x of earthquakes is proposed, where 


Oe number of earthquakes having a magnitude 0, 


This formula represents observational facts very nicely. 


Agee raat magnitude of the largest conceivable earthquake. 


It involves only one constant 


log »,, in a contrast with the formula by GUTENBERG and RICHTER 
log N=a+b (8—M) 
which involves two. a and 0 in the above formula have beee found not to be quite independent. 
The new formula implies that 2(M) is equal to m/, instead of being equal to 


n,-¢(M), as is usually taken for granted. 


§1. B. Gurenserc and C. F._ RICHTER 


introduced a formula 
log N=a+b(8—M) 

for expressing the relation between N and 
M, where N is the number of earthquakes 
which took place in a certain region dur- 
ing a certain length of time and which 
were M+t4M/2 in magnitude. 4M is the 
interval of steps by which ™ is classified. 
An important feature of this formula lies in 
the fact that log N isa linear function of M 
and many subsequent studies have shown 
that this is true for most earthquake regions 
of the world. 

In this formula, the constant @ obviously 
represents log Ns, which is the logarithm of 
the number of earthquakes having a magni- 
tude 8+4M/2. The constant 6 is equal to 
log (Nu/Nu+i) and represents the slope of 
log N with respect to M. The constant 0 
is generally considered to be independent of 
both a and time and also to be a parameter 
which is characteristic to each seismic region 
concerned. If this is really so, this amounts 
to mean that the number of earthquakes of 
every magnitude increases with time in the 


log N 


Fig. 1. 


same proportion, so that the results will bea 
parallel translation of the line in the (log NV) 
-M diagram as shown in Fig. 1. 


log N 


Fig. 2. 
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Now suppose the seismicities in two regions, 
say X and Y, are represented by the two 
lines x and y in Fig. 2, of which the slope is 
gentler for x than for y. Suppose also that 
in a long run of time, the earthquake number 
in X increases, and the line x is shifted to x’ 
until N) which is the number of earthquakes 
having zero magnitude in X, will coincide 
with N, in Y. If the parallel translation of 
the line is true, the natural consequence will 
be that the line x’ will be lying above the 
line y everywhere, indicating that earth- 
quakes having greater magnitudes have been 
more numerous in region X than in Y. 
Although the time involved is of different 
length, can this statement be easily ac- 
ceptable ? 


§2. In their book ‘‘ Seismicity of the 
Earth’’, GutTensperc and Ricutrer (1949) de- 
scribed the results of their investigation of 
the number of shallow earthquakes in various 
regions of the world in relation to their 
magnitudes. The values of a and 0 in the 
formula 


log N=a+b(8—M) 
as given by them are listed here in Table I. 


Table I. Values of @ and 6b (GuTENBERG 


and RicHErR) and a. 


Religion a | 8 t eee 
Alaska = el dee 41 8.9 
Sr yeah ie 21240 Sadqies Tat 
South America = tod 1.0.45.) 936 4.1 
Kermadec Islands —2.3 | 1.3 34 9.6 
Solomon Islands —1.23) 1.01 | 36 8.4 
Japan —0.90) 0.80 | 42 fall 
Sunda Islands p17 OF9) 39 flee 
Pamir to Eastern Asia —1.7 | 0.6 41 4.7 
Turkey —2.1 | 0.9 41 6.7 
Atlantic —2.4.| 1.4 25 10.2 
Indian Ocean —2.4/1.3 | Al 9.6 
eee ee 


On investigating the values listed in Table 
I, one will find that @ and 6b are not quite 
independent. The larger the absolute value 
of a, the larger is the value of 6 generally. 


The values of a and b corresponding to the 
11 regions are plotted in Fig. 3, in which the 
suspected dependency may be seen. It will 
be understood at least that there is a reason 
to suspect if a@ and b are dependent rather 
than are not. 

It is to be noted here that the values given 
in Table I apply to the mean annual number. 
of earthquakes in each of the 11 regions and 
also than 4M was taken to be 1/4. 


| 
| 


b 


2.0 


1.0 


Bigce. 


In deducing the mean annual number of 
earthquakes, the two authors must have 
divided the total number of earthquakes in 
each of the regions by the number of years 
for which the statistics was available for 
that region. The column under a heading ¢ 
in Table I gives the number of years which 
was presumably taken in the statistics of 
GUTENBERG and RICHTER. 

If represents the total number of earth- 
quakes in ¢ years, it must be related to N as 
follows : 


n= Nt 
or 
log n=log N-+log t. 
If 
log N=a+b(8—M) 


as GUTENBERG and RIcHTER suggested, then 
log n=a+8b+log t—bM, 
=a—pM, 
where 


a=a-+8b+log t 
and 


B=b, 


The values of a@=a+8b+logt are also 
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0 50 


10.0 


x 11 Regions in the statistics of GUTENBERG 


and RICHTER (Table I) 
@ Japan (Table II) 
© 3 Parts in Japan. 
Fig. 4. 


given in Table I. It will at once be noticed 
that the larger the value of a, the larger is 
the value of 8 also. The values of a and B 
are plotted in Fig. 4 which shows a beautifully 
linear relation between them. The relation 
may approximately be expressed by an 
equation : 
a=8B. 
§3. Such a dependency as this can happen 
when 
a+ logt=0. 
Since a is the mean annual number of earth- 
quakes having a magnitude 8+4M/2, we may 
write 
log Ns+log t=0 
or 
Ns=Nz Ls. 
This implies that the number of earthquakes 
having this magnitude is only one, irrespec- 
tive of region and time. This cannot of 
course be exactly true, but this would seem 
to mean that the number of earthquakes 
having this magnitude does not differ much 
from a region to another. As a matter of 
fact, the (log m)—M lines for the 11 regions 
which are drawn in Fig. 5 do appear to 
converge to a point at about M=8 on the 
M-axis. Although the number of earth- 
quakes having zero magnitude differs con- 


log n 
10.0 
5.0 
bathing SS 
> 2 
2 4 6 8 
Bigs. 


siderably from a region to another, that 
having a magnitude about 8 does not. 
If the relation 
OO 
is put into the formula 
log n=a—BM, 
we get 


This formula, which works very nicely for 
all the regions, involves only one constant 
log m, in a striking contrast with the classical 
formula 

log N=a+b (8—M) 
which involves two constants @ and 0. 

In the new formula, 8 is approximately the 
magnitude of the largest conceivable earth- 
quake so that it will be written as Mz. 
Thus we get finally 

log n=(1—M/M;) log mo, 
where both m and are functions of time, 
of course, and the ratio of m and m will also 
be a function of time such as 
log (/10)= —M/M,z log no. 
According to this new formula, the parallel 
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translation of the (log )—M line will not 
take place. Instead, the (log ~)—WM line will 
be shifted with its one end being hinged at 
one point on the Maxis. This means that 
smaller earthquakes will increase more and 
more in percentage with time. 


§4. In 1952, the present writer (Tsusor: 
1952), referring to earthquakes in and near 
Japan, investigated how the constants @ and 
6 in the formula 
log N=a+6(8—M) 
tend to change when more and more years 
are added into the statistics. The results of 
the study are listed in Table H. They can 
at once be used for deducing the values of 
the constants a and § in the new formula: 
log n=(a+8b--log t)—bM 
=a—BM. 

The values obtained are given in Table II. 

The values of @ and @ are plotted in Fig. 


4 also. Again, the points lie on a straight 
line 
A—Ssi 
Table II 

ree a aege a=a+8b 
+log t 
1931 — 0.40 0.52 Sykes 
1931-32 — 0.97 79 5.65 
31-33 — 0.67 60 4.61 
31-34 — 0.80 65 a0) 
31-35 — 0.86 69 5.36 
31-36 — 1.05 79 6.05 
31-37 — 1.25 90 6.80 
31-38 — 1.10 85 6.60 
31-39 = 1.15 86 6.68 
31-40 — 1,28 90 6.97 
31-41 = 1.25 91 7.07 
31-42 — 1.30 95 7.38 
31-43 = 1.35 98 7.60 
31-44 — 1.36 99 Chesiial 
31-45 — 1.44 1.04 8.06 
31-46 ao dhe ote, 1.00 7.89 
31-47 =i 1 aul 1.01 7.90 
31-48 — 1.43 it OR 8.07 
31-49 — 1.46 1.04 8.14 
31-50 — 1.46 1.04 8.16 

ar ee Oe 


Also in the same previous paper, the whole 
area of Japan was divided in three parts A, 
B, and C and the constants @ and 6 were 
determined for each of the three parts 
separately, using a 20 year statistics of earth- 
quakes. The results obtained were: 

A: log N=—1.60+1.06 (8—M), 
B: log N=—1.57+0.72 (8—M), 
C: log N=—1.61+0.66 (8—M), 

The values of a which are deduced similarly 

are as follows: 


A: 8.18 
B: 5.49 
Ce 4.97 


They agree very well with the correspond- 
ing values of 88 which are 8.48, 5.76 and 5.28 
respectively. The three points plotted in Fig. 
4 show this relation clearly. If we take a 
value of 7.6 instead of 8, the agreement is 
still better. 

In the same previous paper, it was stated 
that “....if the total number of earth- 
quakes in a certain area is twice as much as 
those in another, for instance, it does not 
mean that earthquakes of various magnitudes 
are equally double in number. Earthquakes 
with smaller magnitudes increase more 
rapidly than those with greater magnitudes 
dOs -s5ahe 

All the above considerations equally point 
out that the distribution of earthquake 
numbers among various magnitudes is a 
feature which is neither stationary with time 
nor characteristic to a seismic region. 


If the formula 

log n=(1—M/M,) log ny 
is once accepted, the next problem will 
naturally be to understand what it really 
means. The classical formula 

log N=log Ns+6(8—M) 

can be written as 
log N=log Ns+8b—bM 

=log No—bM 


§ 5. 


so that 
N= Nod(M). 
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N and N, will be proportional and the pro- 
portional factor #(M/) depends of M. On the 
other hand, the new formula can be trans- 
formed into 


n= No M/ Mx 


nN=No tC") 


In this case, log » and log m will be pro- 
portional and the proportional factor f(M) 
depends on M. 

Although the conclusion reached appears to 
be a rather strange one and is hard to ap- 
prehend from our usual understandings, it is 
a logical and perhaps an inevitable conse- 
quence of facts. 

Here we will recall that a linear formula 
between log E and M was suggested by 
GUTENBERG and RicHTER, such as 

log E=p+qM 
where E is the energy of an earthquake 
having a magnitude M. p represents ob- 
viously log &) and if we write Ez for the 
energy of the greatest conceivable earthquake, 
we get 
log E,=log Ly) +qMz, 


and 
M=(log E—log E,)/q 
M,=(log Ez—log Ey)/q 
and therefore 


y—M _1_ log E —log Ey 
_ log (E /Ex) 
~ log (Ey/Ex) 


Combining this with the new formula gives 
finally 


_ log (E/Ez) 
~~ log (Ey/Ex) 


and we see that 


f(M) 


log n log no 


_ log (E /Ez) 
log (Eo/Ez) ° 
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On the Sensitivity Ratio of Clay 


By 


Shin’ichi YAmacucui 
Physical Department, Defence Academy, Yokosuka 


Summary 


By taking as sample Osaka alluvial clay which is called Kaolinite and the clay gathered 
from Okuramura land creep region, Yamagata prefecture, the author studied the sensitivity 
of clay. As the result, we could comprehend that it is due to the special water in case of 
the former clay, while in case of the latter clay it is due to the difference in the arrange- 


ment of the scale-like particles. 


Introduction 


§1. 


The clay generally lessens its strength only 
by remolding, without changing its water 
content. The ratio Sv/Sr, in which Sv is the 
strength of the undisturbed clay and Sr is the 
strength of the same clay fully remolded, 
with its water content unchanged, is called 
the sensitivity ratio. Since by remolding 
the clay changes its nature and the critical 
point of shear failure becomes obscure, the 
strength of the remolded clay Sr is deter- 
mined from the stress for that strain which 
causes shear failure of the undisturbed clay. 
The sensitivity ratio varies considerably with 
the kind of clay, ranging from 2 to 500. So, 
the clarification of this property is important 
not only for soil mechanics, but also for geo- 
physics, because several natural phenomena 
are associated with the disturbed soil. Two 
theories have been advanced to explain the 
_differences in the property : the theory of A. 
Casagrande and that of K. Terzaghi. 

Casagrande’s theory is as follows: in a 
clay/silt mixture those clay particles which 
enter the lesser spaces between the silt 
particles are consolidated by higher stresses 
than those which lie in the larger spaces. 
This ‘bond clay’ possesses considerable 
strength. When such soils are remolded, 
the total strength decreases, because the 
strength due to bonding is destroyed. 

K. Terzaghi’s theory states that: if the 
soil is consolidated for a sufficient time, the 
stresses squeeze out the absorbed water, so 


that the particles come into contact with one 
another and at such point the true cohesion 
force functions. If any pressure applied 
from the outside is removed, the absorbed 
water returns gradually to the contact points 
of the particles, but it takes a fairly long 
time to do so completely. When the soil is 
remolded, the total strength becomes less, 
because the distance between the particles 
becomes greater. 

Recently, O. Morretto“) demonstrated that 
the remolded clay put in a state of rest 
gradually recovered its strength, and the 
strength ratio of the remolded clay after 
resting for the remolded clay became larger 
as the liquidity index of the clay sample 
became larger. Newland et al) showed that 
the sensitivity ratio becomes smaller as pre- 
consolidation pressure is larger. Recent 
experiment by the author indicates that a 
paste which is consisted of the microparticles 
of clay below 1 when in a state of rest 
increases its strength just in the same way 
as does the remolded clay containing all clay 
and silt particles. 

According to Terzaghi’s theory, the strength 
of the remolded clay declines as the time of 
resting increases. And according to Casa- 
grande’s theory, the strength of a paste 
which consists of the microparticles of clay 
below 1 does not change when in a state of 
resting. Therefore both theories cannot 
always explain those experimental results. 

Thus, the author attempted to clarify the 
problem of sensitivity. The samples he 
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used were taken from Osaka alluvial clay 
which is chiefly consisted of kaollinite and 
from Okuramura clay of Yamagata prefecture 
which is chiefly consisted of montmorillonite. 


§2. Osaka Alluvial Clay 


() Sensitivity 

Many completely remolded specimens, 3.4cm 
in diameter and 8cm in height, were made 
and each was placed in a rubber sleeve 
closed with latex. Then, varying the period 
of resting for each specimen, their strengths 
were tested. Unconfined compression test 


-3 
x/0 cm 


Displacement — 


0 / 2 3 4 le 

Load —~> 

Fig. 2.1 Stress-strain relation in stress control 

test (6A6-6, 916-S) Figures in the diagram 
indicate resting days. 


Strength > 
Db 
Electric Resistance > 


1 Lo 2 
/ 10 1? Days 10° 


Resting Day —> 
Fig. 2.2 Recovery of strength and of electric 


resistance. 
ee aw WANS Rouge Resistance 
| Chadiseurbed Clay @ y% 
| Remiia Clay _ ©) TX 


was carried out in the controlled stress 
manner. The rate of increase of stress was 
200gr/cm? min. The sensitivity ratio deter- 
mined in this way was about 10. 
recovery of strength proceeds gradually, 
merely about one half of the strength was 
recovered within two months. Part of the 
results of the controlled stress unconfined 
compression test is shown in Fig. 2.1; and the 
recovery of strength, in Fig. 2.2. 


(2) Void Ratio and Electric Double Potential 
of Particle Surface 

Of the clay just before being remolded, of 
the clay immediately after being remolded 
with its water content unchanged and of the 
same clay after various period restings, the 
strength, void ratio and electric double sheet 
potential of particle surface were measured. 

Void ratio should be calculated indirectly 
from the measured values of water content, 
specific gravity of soil particle and wet 
density of the clay. However, in the 
mensuration of wet density some errors may 
easily be accompanied due to the measure- 
ment of volume. Therefore, to obtain the 
more accurate value in the measurement of 
volume, the author adopted the following 
method ; a loaf of clay wrapped with paraffine 


was steeped into water and then its bouyancy 
was measured. 


Electric double sheet potential of particle 
surface was calculated from the electro- 
osmotic equation, 

V=lrAD/4n 
in which V is the volume of water, passing 
through a cross section of the cylindrical 
tube with the uniform diameter, per unit 
time and € is electric potential difference 
between the fixed plane and the movable 
plane of the double sheet, 7 is the viscosity 
of water, 7 is the inner radius of the tube, 
D the dielectric constant of water and H is 
voltage gradient between the electric poles. 

As both D and 7 being already known and 
7 possible to be estimated from the porosity 
by considering the void in the soil as a 
vacant, cylindrical tube, ¢ can be calculated 
if H and V are measured. Of the values of 
¢ thus obtained, there was hardly recognized 
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any difference between the undisturbed clay 
and the remolded clay. 


The void ratio of the remolded specimen 
proves a little larger than that of the un- 
disturbed specimen. 

On the other hand, from the linear rela- 
tionship of logarithums both of the strength 
and of water content of the saturated un- 
disturbed specimen,) the relation between 
the void and the strength can be determined. 
As the result, it was comprehended that the 
decrease of strength by remolding couldn’t 
be explained according to such a small dif- 
ference of the void ratio. 


(3) Electric Resistance 

Preparing the specimen in the same way as 
in the unconfined compression test and placing 
the electric poles, which are made of stainless 
steel and are insulated from the compression 
test apparatus, at the both ends, its resistance 
was measured by means of a Kohlraush 
bridge. As resistance does not change with 
frequency of the alternating electric current 
in the range from 50 to 10,000 cycles, 10,000 
cycle alternating current was used. When 
the loads were put on the specimens during 
the experiment, some change appeared in the 
electric resistance. This may be interpreted 
due to a change of the contact area between 
the electric pole and specimen, so in carrying 


Electric Resistance > 


Water Content —~> 
Fig. 2.3 Relation between electric resistance 
and water content (6A6-6, 920-S) 
Undisturbed clay @ 


Remold clay © 
Volume is constant —— 


Air is not included ~—..---. 


Theoretical curve { 


out the experiment, a load placed on the 
specimens was increased until the apparent 
electric conductivity reached saturation. 

Many measurements were taken thus, and 
one example is shown in Fig. 2.3. Since the 
specimen measured in the loaded state differs 
from the original one both in diameter and 
in height, the adopted value is corrected to 
the original dimensions. If the size of the 
measured specimen is de™ in diameter, Icom 
in height and its electric resistance is R, the 
converted electric resistance R’ in terms of 
the standard size of the specimen must be 
shown as follows ; 

8 
eve 

In considering the result, we assume at 
first that the electric conductivity is concerned 
only with the contained water and that the 
clay particle is nonconductive body. By 
doing so, we determine the electric resistance 
for the volume-ratio of water in the clay 
alone. 
Thence, the following equation can be ob- 
tained ; 


rR’ 


esx Cane” (2.1) 
where R;,, and Rs. are respectively the electric 
resistances, per unit area of the vertical 
section and per unit length, of the contained 
water and of the clay, and 7, is the volume 
ratio of contained water in the clay. When 
the water content is 65%, the clay consists 
merely of the particle and water, but if the 
water content becomes less, two extreme 
cases ought to be considered: one in which 
the volume shrinks as much as the water 
content decreases and air does not enter into 
the specimen, and the other in which as much 
air enters as water is lost and the volume 
does not shrink. The relation between the 
water content and the electric resistance for 
these two extreme cases calculated from the 
equation (2.1) is entered in Fig. 2.3. 

The actual measured values when desiccated 
naturally in the atmosphere, indicate that 
the result accords, in case of the remolded 
clay, with the former instance where air 
does not enter into the specimen and the 
volume shrinks and, in case of the undisturbed 
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clay, with the latter instance where the 
volume does not shrink and air enters. 
Measuring the relation between the water 
content and the void ratio for the case in 
which the specimen is desiccated naturally in 
the atmosphere, we find it to be as stated. 
(Fig. 4:1): 

Here, we discover a fairly distinct dif- 
ference between the remolded and the un- 
disturbed clays. That there is some dif- 
ference in the electric resistance between 
the disturbed specimen and the remolded 
specimen, even if the water content or the 
void ratio of both specimens may be equal, 
is probably resulted from the fact that there 
is difference in the specific electric resistance 
of contained water between the two. But, 
scrutinizing Fig. 2.3 more closely, we may 
perceive that, if the water content is trans- 
ferred about 20%, the volume—constant 
curves of the remolded specimen and of the 
undisturbed specimen overlap. This suggests 
that in the undisturbed sample there is about 
20% water which evaporates in the desicca- 
tor, yet acts as a non-conductive body to 
electricity and that between the undisturbed 
clay and the remolded clay there is no dif- 
ference in the electric conductivity of water 
excluding 20% portion afore-mentioned. 

Of these two conclusions thus far delivered, 
the author advanced to investigate which is 
more appropriate, from another angle in the 
following chapter. 


(4) Specific Heat Capacity 

The specific heat of clay is measured in 
the following ways: the samples, which are 
all made 7cm in diameter, 7 or 8cm high and 
o00gr in weight, are covered with latex and 
placed in 1,000gr of water and heated. For 
heating, the electric heater (100v, 60w) was 
used. And the supplied electric voltage, 
being increased to 100v through the stabilizer 
beforehand, then dropped to 50v by the trans- 
former. By doing so, the rate of the 
increase of temperature becomes 2-3°C per 
20 min., and both clay and water are heated 
almost equally. 

Then, the following equation, where Cu 
Cw and Ce are each respectively the specific 


heats of the vessel, water and the clay, and 
the given thermal energy is En, the rise of 
temperature produced is T;, is obtained ; 


En,/T,;=Cv+Cw+Cc (2.2) 


And in the experiment performed exactly 
in the same way without the clay specimen, 
the equation 


En2/T2=Cv+Cw (2:3) 


is derived, in which J, is the rise of tem- 
perature. Here, from (2.2) and (2.3) Ce can 
be obtained. 

Next, by determining the water content of 
the clay and subtracting the heat capacity of 
the contained water, provided that water 
contained in the clay is invariably common, 
we can obtain the specific heat capacity of 
the clay particles alone. As is shown in Fig. 
2.4, the specific heat of the undisturbed 
sample is twice or thrice as much as that of 
the remolded sample in spite of the fact that 
there is no difference in the specific heat of 
the clay particles. 
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Fig. 2.4. Specific heat of clay particle (Taisho- 
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These values seem to be too high in the 
case of the undisturbed sample and also, 
when compared with the specific heat of rock 
or of sand, 0.17-0.23. 

These results come forth from estimating 
the specific heat capacity of all contained 
water in the undisturbed clay as 1.0, and yet, 
in reality, the special water as was already 
mentioned in (2.3) and the common water of 
which specific heat capacity is 1.0 are both 
comprized in the contained water of the 
undisturbed clay. Therefore, by presuming 
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that the special water turns into the common 
water as temperature rises and that in this 
process much latent heat is needed, we can 
elucidate the difference of heat capacity 
between the undisturbed clay and the remolded 
one. 

If it be the case, the part of the special 
water must revert to common water when 
temperature rises, and, accordingly, the 
electric resistance must decrease inversely 
with the rise of temperature. In Fig. 2.5 the 
experimentally obtained relation between 


temperature and the electric resistance is 


shown. 


S 
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Fig. 2.5 Relation between electric resistance 
and temperature 
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It is found, as we expected, that the 
electric resistance decreases as temperature 
rises. Yet this process is reversible, so the 
electric resistance increases again if tem- 
perature falls. Of course, it is well known 
that the electric resistance of solution itself 
“decreases as temperature rises. So, the 
specimens of the undisturbed clay, of the 
remolded clay and of the paste which is 
made from the dried and powdered clay 
particles and pure water, were investigated 
as for the relation between temperature and 
the electric resistance. The relation between 
the ratio of electric resistance at various 
temperature and that at 18°C, is shown in 
Fig. 2.6. Except in case of the undisturbed 
sample, every line lies almost on the same 


curve. 
Examining this relation by taking the 


Ratio —~> 
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Temperature —> 
Fig. 2.6 Relation between electric resistance 
ratio R/R) and temperature. Where R, is the 
resistance at 18°C. 


Undisturbed @ Water Content 40% 
Remolded © Water Content 40% 
Paste [] Water Content 67% 
Paste Z\ Water Content 100% 
Dispersed system = Water Content 500% 


logarithm of electric resistance as the axis of 
ordinate and the reciprocal of absolute tem- 
perature as the axis of abscissa, we obtain a 
straight line relation and the value of its 
gradient. 

Since, except in case of the undistrubed 
clay, the straight line relation and the value 
of its gradient are in accord with the well- 
known relation between the viscosity of 
water (7) and temperature, 7 a exp (-Ew/kT) 
in which k is Boltzmann constant and Ew is 
is activation energy of water, we can assume 
that the change of electric resistance by 
temperature is caused by the change of 
viscosity coefficient of the soil moisture by 
temperature. 

Next, we obtained the ratio of electric 
resistance at 8°C to that at 18°C for each 
water content, as shown in Fig. 2.7. As 
clearly shown in this figure, the case of 
undisturbed clay is distinguishable from the 
other cases. 

The change of electric resistance by tem- 
perature in the paste or in the remolded 
clay, which is observed from the figure to 
occupy about 80% of the change of electric 
resistance by temperature in the undisturbed 
clay, was ascertained as caused by the change 
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Fig. 2.7 Relation between electric resistance 
ratio (at 18°C to at 8°C) and water content. 
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of electric resistance by temperature in the 
soil moisture. Accordingly it is evident that 
some one fifth of the change of electric 
resistance by temperature in the undisturbed 


clay is owing to the change of special water 
into common water. 


Then, the latent heat of special water can 
be calculated as follows : if the fresh saturated 
undisturbed clay, which keeps a water con- 
tent of 40%, is heated from 8°C to 18°C, the 
electric resistance of the specimen falls from 
2500 to 1900. (cf. Fig. 2.5) 

However, this change is not only due to the 
decrease of special water, but to the increase 
of electric conductivity of common water 
also. 

Nevertheless, it may be considered that the 
change of resistance due to the decrease of 
special water alone is approximately (250- 
190)Q/5=120. This change, 120, is equiva- 
lent to the change of 2% of special water 
into common water. (cf. Fig. 2.3) From the 
above consideration it is estimated that a rise 
of 1°C in temperature produces a 0.2% change 
of special water into common water. Since 
the difference of the specific heats between the 
undisturbed and clay the remolded clay (water 
content both 40%) is about 0.3-0.1, not merely 
from Fig. 2.4 but also in consideration of the 
experimental result in (2.9) we can estimate 


that the latent heat of special water is about 
150-50 cal/gr, [(0.3~0.1)-0.002]. 


(5) Lufluence of Temperature on the Strength 
of Clay 


In the previous chapter, that there was 
found difference in the quantity of special 
water between the undisturbed clay and the 
remolded clay was delivered. 

So, the author experimented if the dif- 
ference between the undisturbed clay and 
the remolded one could be interpreted ac- 
cording to the quantity of special water. 

As the temperature of the undisturbed clay 
increases, the contained special water turns 
into common water and displays the same 
effect as the increase of water content, so 
that the strength has to diminish. While in 
case of the remolded clay, however, it can 
be expected that the strength should rarely 
have any relation to temperature. 

By leading water in the pressure chamber 
of the triaxial shear equipment into the 
capacious thermostat, the clay specimen 
equipped in the triaxial shear equipment was 
made surrounded by water at some constant 
temperature and was keeping the same tem-. 
perature as that of water. About the clay 
specimen thus settled, stress controlled un- 
confined compression test was performed 
without giving any pressure to water in the 
pressure chamber. The size of the specimen 
and the rate of stress increase is quite the 
same as were refered to in 2.1. The thus 
obtained relation between temperature and 
the strength is, in case of the remolded clay, 
shown in Fig. 2.8 and in case of the un- 
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Fig. 2.8 Relation between strength and tem- 
perature of remolded clay (Taishobashi 208, 
16.m16-17.m06A, Water content 49.5%) 
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Fig. 2.9 Relation between strength and tem- 
perature of undisturbed clay (Taishobashi 
208, 16.16-17.06, Water content 49.5%) 


disturbed clay is shown in Fig. 2.9. From 
Fig. 2.9, we can observe that the strength of 
undisturbed clay lessens from 2.1kg/cm? to 
1.3kg/cm? when temperature increases from 
6°C to 36°C. On the other hand, it is also 
estimated from Fig. 2.5 that the electric 
resistance of clay specimen diminishes from 
2600 to 1400 by the same amount of the 
increase of temperature. As was elucidated 
in 2.4, we may suspect that, out of the whole 
decrease of electric resistance of the un- 
disturbed clay due to the rise of temperature, 
one fifth of its decrease is caused by special 
water tranforming into common water. There- 
fore, it can be considered that, in this case, 
the increase of common water is worth one 
fifth of (260-140)=1200, i.e. 240. However, 
that the increase of water content worth 240 
of the decrease of electric resistance is 4% is 
attained according to the relation between 
the water content and electric resistance of 
Fig. 2.3. The estimation thus indirectly 
taken that, by 4% of the increase of water 
content, the strength changes from 2.1kg/cm? 
to 1.3kg/cem? is well coincident with the 
strength and was content relation taken from 
the direct measurement of the effect of water 
content upon the strength. 

And from the relation between the water 
content and strength we may perceive that, 


if at 25°C, in order to fall 1.6kg/cm? of 
strength of the undisturbed clay of 0.4kg/cm? 
of strength of the remolded clay at the same 
water content, it is necessary to increase 
18% of water content. This fact is equi- 
valent to that, as was related upon the 
electric resistance of Fig. 2.3, if the electric 
resistance and watercontent relation curve of 
the undisturbed clay is shifted towards an 
origin by 20% of water content, it lies upon 
the same relation curve of the remolded clay. 
Assuming that the special water is non- 
conductive body and comparing as for the 
quantity of common water alone, we can 
estimate from the hitherto calculation that 
as to both the strength and the electric 
resistance any difference does not appear 
between the undisturbed clay and the remolded 
clay. The specimen of same shape was made 
of the fully remolded sample of undisturbed 
clay, the electric resistance of which is Ru, 
with its water content unchanged. If its 
electric resistance is kz, the quantity Ws of 
common water in the undisturbed clay is 
expressed followingly ; 

Ws= We (R7/Ru)3/2 (2.4) 
where We is the water content measured by 
the ordinary method. As to the clay which 
was sampled from the extensive area, the 
author noticed that the strength was deter- 
mined uniquely with the common water 
content taken by subtracting the special 
water content from the water content. It is 
shown in Fig. 2.10. 
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Fig. 2.10 Relation between strength and 
common water content. 
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(6) Relation between Strength Recovery and 
Variation of Electric Resistance 


By the experiment so far, we apprehended 
that the strength of this clay was decided 
uniquely with the common water content 
obtained by subtracting the special water 
content from the water content measured by 
the ordinary method, and that according to 
the equation (2.4) the common water content 
was estimated from the values of electric 
resistance both in the undisturbed state and 
in the remolded state. Then, provided that 
the causes of the phenomenon of strength 
recovery of the remolded clay in resting are 
due to the increase of special water and to 
the decrease of common water during resting, 
the electric resistance of remolded clay must 
increase during the resting period. And 
that, in reality, the electric resistance, similar 
to the strength, increases in accordance with 
the proceeding of resting day is shown in Fig. 
2.2. Therefore, it may be supposed that the 
strength recovery due to resting should be 
caused by the increase of special water con- 
tent and by the decrease of common water 
content during the resting period. 


(7) Relation between Suction Force and 
Special Water Content 

It was corroborated that, as was afore- 

mentioned, the change of nature of clay 


either by remolding or not was owing to the 
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Fig. 2.11 Consolidation curve. 
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change of quantity of special water, and 
consequently to the change of common water 
content. In this chapter, the relation between 
such special water content and suction force 
was studied. 

Since consolidation pressure is considered 
theoretically as suction pressure, the relation 
between the consolidation pressure and water 
content can be regarded to be the same as 
the relation between the suction force and 
water contained. As shown in Fig. 2.11, the 
typical example of the relation between the 
void ratio and the consolidation pressure is 
expressed DEC in case of the remolded clay 
and ABC in case of the undisturbed clay. 
As the consolidation pressure becomes greater, 
or as the water content becomes diminished, 
the difference of void ratio between the un- 
disturbed clay and the remolded clay becomes 
smaller. This suggests the fact that the 
quantity of special water in the remolded 
clay increases, as consolidation proceeds. 
Decreasing of the sensitivity ratio in ac- 
cordance with proceeding of the consolidation 
is ascertained by the experiment.2) And 
it is shown as well that the strength ratio 
between the remolded clay and the clay 
immediately after remolded goes larger as the 
degree of liquidity index is higher,* which 
may be interpreted as follows. 

Between the strength S and the water 
content We, there is a relation, 

log S=a—B We 
in which @ and § are constants.“ Accord- 
ingly, if S changes to (S+4S) when We is 
diminished by We, we can get the following 
equation, 
log (S+4S)/S= 84 W 

The quantity of the decrease of common 
water in order that the special water may 
increase during resting is considered generally 
to be in proportion to the water content or 
to the liquidity index. So, we can be sure 
that the rate of increase of strength becomes 
larger as the liquidity index of the specimen 
is greater. In the experiment of micro- 
calorimeter which will be mentioned after- 

* Liquidity index=(Water Content-Plastic Limit)| 
(Liquid Limit-Plastic Limit) 
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wards, we see the fact that the difference of 
special water content between the undisturb- 
ed clay and the remolded one becomes small 
as the consolidation proceeds is no other 
than the fact that the special water content 
of the remolded clay increases whereas that 
of the undisturbed clay remains constant. 
The author corroborates in the following 
experiment that the remolded clay, when its 
water content decreases so that its special 
water content increases, approaches to the 
undisturbed clay. Though, in case of this 
remolded specimen, its water content is 
comparatively small, air does not enter into 
it because the evaporation speed of its soil 
moisture is slow. The fact that the spacimen 
contains no air can be corroborated by the 
measurement of its void ratio. Of the speci- 
mens at various water contents air-dried 


naturally in the atmosphere, the author 
measured the electric resistance and the ob- 
tained result is shown in Fig. 2.12. The 
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Fig. 2.12 Relation between water content 
and electric resistance of remolded clay 
which has low water content. 


measured electric resistance shows greater 
value than that obtained according to the 
equation (2.1) by assuming that the spacimen 
should absorb as much air as the volume of 
water decreased from the initial water con- 
tent, 46%. Since such a phenomenon means 
that the decrease of common water which 
holds electric conductivity is greater than 
that of water content ordinarily obtained, we 
ought to consider it as is owing to increasing 
of the special water content in accordance 
with decreasing of the water content. How- 
eve, if the electric conductivity of contained 


water suffers any change due to the variation 
of water content, the above-mentioned discus- 
sion can not be consisted. 

Hence, taking the relation the electric 
conductivity of contained water and the water 
content obtained by the experiment, the 
author confirmed the following fact. 

In suspension state where the water con- 
tent of clay is great and the fraction of clay 
is small, as has been well-known in the field 
of electric chemistry of liquid, the electric 
conductivity of liquid part is in proportion to 
root of the volume concentration of clay 
particle. But the author took the result of 
experiment that, if water content is below 
100%, the electric conductivity of liquid 
part might be considered as constant without 
having any relation to the water content. 


(8) Study on Specific Heat Capacity of Clay 
by Micro-calorimeter 

In the previous calorimeter-experiments 
much clay was used, but the quantity ob- 
tained by a thin wall sampler from a single 
location was small. Thus, the measurements 
of the specific heat capacity of clay at various 
water contents and the periods of rest were 
difficult because of the lack of much clay of 
one kind. To overcome this difficulty, a 


miniature calorimeter was adopted as shown 
Wi Big lo. 


® 


Fig. 2.13 Miniature Calorimeter 
(1) Vessel (2) Heater of manganine line 
(3) Sample (4) Rubber Sleeve (5) Water of 
120cc (6) Stilole resign (7) Vacuum bottle 
(8) Heater (9) Thermometer 
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The apparatus is the same as before used, 
but is of very small dimensions: the diameter 
of sample (3) is 3.4cm and the length is 8cm, 
the diameter of vessel (1) is 5cm and the 
height 25cm. The heater (2) attached to the 


particles dessicated at 110°C by using this 
micro-calorimeter. One example of the results 
of this experiment is shown in Fig. 2.14. The 
plotted values in Fig. 2.14 express the mean 
of three measured values. Clearly the specific 


bottom is made of manganine line in order to heat capacity of contained water increases 


avoid the variation of electric resistance 
which is apt to arise due to the change of 
temperature. The electric voltage given in 
the heater, once having been made constant 
by the stabilizer, passed through the trans- 
former to 6V. 

The specimen, covered with a _ rubber 
sleeve, was submerged in 120cc of water ) 
to ensure equal heating throughout. This 
miniature vessel was covered with a stirole 
resin, the outside of which was attached to 
a thermos.“ 

With such an adiabatic method, the flow of 
heat in the outside cannot be neglected if the 
temperature of sample rises 5°C higher than 
the room temperature because of the small 
heat capacity of the vessel and its contents. 
Consequently, another electric heater“) was 
inserted just in the inside of thermos and 
temperature inside the thermos was kept 2°C 
lower than the temperature of clay specimen, 
by controlling the voltage. 

In such a way, the specific heat capacity of 
the remolded clay for various resting periods 
could be taken. The specific heat capacity of 
water contained in the clay could be deter- 
mined by deducting the specific heat capacity 
of clay particles, which had been estimated 
from the specific heat experiment of clay 


Specific Heat Capacity 
of Contained Water 
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Fig. 2.14 Relation between specific heat 


capacity of contained water of remolded 
clay and resting period, 


with the period of rest. This indicates that 
the special water in the clay mentioned 
above increases and the remolded clay itself, 
as the time of resting proceeds, approaches 
to the undisturbed clay. 

Next, we measured the specific heat capacity 
of contained water of the remolded clay and 
the undisturbed one of various water contents. 
The result is shown in Fig. 2.15. The rela- 
tion between the specific heat of the contained 
water and the water content in the clay state 
was calculated followingly. 
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Fig. 2.15 Specific heat capacity of contained 
water in various water content of remolded 
and undisturbed clay. 
@ Undisturbed © Remolded 
Assuming that the specific heat capacity of 
contained water, at 70% of water content in 
the clay state, is 1.3, 1.4, 1.5 or 1.6, and that 
the water lost when the water content 
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decreases from 70% is common water, of 
which specific heat is 1.0, in case at any 
water content We the mean specific capacity 


of contained water is accordingly expressed 
as 


<0.75—(0.7— We)}/ We 


where 6 is 1.3, 1.4,1.50r1.6. The measured 
values were compared with the calculated 
curves, and the water lost in the undisturbed 
clay was estimated to be common water. In 
the remolded clay, however, the increase of 
specific heat capacity of contained water 
corresponding to a decrease of water content 
cannot be explained away only as a result of 
the loss of common water. This means that 
in the remolded clay the special water is 
increasing corresponding to decreasing of 
water content. 

This is also confirmed by the experiment 
that the strength variation of remolded clay 
the special water is increasing corresponding 
to decreasing of water content. 

This is also confirmed by the experiment 
that the strength variation of remolded clay 
by temperature increases with the proceeding 
of resting period and approaches to the be- 
haviour of undisturbed clay, as is shown in 
Fig. 2.16. 
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Fig. 2.16 Relation between strength and tem- 
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§3. Okuramura Clay of Yamagata 
Prefecture 
(1) Sensitivity 
Taken from the the part of creep type 
slide region, a sample of this clay might 


suffer considerable disturbance even in 
sistu state. Therefore, in the process of 
Sampling, it was cut off by sawing so as not 
to break its skelton and was conveyed 
wrapped with vinyl cloth. The author con- 
sidered the sample which had been rested in 
the laboratory for two years as the undisturb- 
ed clay. In reality, the experimental results 
of its strength recovery corroborate the ap- 
propriateness of his assumption. 

Many specimens were made from the fully 
remolded clay with its water content un- 
changed and were conserved in polyethilene 
sack. The shape of sample is cylinder, 2cm 
in diameter and 5cm in height. By the strain 
controlled unconfined compression test, the 
temporal variation of stress and strain was 
self-recorded. The rate of compression 
displacement is about 6.7 x 10-*cm/sec. From 
this record, the stress-strain relation curve 
can be obtained. One example is shown in 


Fig. 3.1. Since there is found some fluctua- 
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Strain 41/1 —> 
Fig. 3.1 Stress—Strain curves 
Numbers in the diagram indicate the resting 
days. 


tion in the water content of sample, the 
relation between the strength and resting 
day -was taken according to the respective 
classes of its water content. And the result 
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vUndisturbed 


Strength/2 > 


9g °e* 


0 20 1c ee ae 100 Dayp 
Resting Period —> 
Fig. 3.2 Recovery of strength 
Water content 36-40% © 
41-46% x 
46-50% ®@ 
51-55% © 


is shown in fig. 3.2. From it, we can per- 
ceive that the strength diminished by re- 
molding completely recovers in about two 
months and Tshebotarioff’s sensitivity ratio 
of the clay after resting is very small. 
(about 3). 


(2) Result of Creep Measurement 

In order to study the reason of sensitivity 
of this clay as delivered in 3.1, the measure- 
ments of void ratio, of the electric potential 
of particle surface double sheet were taken. 
The method of measurement is the same way 
as was performed in 2.2. As its consequence, 
it was observed that the void ratio and the 
electrical potential of particle surface double 
sheet are both alike almostly unchanged by 
remolding. 

Moreover, the electric resistance and the 
specific heat capacity were measured by the 
same method as the one mentioned in chapter 
2, and as well any change by remolding was 
not observed. 

Lastly, the author could explain considerably 
clearly the difference between the undisturb- 
ed clay and the remolded clay by the creep 
measurement. This clay, if given shear 
stress, occurs stationary flow after some 
period. (see fig. 3.3) 

Whereas the strain of common clay or 
metal is in proportion to the logarithm of 
time under the constant shear stress, the 
strain of this clay is, as is noticed in case of 
the frozen clay, ice or the remolded clay, in 
proportion to time. 


Strain > 


0 20 40 60 80 100 
Hour —~> 


Fig. 3.3 Creep diagram (Normal stress 1.59kg/cm?, 
Shear stress 0.318kg/cm?) 
oa: transient strain after shear stress is given 
ab: consolidation region 
be: steady flow region 
cd: transient strain after shear is removed 
def: retardation region 


As the applied stress is varied, the strain 
rate changes. So, arranging the results of 
experiment, the author obtained the following 
equation as the relation between the strain 
rate y and the stress ft: 


7=A exp (Br) 


As this creep occurs caused by the mole- 
cular rate and the flow process, by assuming 
that we may neglect the probability of the 
molecule removing from the lower position 
of potential to its higher position, A and B, 
both of which are constant with physical 
meanings, are expressed as follows: 

A=2in@®/h exp (—E/®) 

B=;is/2®@ 
where @ equals, KT, K is Boltzmann con- 
stant, T absolute temperature, 4 is average 
distance projected in the direction of stress 
between equilibrium positions, m is the 
number of such process per length, # Plank’s 
constant, E the free energy of activation for 
the process, and s is the section area of such 
process vertical to the direction of stress. 

Of the remolded clay and the undisturbed, 
clay, the experimentally obtained A and B 
are shown in table 3.1. 


As is recognized in table 3.1, B differs 
whether in the remolded or in the undisturb- 
edclay. Within B=s/2@, @ does not change 
by remolding. Accordingly, that which 


(3.1) 
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changed is 4s. That the surface area of clay 
particle will not be changed is presumable, 
and that the void ratio of clay will be 
almostly unchanged too is ascertained by the 
experiment. Therefore, the values of 2 and 
s until multiplied by direction cosine, i.e. 
the length of a latus of the three-dimentional 
net work made of the scale-like clay particles 
and the surface area of the scale-like particle 
are both unchanged by remolding. So, the 


fact that the value of 4s changes means that 
the cross-angle of scale-like clay particles 
suffers some change. 

If the scale-like particles of the undisturbed 
clay really form rectangular cross-angle, 
while those of the remolded clay crosses at 
random, the calculated ratio of B of both is 
well consistent with the values obtained by 
the experiment. 


Table 3.1 Values of A and B 


end | Remolded does i . Undisturbed 
eee ee tcent | OE | 28. | 0.54 1.08 1.62 2.13 
A | ae ree | ee Ve Om? cic rat eon bie Ze 
B | os ties on ae tehseayee IMO =e | ee i; x 10-5 1.3x10-5 are 
| 
(3) Re son of Sensitivity other in which 7 alone changes by remolding. 


According to Table 3.1, it can be observed 
that the value of A is different whether for 
the remolded clay or for the undisturbed 
one, being small in the latter case, and that, 
as the normal stress increases, it becomes 
smaller togather in both cases. What is 
possible to be the cause of the change of A 
is estimated as is no other than thé change 
of E or 4 from the equation (3.1). As to 
' what degree v or 4 can change, the author 
speculated followingly. 

By the measurement of electric resistance, 
we could ascertain that this clay is isotropic, 
so that we may consider s=1/n?. Accordingly, 
B=2/2n*®. 

Then, as to the fact that B differs either 
in the remolded sample or the undisturbed 
one, two cases can be considered; one in 
which only 4 changes by remolding and the 


Table 3.2 Values of 2./4,, assuming 4 
is variable and ” is constant. 


er 
| 


Normal | 0. 54ke/em? | 1.08—2.00 | 2.13—2.66 
Stress | 
WT. | Pesto 1.44) 


i 


The variations of both cases were obtained 
from Table 3.1, and the results are shown in 
Table 3.2 and Table 3.3, where n., n, and 4.., 
2, represent respectively m and 2 of the 
undisturbed specimen and the remolded speci- 
men. 


Table 3.3. Values of ,/n., assuming ” 
is variable and 2 is constant. 


Normal 
Stress | 


Nr | ten | 1.45) | 1.2(5) | 1.2 


| 
0.54kg/em? | 1.08—3.00 | 2.13—2.16 


Provided that the case of the undisturbed 
clay is standard and 4.2,.=1, the equation 
(3.1) is expressed as follows : 


log A=13.08—1.08 x 10°F (Sra) 


where @)=4.0x10-", and A=6.6x10-?’ are 
used. 

In calculating log A of the remolded clay 
for the two cases, one in which 4 changes in 
consequence of remolding and the other in 
which changes in consequence of remolding, 
by utilizing table 3.2 and Table 3.3 respec- 
tively, some correction shown in Table 3.4 is 
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Table 3.4 Correction value of 13.08 


EEE 


needed to the value 13.08 in the equation 
(3.2) 


Normal | ~200 ~2.66 ay: ‘ 
Stress Deere Les wine By utilizing both table 3.4 and the equation 
Rote ee 8) eee, uaa Ot ae (3.2), we can obtain the value of E from the 
a; variable 20.33. | /=0.90 lex Osl6 , 
nm: constant | : | ; | value of Ain Table 3.1. The result taken is 
| i | . 
ca reaans ea ibe’ Pes shown in Table 3.5. 
A: constant i 
Table 3.5 Values of E 
a 
Kind | Remolded | Undisturbed 
] a ie i | | 
Normal Stress aa 1.08-2.00 | 2.13-2.66 | 0.54 | 1.06 1.62 | Zee} 
A: variable 4 gg.19-12/ 1.96x10-" 2.02, x 10-# 
m: constant | 
= ! 1.97 x 10-12 | 2.00 x 10-* | 2.02 x 10-14/2.06 x 10-2 


m: variable 
A: constant 


The value of & obtained as for the creep 
would represent the quality of micro-particle, 
and the author will corroborate later in 3.4 
that the microparticle plays a main role in 
the strength of this clay. Therefore, con- 
sidering that the strength is defined with E 
and that it is in proportion to exp (£/@), we 
can express sensitivity ratio Sr followingly : 

Sr=exp [(Eu—Er)/®] 
in which Lu and Er are the values of E of 
both the undisturbed specimen and _ the 
remolded specimen respectively. The cal- 


culated values of Sv by using the values of 
Table 3.5 are shown in Table 3.6. 


Table 3.6 Values of Sp 


| 
Normal Stress pease 1.08—2.00 2.13 Measured 


| | Values 
A: variable | | 
ie: constant | 9.3 | 3.5 | 2.4 
Sy,— —— bs 7 a 3 
m: variable 
a: constant 9" | 1.6 eas) 


The result that the volues of S, becomes 
smaller in accordance with the increase of 
normal stress is well consistent with the 
other reports. @ @ 

From the calculation hitherto, the author 
speculates as to the cause of sensitivity as 
follows. 


| 1,92, 10-74 1.99 x 10-12 | 2.04, x 10-# 


chalet 


The fact that the strength of this clay 
becomes different either by remolding or not, 
even at the same water content, or with the 
same void ratio, may be interpreted due to 
the reason that in the undisturbed state the 
arrangements of the scale-like clay particles 
are in good order and the removal of clay 
particles are performed under the condition 
of much activation energy, whereas in the 
remolded state the arrangements of the scale- 
like particles are at random, so that they can 
turn freely and their removals are possible 


even in the condition of less activation 
energy. (see Fig. 3.4) 


Fig. 3.4 Three dimensional network of clay 
(after Tan-tyong-kye’s paper). 


(4) Strength Increase in Resting of the Paste 
made of Microparticles of Clay alone 

As was discussed so far, the sensitivity 

was explained according to the activation 
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energy obtained from the creep phenomenon. 
However, it might be possible that the larger 
particle of clay should play the main role in 
the strength, while the microparticle per- 
forms the main role in the creep phenomenon. 
And, according to Casagrande’s theory it is 
also probable that the clay without silt should 
not give rise to bond clay so that it should 
hold no sensitivity. 

Then, the author, utilizing the wet me- 
chanical analysis, produced the paste consisted 
of the microparticles, the particle size of 
which is below ly, which were made by 
precipitating the larger clay particles and 
then by evaporating the remnant fluid. And 
he examined the aspect of the strength 
increasing of this paste during resting. The 
examined aspect that the strength of this 
clay reaches saturation in some two months 
of resting is similar to that of strength 
recovery of the original clay. (see Fig. 3.5) 
This proves nothing but the fact that the 
microparticle performs an important roles as 
well in the strength recovery of the original 
clay. 
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Days 
Resting Day > 
Fig. 3.5 Strength recovery of samples which 
are constituted with microparticle 


Therefore, it may be appropriate that we 
should explain the sensitivity by the activa- 
tion energy obtained from the creep pheno- 
menon. 

And the author studied experimentally, too, 
that the strength recovery of the clay paste 


mentioned is not caused by the change of 
void ratio, and showed the result in Fig. 3.6. 
According to this figure, the void ratio is 
regarded as constant without any relation to 


Void Ratio > 


a 10 10 10 
Resting Day —~> 
Fig. 3.6 Void ratio-resting day relation (Water 
content 66-71%) of the paste which consists 
of microparticles smaller than 1. 


Dox 


resting day. Then, in Fig. 3.7 and in Fig. 3.8 
are shown the water content and the void 
ratio of each clay paste respectively classified 
by the particle size by the wet mechanical 
analysis, when put in the atmosphere and 
the weight of which keeps equilibrium. Since 
the condition in which the weight of clay is 
in state of equilibrium in the atmosphere is 
below the shrinkage limit, its void ratio 
ought to increase generally as the particle 
size becomes small. 


t Hl 

g g 

8 a 

2 Sa 

S wo 10P . 

wu ° 09 ° 

= = 0B . 

= 01 é 
10 (oaonie ne 


AB RCSD VE 


Fig. 3.7 Relation between 
water content, which is 


Fig. 3.8 Relation between 
critical void ratio and 


balanced with atmos- PALvIClem SIZE essen, 
pheric vapour pressure, B-0.74, C-0.65', D- 
and particle size A-1.1z, 0.554, E-0.42u. 
B-0.784, C-0.654, D- 

0.554, E-0.42u. 


But Fig. 3.7 represents that the void ratio 


decreafes as the particle size becomes smaller. 
Such a fact is just consistent with the result 


consisted merely of the microparticles above- estimated due to the assumption that this 
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clay paste forms three dimentional network 
of the scale-like clay particles of which 
thickness is uniform. And the fact that, 
as shown in Fig. 3.8, the water content 
increases as the particle size becomes small 
proves that the surface area of clay particle 
per unit volume of clay increases as the 
particle size becomes smaller, since the 
thickness of absorbed water of clay particle 
which keeps equilibrium with the air can be 
regarded as uniform. 


§4. Study on Shrinkage Process 


In the previous chapter 3, we presumed 
that the strength of consolidation test was 
decided by the activation energy with which 
a clay particle passes by the neighbouring 
clay particle. This fact can be ascertained 
by the shrinkage process of clay when the 
soil moisture of which evaporates. Namely 
it means that, as the strength of clay is 
larger, its activation energy becomes greater 
and the shrinkage rate accordingly diminishes. 
Therefore, even when the evaporation rate 
of soil moisture remains constant, the clay 
may differ in the shrinkage rate, in the way 
of entering of air and in the process of 
variation of void ratio due to the difference 
of strength. 

Then, the following experiment was _per- 
formed. Preparing a lot of undisturbed 
specimens and remolded specimens, we desic- 
cated them naturally in the atmosphere. By 
varying the desiccation period, we can obtain 
the specimens in various stages of the process 
of desiccation. The water content and the 
void ratio of these specimens in every stages 
were measured and shown in Fig. 4.1 and Fig. 
4.2. As was expected, of Osaka alluvial clay, 
the void ratio of undisturbed specimen is 
greater than that of remolded specimen at 
the same water content. Of Okuramura clay 
of Yamagata prefecture, however, the void 
ratio of remolded specimen proves a little 
greater than that of undisturbed specimen at 
the same water content. Whereas the strength 
of Osaka alluvial clay, for instance when its 
water content is 40%, is about 6.0kg/cm?, 
that of Okuramura clay, when its water 


Void Ratio > 
a , 


50 % 


10 20 30 40 
Water Content —~> 


Fig. 4.1 Drying process of Osaka alluvial clay 
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| Natural dry a O 
| Slow dry a EN 
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Fig. 4.2 Drying process of Okuramura clay 


| | Undisturbed | Remold 
leg | 

| Natural clay | cs) | O 
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content is the same as the former, is only 
about 0.8kg/cm? and its activation energy is 
small too. So, in case of Okuramura clay, 
since the shrinkage rate can follow the 
evaporation rate of the soil moisture even in 
the undisturbed state, the rate of air entering 
into the clay in the process of desiccation 
does not differ whether for the remolded 
specimen or for the undisturbed specimen. 
Consequently, it may be considered that the 
void ratio of remolded specimen becomes 
greater by the air content entered by remold- 
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ing than that of undisturbed specimen. To 
corroborate this idea furthermore, the author 
put the Osaka clay in the saturated-vapour 
chamber in the desiccator which contained 
water in stead of the drying medicine in the 
bottom, and desiccated it slowly, while Okura- 
mura clay was desiccated rapidly in the 
drying owen. In consequence, as shown in 
Fig. 4.1, in the former case there scarcely 
appears any difference between the remolded 
specimen and the undisturbed specimen in 
their void ratio, but the void ratio of un- 
disturbed specimen becomes greater again 
than that of the remolded one, if below a 
certain water content. It can be interpreted 
that, as the water content diminishes, both 
the strength and the activation energy become 
greater together and the shrinkage rate of 
the undisturbed specimen correspondingly 
decreases too much to follow the evaporation 
rate of such a degree. 

On the other hand, in case of Okuramura 
clay, the void ratio of undisturbed specimen 
becomes surely greater. This fact proves 
that the shrinkage rate is exceeded by the 
evaporation ratio of soil moisture. Though 
the shrinkage rate can accompany the eva- 
poration rate in the exterior part of remolded 
specimen, that rate is so large that the 
shrinkage in the interior part of the specimen 
cannot follow. As its consequence, many 
cracks are produced. 

Then the author made the dimension of 
the specimen smaller or flatened, or keeping 
the door of the drying oven open to check 
the evaporation rate and yet such cracks 
couldn’t be prevented. So, by taking the 
inner, more-or-less uniform shrinking part of 
specimen, the void ratio and the water con- 
tent were measured. Though every aspects 
of the change of water content could not be 
grappled, it can yet be observable that the 
void ratio of remolded specimen becomes 
smaller, compared with that of the undisturbed 
specimen at the same water content. (see Fig. 
4.2) 


§ 5. 


By taking as sample both Osaka alluvial 
clay which is called Kaolinite and the clay 
gathered from Yamagata prefecture, Okura- 
mura land creep region, the author studied 
as to the sensitivity. As the result, we 
could comprehend that it was explained due 
to the special water in case of the former 
clay while in case of the latter clay it was 
explained according to the difference of the 
arrangement of the scale-like particles. 


Conclusion 
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Gravity Survey on Kuttyaro Caldera Lake. 


By 


Izumi YoKoYAMA 
Geophysical Institute, Hokkaido University. 


Abstract 


A gravimeter survey was made on the frozen surface of the gigantic caldera lake, 
Kuttyaro. A Gravity low, reaching —46 milligals, was found, with iso-anomaly contour 
lines conforming the caldera shape. A layer of coarse pumice thickening toward the centre 
is required to account for the observed anomalies. This is compatible with volcanological 
inferences. Gravity anomaly distributions at several other calderas are discussed also. 


$1. Introduction 20km in diameter, this is one of the most 
; gigantic calderas on the earth. A small 
ae se Severo, located in the eastern part jsjand at the centre of the lake and parasitic 
of Hokkaido, (Fig. 1), Japan, represents the yojcanoes in the eastern half of the caldera 
western half of a circular caldera which was erupted after its formation. Pumice which 
formed in Quaternary age. Measuring about was ejected at the time of the catastrophe of 
the caldera-formation is found widely spread 
around the lake, and in fact, even at the 
shore of the Okhotsk Sea as distant as 50km 
from the volcano. As for geophysical pheno- 
mena in this district, it is known that a 
moderate earthquake (M=6.0) occurred in 
May 1938 being accompanied by small crustal 
deformations. The hypocentre of this earth- 
quake was located at the south-east corner of 
the lake and 20km in depth, although no 
more detailed studies were prosecuted. Being 
interested in the underground structure of 
this typical caldera, the present writer wished 
to get gravity anomaly values on the lake. 
At present, gravity measurement on water 
is of course very difficult. But fortunately 
the lake water freezes hard in winter, so that 
a Worden gravimeter could be used on the 
surface of the frozen lake. A gravity survey 
on the ice surface was carried out in March 
1958. The Bouguer anomalies on and near 
the lake are shown in Fig. 2, which also 
contains the anomalies obtained by the Geo- 
graphical Survey Institute by another series 
of gravimeter survey along the route of 
precise levels along the eastern margin of the 
Akan calderas. Contour-interval is 200m. caldera. Since our gravity points are situated 
K: Lake Kuttyaro, A: Lake Akan, on the level surface of the lake, the free air 
M: Lake Masyu. reduction and the Bouguer reduction are 


=== Okhotsk 2 Sea 


Fig. 1. Topographic sketch-map of Kuttyaro and 
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uniform throughout the lake. This fact 
simplifies the calculation very much. Topo- 
graphy corrections are very small and can 
safely be ignored in the following discussion. 
A conspicuous fact to be noted in the figure 
is that a relatively low anomaly which 
reaches about —46 milligals, is found, the 
lines of equal anomaly being concentric with 
the centre of the caldera. Comparing the 
distribution of Bouguer anomalies with that 
of depth of the lake, shown in Fig. 3, no 


Om 


ss 
Fig. 2. Distribution of Bouguer anomalies on the 
Kuttyaro caldera (in milligal), 


relation is found between them, and therefore 
the Bouguer anomalies observed may be 
regarded as due to the subterranean structure 
of a larger scale. 

Here, the writer intends to interpret the 
Bouguer anomalies found on this Kuttyaro 
caldera and to discuss its subterranean 
structure, together with those of other cal- 
deras from the standpoint of volcanology. 


Fig. 3. Distribution of the depth of Lake 
Kuttyaro (in meter) after H. TANAKADATE. 


§2. Subterranean structure of the Kuttyaro 
caldera from the gravity anomalies 


The gravity potential due to a certain 
anomalous subterranean mass distribution 
can be expressed in the following Bessel 
Fourier series form: 


i) 


4U(7,0,2)=3o| Oak J.Akr) exp (—kz) an nd 
(1) 


using customary notations. From this ex- 
pression for the potential, we get the gravity 


anomaly at the earth’s surface z=0 as 
follows : 


oo 


Ag(r,0)= SY |BanSu(hn dk 


ens 


(2) 
sin MS 
Assuming that the distribution of gravity 
anomaly on the Kuttyaro caldera shown in 
Fig. 2 is nearly concentric, (2) can be con- 
verted into a simpler form : 

a t: 

$e \4 g(r, 0) d0= [Bate (er) dk. 
The left-hand side is the average of 4g(r, A) 
taken over the whole azimuth @ along a 
circle r=r, drawn around the centre. Denot- 
ing this average by 4g(r), we get 
4g(0)= [Bus dk, 


0 
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and 
TE Gye |B. i tER db. 


By the method of least squares, the five 
coefficients Bu, Bo, --, Bo; in the above ex- 
pression are determined over the range from 
r=0 to r=10km in Fig. 4 by averaging the 
gravity values in eight directions with 30° 


4g 
meal 
60 
40 
<a 
5) 5 To km 


Fig. 4. Radial distribution of Bouguer anoma- 
lies on the Kuttyaro caldera. 


intervals around the centre of the caldera. 
The outermost circle is so large in radius 
that the gravity anomaly on it becomes 
negligibly small. Here, the centre of the 
caldera is taken at the point where the 
contours of the Bouguer anomaly appear to 
converge in Fig. 2 and this point coincides well 
with the geometrical centre of the caldera-rim. 
~The coefficients thus deduced are as follows; 


[ERE BoA Sie): Bou=-2.3, 
Bos=2:1. (milligal) 
In the present case, the gravity anomaly at 
the depth z=z is given by 
Ager, 2)= | Bo Jk) exp (—kz) dk. 


0 


(3) 


Using the above coefficients Bu, Bo, - -, 
B.;, the approximate distribution of gravity 
anomalies at various depths are calculated as 
shown in Fig. 5 until the series ceases to be 
convergent. The higher harmonics increase 
in prominence as the anomaly producing 
mass is approached, since the value of k is 
larger for the higher harmonics than for the 
lower, but they are unlikely to exceed 
materially the lower harmonics in amplitude 


Fig. 5. Gravity anomalies at the 
various depths. 


unless the mass is reached. The depth at 
which the above integral ceases to converge 
does not affirm that there are no sources 
down to that depth, but rather implies that it 
is more probable that the sources are to be 
found at that depth and greater. In Fig. 5, 
the curve of the field at a depth of 2km 
becomes much steeper than the original one, 
varying by about 110 mgals within a hori- 
zontal distance of 10km. This remarkable 
variation in anomaly may indicate that the 
anomaly producing mass is distributed at a 
comparatively shallow place from the horizon 
of the depth 2km with a probable density 
contrast. If a greater depth is assumed in- 
conceivably high density contrasts are re- 
quired to reproduce the field. 

In the above discussion, we reached the 
conclusion that the gravity anomalies observed 
on the Kuttyaro caldera are due to masses 
that are shallow, relative to the least horizon- 
tal distance in which the anomaly changes 
in its whole amplitude. The absolute value 
of the depth is put out of the question 
because, anyhow, the problem is clearly in- 
determinate. In this case an excellent and 
simple approximation to the mass distribution 
is given by 
MN (4) 
2nG 
At the centre of the caldera 7=0, 4g(0) 
amounts to -46 mgals. Then 2(0) is ob- 
tained as -1.1x10* gr./em. If we assume the 
density contrast to be 0.3 gr./cm*, the thick- 


m(r) = 
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ness of the subterranean mass at the centre 
becomes about 4km. Similarly, the distribu- 
tion of the coarse material which is respon- 
sible for the negative gravity anomalies is 
obtained as shown in Fig. 6. The selection 


Lake Kuttyaro 


2 Km 
5, 


Kook 


Fig. 6. Distribution of the coarse material at 
the Kuttyaro caldera. 


of the value of density contrast 0.3gr./cm* 
admittedly holds some ambiguity, though the 
value around this may be a probable one 
corresponding to go=2.0gr./em* and gi=1.7 
gr./cm* in Fig. 6. 

On the other hand, a cylindrical structure 
of which the axis extends vertically with a 
low density is also conceivable beneath the 
lake for approximately interpreting the ob- 
served gravity anomalies. In this model, the 
boundary surface of the cylinder needs to be 
situated far inward from the present caldera- 
rim: The writer believes that the subter- 
ranean structure as shown in Fig. 6 is the 
most probable, because the above-mentioned 
cylindrical structure is to be accompanied by 
a discontinuity in the profile of the gravity 
field but the actual results do not show any 
such discontinuities with in the accuracy of 
the survey. In other words, the caldera-rim 
itself does not represent the vertically deep- 


rooted fault from the view-point of density 
contrast. 


§3. Mass-defect caused by the caldera- 
formation 


The total anomalous mass or mass-defect is 
given by the theorem of GAUSS which 
relates the integrated anomally over a hori- 
zontal plane to the limit of the detectable 
gravity anomaly in the following form: 


00 00 


| [ee y) dxdy 


co — 00 


M= 

2nG 6) 
where G is the gravitational constant. Thus 
the total anomalous mass M can uniquely be 


determined from the gravity values observed 


at the surface while its distribution is in- 
determinate. If the discussion is confined to 
the caldera-formation, the distribution of 
gravity anomaly on the Kuttyaro caldera 
may safely assumed to be symmetric around 
its centre including unobserved areas where 
the parasitic cones erupted after the caldera- 
formation. The total mass-defect for this 
caldera is calculated by (5) using the radial 
distribution of gravity anomaly shown in 
Fig. 4 and is found to be 


M=7.8 x 10'° tons. 


For reference, another estimation is possible 
on the basis of the subterranean structure 
shown in Fig. 6. The volume of low density 
material in Fig. 6 is roughly 3x10''m* and 
if multiplied by the density-contrast 0.3gr./cm', 
the total mass-defect reaches 9x10" tons. 
The change in the density-contrast does not 
affect this estimation very much. 

The above values of the total mass-defect 
beneath the caldera should be almost equal 
to or less than that of the total material 
ejected from the caldera, because some parts 
of ejecta might have their origin at the 
deeper places than the region which was 
discussed above. In fact, Y Katsur (1958) 
has studied the distribution of ejected pumice 
around the caldera and estimated their total 
volume at 10''m® in the order of magnitude. 
Considering that the density of pumice is 
approximately 1.5~2.0gr./cm', the mass-defect 
found by the gravity survey is compatible 
with the total mass of ejected pumice and is 
therefore evidently due to the pumice-erup- 


tions accompanied by the formation of the 
caldera. 


This volume of ejecta from the Kuttyaro 
caldera, which is 10''m®‘, is a match for that 
in the greatest volcanic eruption in historical 
times, namely, the 1815 eruption of Tambora 
in Indonesia although the description of this 
eruption does not afford a precise estimation 
of its ejecta. 


§ 4. 
The mechanism by which enormous amount 


of pumice was ejected and by which the 
Kuttyaro caldera was formed, cannot definitely 
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be understood solely by the results of the 
gravity survey, but the eruption, accompanied 
by the accumulation of pumice, may clearly 
be said to have been so violent, that the 
region from the surface to a depth 3~4km, 
was completely devastated. In order to 
gain a clear understanding of the caldera- 
formation and the subsequent volcanic activi- 
ties, the following few references may be 
interesting. 

H. Kuno (1953) classified Japanese calderas 
as of Krakatoa- and Glencoe-type of H. 
Wiuuiams (1941) according to their causes of 
formation. According to his discussion, in 
the calderas of Krakatoa-type, the post- 
caldera lavas comprise a rock series formed 
through contamination of the same magma 
by granitic material, whereas in those of 
Glencoe-type, the post-caldera lavas represent 
an advanced stage of crystallization of the 
pre-caldera lavas. An example of the former 
is Kuttyaro caldera, which was devastated to 
a depth 3~4km during the eruption of a 
tremendous amount of pumice as mentioned 
before, and that of the latter is the Mihara 
caldera, Oosima Volcano, which was formed 
by the subsidence of a massive mass of rocks 
separated from the surroundings by a com- 
paratively simple system of ring fissures, 
without a pumice eruption. From this geo- 
logical standpoint, low’ Bouguer anomalies 
should be observed on the calderas of Kraka- 
toa-type and the high Bouguer anomalies on 
the calderas of Glencoe-type. Actually, the 
result of a gravity survey carried out on 
Volcano Mihara by the present author and H. 
Tazima (1957) shows the high anomaly 
amounting to about 15 milligals. In addition 
to the above example, C. Tsusor and his 
colleagues (1954) observed a low anomaly 
reaching about 20 milligals on the Aso caldera 
in Kytisyi which is of Krakatoa-type. On 
Lake Akan which is a part of the Akan cal- 
dera, Krakatoa-type, a gravity anomaly 
showing the same tendency as the Kuttyaro 
caldera was also observed in March 1958. 
These results of observations, though rather 
few in number, seem to justify the said 


geological classification concerning the forma- 
tion of the calderas, but the mechanism of 
subsidence or collapse is another problem. 
Although several theories have been put 
forward by geologists in regard to this 
problem, the present writer does not consider 
them to be well established, especially from 
the geophysical standpoint. 

If the subterranean structures of many 
calderas in the world are studied by geophy- 
sical methods, the old and much-disputed 
questions as to the caldera-formation will be 
revived and a valuable contribution to physics 
of the earth-crust may be expected. At the 
same time, physical properties of rocks in 
the deep and the depth of the magma- 
reservoir need be studied in relation to the 
caldera-formation. 

In conclusion, the writer wishes to express 
his sincere thanks to Prof. C. Tsusor who 
encouraged the writer throughout the course 
of this study and to Messrs. H. Tasma, T. 
Murase and I. Onpa for their kind coopera- 
tion in the field observations. His sincere 
thanks are also due to the Director of the 
Earthquake Research Institute who kindly 
lent the Worden gravimeter to him. 
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Ray-theoretical Construction of Dispersive RAYLEIGH Waves 


By 
Kyozi TazIMp 
Department of Geophysics, Faculty of Science, Hokkaido University, Sapporo. 


Summary 


Many studies have been published on dispersive RAYLEIGH waves. 


However they 


were mainly confined within either dispersion-curves or ratios of vertical displacement to 


horizontal one where POISSON’s ratio of the layer was usually assumed to be 0.25. 


The 


present author has attempted to carry out more general studies than before on absolute 


amplitudes for various values of POISSON’s ratio. 


In order to investigate the superficial 


waves near the surface of the earth, scholars must undertake these general studies, though 
theoretical treatments will be considerably troublesome. 
For this purpose, displacement-potentials have been expressed with reflecting coefficients 
in this paper where a useful operator has been derived from the characteristic equation. 
The characteristic equation in general has been expressed by a hyperbola. 


§ 1. General expressions of multiply reflected 
rays. 

“Solid layer over solid half space” (Ewinc 
et al., 1957) will be studied here. Displace- 
ment-potentials (Lams, 1903) of P and S-waves 
in each layer can be expressed respectively 
as follows: 

i= Arel1? + Bye-*1? dh) = Cye"Pi? + De-#8 1? ; 
(1.1) 
@2= Bre *#2? ; y= Dye *P 2? A (1.2) 
where the common factor of exp {d(wt—E&z)} 
is omitted and the original ray is excluded. 

In the above equations notations are em- 
ployed as follows: 

0;; density, 
4;, 45; LAME’s elastic constants, 
w; circular frequency, 
&; circular wave-number in the direction 
of x which is taken as horizon, 


P=hP—ayp=ke—B;?, hy=0/Up; , “Rj =00/03; , 
Vpj={(As+2u5)/05}? ,  Vs5= (43/05)? , 
j; the number of the layer counted from 
the most superficial layer. 
Recalling “SOMMERFELD’s radiation condition” 
(SOMMERFELD, 1949) to mind, the writer chooses 
the signs of the wave-number as next: 


£=E-iE , as=aj—ia;, Bi=Bj—18) (1.3) 


where &, E, tee Bs are all positive and real. 


Ai, Bi, +++, Dz in (1.1) and (1.2) are arbitrary 
constants which are to be decided by the 
boundary conditions. It is very troublesome 
to decide these six unknown factors from six 
simultaneous equations. Therefore another 
method will be attempted here. 

Reflecting coefficients as to displacement- 
potentials on the upper and the lower bound- 
aries of the superficial layer will be indicated 
by 

(Al 6, C,D) and (A’,. 2 CoD is) 
respectively in the order of PP, PS, SP and 
SS-reflections. 

Considering the line source of P-wave at F 

in Fig. 1, the P-ray going downwards from 


aaN 


Fig. 1. Space-coordinates under consideration. 


the origin will be reflected once on the lower 
boundary. Then the secondary P and S-rays, 
indicated by P:’ and Si’, will be generated. 
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The ray having the notation of ’ is in agree- the ray of P’-type must be proportional, as 
ment with the ray going upwards, correspond. it was, to exp(éaiz) and the ray of S’-type 
ing to the first term of (1.1). By the reflec- to exp (2/12). Thus P’ and S’-rays twice re- 
tions of P,’ and Si’ on the upper and the flected on the lower boundary can be totaled 
lower boundaries, number of rays will again as P,’ and Sy’. 

be doubled at each reflection. In this case, Taking account of phase-lags, (P2’, S:’) and 
(Py’, Si’) can be related by the next equations. 


Py = PY {AA e-#%14 + BC e-t PF} 4 Sy et BY2{ DC et PDE + AC’e 14} : in a | 
So! = Py e-*@1-8D4{ A B’e4@1 *8DE + BD’ e817} + Sy {DD e-*P1¥ + CBee titbaa} 
Putting 
a=AA’e-#%1F + BC’e-Mt BDz f b= DC e*@ 18D ELC A’e-2te18 k | q 6) 
c= AB ce tP 44 BD e-*A4 , d= DD/e-*F1¥+CBe ts BE , 


(1.5) will be extended to more general forms: 


Py! e-**1? = aPy/e-!*17 +S e781 — 
Psi e-*%1? = a Po’ e~'1? +S 0/ e-*F , 
rs 
Pp 1e7!"1? = aP np! e7'*1? 4S n/ cP 
(a Pr! +P. Piet aa x Py e7%12 +H > Sn’e7iB , (n=1, 2, 3, -+-c0), (1.7) 


in which either P,/,,<Pi’ or P,,,<>) Pr’ when n becomes large. Thus one has 
Pes < > a= 


Neglecting the left hand member in comparison with the right, (1.7) will be reduced to 


(a—1)> Pr’ )e*17 + B(D) Sn Je Fi? = — Pye , (1.8) 


In the same manner as the above, one gains, from the second relation of (1.5), 


CCX) Pn’ )e*1* + (d—1)(3) Sn’)etF1# = —Sy’e*Bas (1.9) 


Letting (1.8) and (1.9) be simultaneous equations, one obtains 


Putting, then, the P-ray going downwards from the line source as 
exp {—iai(z—E)} , Gla} 
one sees 
Py =e 18-2) A’e-ty (2) and Sy =e-t#1( 4-2) B’ e181 H—2) | 
Petr = A’e'Me-*%H and Sy/e-#it= Belt 8e-t tBu (112) 


On the contrary, putting the P-ray going upwards from the origin as 
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exp {iai(z—E)} , (URLS) 

one sees 
Py’e**1#=ae**1" and = Sy/e~*P1#=ce-!*14 (OL 


Because (1.11) and (1.13) are generated from the origin at the same time, P,’ and S,’ must 
be expressed by the sum of (1.12) and (1.14): 


Pye7t4127= A’ e-2 aA eiay8 + qe7te148 2 S1/e7B12 = B’ et B YH gle FE 4 cg-iayE : lca) 


Putting (1.15) into (1.10), one gets 


‘Oss Pao? =[{bB’ e184 —(d—1) A/e-#1 4} eit #4 (he—ad-+-a)e-*%1#/M 
(1.16) 
and (1 Sno"? =[{c A e214 — (a—1) Bet BYE} git 8 4 ¢g-t018)/ 
where 
M=(a—1)(d—1)—bce Cala 


and the right hand members ought to coincide respectively with A: and C: in (1.1). 

The ray going downwards in the superficial layer can be easily calculated from the rays 
going upwards which have already been obtained. Employing notations of (Pi, Si), (P2, S2), «+: 
for rays going downwards, one has the following relations: 

Po =e7**18 Ae-i*12 
Py= Py Aew™*12 + Sy Cet 14 Bidz ; 
Ps=P2 Ae-*%12 + So Ce-t@1t Bz : 
Pr= Pr Aer S, Ce G1 Ewes 
{Co Pn) + Pojet*v= Afe 19 -+(>) Pn’ )e-*17} + CCS Sn’ e741 . (1.18) 
In the same manner as the above, one obtains 
oe Sn) + Soper = Ble" +( 3) Poe} 3 DD Sn ee (1.19) 
Putting (1.16) into the right hand members of (1.18) and (1.19), one gets 
{(>) Pn) + Po}et*1?=[{bA B’e~ 114 8D 4—(d—1) A A’e-2!*14 4 cCA ‘e214 
—(a—1)CB’e~ 14 FD) 4} e419 4 {¢C—(d—1) A} e719 ]/M 
and {(>) Sn)+ So} e*1* =[{B BB’ e-*1*8 2 —(d—1)BA’e-#%12 4+ cDA’e- #14 
—(a—1)DB’ et BD) e418 4 {cD—(d—1) Bye?" ]/M 


(1.20) 


in which the right hand members ought to coincide respectively with Bi, and D; in (1.1). 
gz and #» in (1.2), if need be, can be calculated from (e!%1"+ Bie?" and D,e-*Fi*, multi- 
plying the refracting coefficients on the lower boundary, 
Reflecting coefficients A, B, C and D are easily expressed as follows: 


Cer Care 


ntl B (EY GF 


(1.21) 
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and they have the relation: 
AD—BC=1. (1.22) 
Using (1.6) and (1.22), one obtains the next expressions from (1.16), (1.17) and CEAD)E 
Ar=[{A/e-2"1¥ + D(B’C’ — A’ D’ e+ B DA} gta 8 
+{AA’e-%#%14 + BC’e*t BDH + (B’C’— A’ D’ je *#@14B 7} e191 9) M , 
Bi=[{AA’e721¥ + CB’ eB # + (B’C’— A’ D’ e241 B D2} gia 8 
+(A—D’e-2!812)e-'%1"]/M , (1.23) 
Cy=[(B’ e618 2 = BBC SA’ Deh Gr Pa ei" 
+{AB’e-*@1*8DF + BD e-*F12 }e-'%17/M * 
Di=[{BA’e72!%14 + DB’ et 144} i184 { B+ B’e-1+B DH eit B)/ MY , 
in which 


M=1—{AA’e-¥ 912 + DD e814 + (B’C’— A’ De # G18 8 + (BC + B’Che trt8p2} | (1,24) 


§ 2. A useful operator derived from the characteristic equation of superficial waves. 


If the right hand side in (1.24) is zero, € and w must be related in general by the equation, 


ME, o)=0 (2.1) 
which is called the characteristic equation of dispersive RAYLEIGH waves. 
In this case, one sees that 
dw 0M 0M 
Ge as ae ens) ee 


where U means group-velocity. Denoting phase-velocity by c, the next relation may be 
obtained from (2.2): 


(Tne) Hine ol Go tae ME > ae 
in which one may find that the operator, 
o=( a a | (2.4) 
has two remarkable properties: 
Was, Bs)=a;,B; and Mayj/E, B/E\=0. (2.5) 
On LovE waves, 
FE, wo) =1— Ke-2?14=0 (2.6) 


must hold (Taztmp, 1957). Now, because 


Ons=nj and Xnj/E)=0, (2.7) 


and the reflecting coefficient, K, is composed of the ratio of 92/2.=(H2/E)/(m/E) alone, the next 
calculation will be very easily executed: 
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OF E, o)= —0 Ke 14 = — Ke 14 = 2insHKe-214 (2.8) 
Then (2.6) must hold in (2.8), so one has from (Zed) 


1 1 2t 
a7 )=- Fabia) 


mi (E, SOTA =H) : (2.9) 


which latter is the amplitude-function of Love waves, being the same as that already obtained 
by the other troublesome method (Tazime, 1957). 


§ 3. Short notes on the characteristic equation. 


Exchanging reflecting coefficients as to displacement-potentials by those as to vertical dis- 
placements, one may find that (2.1) with (1.24) is similar to the expression that was found 
by Toustoy and Eucene (1953). 

Moreover, the form of (1.24) will not be changed by replacing air over the superficial layer 
with any medium when reflecting coefficients on the upper boundary take proper expressions 
other than (1.21). 

As far as is discussed in this section, such a general case will be considered. Now (1.24) 
can be rewritten as 


M={1—(AD—BC)/(A’ D’— B’C’)¥/27e-t1 +B) 4 2 {( A A’ )V/2e-t41 4 — (DD) Ve B18 2 
—e7 +B #{ BC’ + CB’ +2(ADA’D’)¥2?—2(AD—BC)'(A/ D’ —B’C/¥?} (3.1) 


Putting 
p=(ait+Bi)H/2 , q=—(ai—B)H/2 , 
X=e'?—(AD—BCO)Y(A’D’—B’C’)Ve-*? ,  Y=(AA’)Ve4—(DD')2e-4 , (3.2) 
P2=BC’+CB’ +2(ADA’D’)¥?—2(AD—BC)¥*(A’D!—B’C’)? , 
(3.1) may be expressed as 
Mae trtbpe( X?—Y2—P) . (3.3) 
If in (3.3) 
i= 0, (3.4) 
the characteristic equation (2.1) must be reduced to 
MEE, 0) =(X—Y)(X+Y)=0 (3.5) 


which is a very simple form. 
Using the last expression in (3.2), (3.4) should be satisfied by the next two conditions: 


BC’=B’C besides A’D’/AD=B'C’/BC . (3.6) 


If the medium above the superficial layer has the same elastic property as that of the 
medium below the layer, one may find 


A=A’, B=—B’, C=—C’ and D=D’. (3.7) 
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In this case, therefore, (3.6) will be satisfied and / becomes zero. This means that disper- 
sive RAYLEIGH waves in a plate being placed in any isotropic homogeneous medium must have 
a characteristic equation of a type similar to (3.5). 


§ 4. Dispersive RAYLEIGH waves in a plate. 

As has been seen in the previous section, of course / in (3.2) becomes zero if air exists 
above and below a plate. In this case, reflecting coefficients have relations (3.7) and expres- 
sions (1.21). Therefore M in (3.3) may be solved into 


ME, 0) =e PAM OE, o)-ME, w) (4.1) 

in which 
MO=X—Y=(e?—e-'”)— A(et—e-) and M@=X+Y=(e'?—e—'?)+ A(e4—e-4) . (4.2) 
MoE, o)=0 and ME, o)=0 (4.3) 


correspond respectively to the characteristic equations of symmetric modes and antisymmetric 
modes of dispersive RAYLEIGH waves. 

In order to construct cylindrical waves generated from a line source of P-wave at z=E, 
the next operation must be done upon plane waves as (1.1): 


ioe << 


Thus displacement-potentials of dispersive RAYLEIGH waves will be gotten by such a calcula- 
tion as 


[¢ or ?]u-o=zi[residues in (4.4)]. (4.5) 


Putting (1.21), (1.22) and (3.7) into (1.23), and performing the operations (4.4) and (4.5), one 
will obtain the following expressions (Tazimn, 1958): 


r me 
[Ai]u0.0=—ni] i Gletste etre) | 
mQ) 0 


4aM © 
ee eee i idee 
inchs iia ee Ae 46) 
in which fe: ciBH _p-ipH 


{et (+B)/2°H__ g-t-(a+B)/2°H) { pt (@—B)/2+H _p-t+(a—B)/2-H) ’ 


and 


= “ 9—t BIT 
[Ci]ao-9=—rif ale tater att) 


20M uQ)oo ’ 
. Be-te# (4.7) 
1) uQaip= 71 == fhe By -iaA piaH 
[Dilu@ <0 char (eft8 Ke ei@ yee ; 
Then putting (4.6) and (4.7) into (1.1), one has 
ie ya Peer tas fe , 
[d]uM-o= —ni) laM® G{ei* A) ee ies ‘ 
uM =0 
(4.8) 


| B : 
[P]uM-0= — ae {oie (A-H/2) semi ise 2/0) 
§ 


uno 
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or applying (2.3) to the calculation of Me™, 


[b]u-9=— one ( tt \( ie a \ cos a(E—" Joos Hie . 


ade NMOS a a sin BH 
[Peco ee (1 _1)/ 1 8 sinaH \"'/ B/E—1 ) cos (wH/2) 
Ter (i c dh a sinBH ) ( 2B/E ae (8H/2) (4.9) 


16 H 
x co E-—— a 
sa ( 7) sin A(z =i 


since 


_ _sinGH_ _/( #/E—1\cos(@H/2) ‘ 
a 2B/E = (8H/2) ’ (4,10) 


when M®=0, 
Following procedures similar to the above, one will also reach 


[b]u(2)-0= — | and ja psa ) sin a(B— > )sin ae) , 


eH UF .e ~ a sinBH 
ee ee Ie Pe B sinaH \~"'/ @/E?—1 \sin (@H/2) 
[P]arc2><0 wala ~)Q ; So ( a i ae (4.11) 


3 Jak tal 
x sin aE 9 eos a(2 = , 


sin BH (fe ie (wH/2) (4.12) 


~ sinf(a—B)H/2}\ 2B/E /sin (BH/2) ° 


since 


when M@=0. 

These results (4.9) and (4.11) have been verified with those derived from the ordinary 
method by solving four simultaneous equations of boundary conditions for finding Ai, Bi, Ci 
ange 7) in’ (11). 


§ 5. Dispersive RAYLEIGH waves ina superficial layer resting on half space absolutely rigid. 


Reflecting coefficients A’, B’, C’ and D’ on the lower boundary of the layer are calculated 
‘in this case as follows: 


B’=2alE)\/A+(a/E(G/E}=—(a/BC , 


having the next relation with each other, 
AD —BC=1. (5.2) 
Between (1.21) and (5.1), further, another relation is to be found, 
BOS B'C . (5,3) 
Putting (1.21), (1.22), (5.1), (5.2) and (5.3) into (3.2), one has the characteristic equation 
My *=P, 
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where 
X=e'?—e-'P , Y=(AA‘)/2(e4—e") , [2=2(BC’+AA/—1). (5.4) 
Because a as well as 8, contained in p and g, may be either real or purely imaginary, 
(5.4) will be divided practically into the next three cases. 


G@) clv,>m) =a, B=. 


an—[ioF £/(E-1)}/ freed FE )}. 4 (0-8 2)/ (48 6) 
Looking back at (1.21), (1.22) and (5.2), one sees 


BC’ =—{(1—A?)(1—A”)}? , 


because c/vs is larger than // 2 in this case. 
Therefore (5.4) may be rewritten as 


x=sin (Ey aca, Y=(AA‘}3 sin {( 8 ye 


’ Dy i 
ji gl AA Hd AAA) F 
(ii) vpc>us; a=—ia, B=B. 
Putting 
B C ) a Bp 
tan d= qo —— and tano’?=—— —— 
EP NS eas 
one has 
A=—e7s and Al == — rtd? 
On the other hand, from (1.21) and (5.1), one has 
BC’ =Bie +94 B ee ) cos 6 cos 0” . 
ay / ef 
Therefore (5.4) may be rewritten as 
x Xai? 
where 
X=e@/2 cog 5 (BH—6d—0’)—e-@n)/2 cos 5 GH+3+0%), 
(5.7) 


u 


X= ettw sin + GH-8—0') 4 e-GOM/ sin + (BH+0-+9)), 


bo 


1,?=sin (0-++0")+{sin? (0-+06’)— sin? (0—d’)}/? => sin 6 cos 0’(c/vs)*. 


(iii) vp>vsc; a=—id, B=— if. 
speed B/Coyy pet Ea} ard DIB) 


wah MEMO} DIG AN 
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Therefore (5.4) may be rewritten as 
Pe ee Banna if eoce ee 
sinh (Ste), Yaoaesinh (Geaarap 
(5.8) 


ae a ae he ae c \4 
cet edu) 44 Tease ‘ea 
AWS eae bie Us 

Keeping (5.4) in mind, one may obtain displacement-potentials through process similar to 
that of the previous section. Results are shown by 


[8u-0=— 27 Te lees a(z—E)(sin 2p+AA’ sin 2q)-+A’ sin {a(2-+ E)-+29} 
—A sin {a(z+E)—2p}}] / {sin 2p+AA’ sin 2q¢+(8/a)(sin 26— AA’ sin 2q)} , 


and (5.9) 


Aes emt) [Bicos (aE +-Rz—2p)+- A’ cos (aE —fz-+2q)}-+B’ {cos (aE Bz) 
eLEN Ue & 


+A cos (ak — g2)}1 | {sin 2p+ AA’ sin 2q¢+(8/a)(sin 2b— AA’ sin 2q)} . 
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Transition from Solid to Liquid Superficial Waves in a Plate. 


By 


Kyozi Tazime 


Department of Geophysics, Faculty of Science, Hokkaido University, Sapporo. 


Summary 


Even if POISSON’s ratio reaches 0.5, neither M@,™ nor M,@ can coincide with the first 


order of superficial liquid waves. 


respond respectively to sinaH=0 or sin BH=0. 


sion-curves of liquid waves. 


On the contrary, amplitudes must be zero on the latter. 


Each branch of dispersive RAYLEIGH waves should cor- 


The former is not other than the disper- 
Not 


always does the first order of dispersive RAYLEIGH waves correspond to the first order of 


sin aH=0. 


It must be noted that notations and expres- 
sions in this paper will follow those employed 
and obtained in the previous paper by the 
present author (Tazime, 1959). 


$1. Superficial waves in a liquid plate. 
It may be difficult for a liquid plate to 
TEE ; 
E 
z's 
H Plate 
zZ 


Fig. 1. Space-coordinates under consideration. 


Ey ee iY Pd 
[d]u.0=— Gs p JQ mi 


exist in actual fact, but one can consider it. 
Taking space-coordinates as shown in Fig. 1 
and assuming a line source at z=£, displace- 
ment-potential of superficial waves can be 
written as 


[d]uco= 


2rw ( il 
a?H \U 
in which the characteristic equation is given 
by 


l )sin Fe ae ala 
(by 


MO=0 that is sinaH=0. (1.2) 


Dispersive RAYLEIGH waves in a solid 
plate. 

In a solid plate, as well known, M® and 
M©2-waves must exist. Displacement-poten- 
tials of them can be expressed as follows. 


ch) 


§2. 


= vy > Vs; a=a, B=B. 


Sina Ve ae a( = a) 
aa) cos a ( 7) cosa aaah 


a sinBH 


od ae pb sinaH \-* 
aera mite (G2) Bay 


lg 1 eos (at cos a "> )sin A(- 2) 


cos (8H/2) 


in which the characteristic equation is given by 


MO=0 that is 


And 


sin 3(a+8)H+A sin }(a—8)H=0. 


@Q.2) 
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Qnw fl 1 B 

Cay ee | a 

[¢] 1-0 ala i ( es 
~2@ {1 1 

2.9 =1— — Hh 
a Sara c NC 
in which 
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s aH = 
ae sin a( 


sin BH 


HY Se H 
a sina (2 Sak 


sin (a H/2) 


M@=0 


that is 


a sinpBH 


ee sin ott) Hee! 


sin 3(a+$)H—A sin }(a—f)H=0. 


sin (@H/2) sin a( E =f cos A(z = a 


(2.4) 


A in (2.2) and (2.4) means PP-reflecting coefficient as to displacement-potential on the 


ae ed C/Vp 
2 


3 


iaeek DS T 


4 


Si _ it ah ao 


— 


<= 
< 
SS 
e) 2, 
2% 
O20 
Q50 
L <s 
O36 
(~ a4 
O4as 
-|~— 048 = = 


oT emer 


cos (wH/2)=0 


PP-reflecting coefficient on the boundaries. 


or 


boundaries of the plate. 

Letting rigidity of the plate tend 
to zero, A above mentioned will, as 
illustrated in Fig. 2, approach to —1, 
whichever value phase-velocity may 
have. In this case, (2.2) and (2.4) 
will be reduced to 


cos (a@H/2) sin (8H/2)=0 

for M“-waves , 
sin (@H/2) cos (8H/2)=0 

for M©®-waves. 


(2.5) 


Therefore the dispersion-curves for 
M“©-waves will be given by 


sin (@H/2)=0. 


The former of them was contained in (1.2), that is M“=0, and one sees that 


[d]u-0=[b]uc-~0 


and 


[P]ac>-9=0 : 


On the other hand when the latter is satisfied, one sees that 


sin (@H/2)=0 


[¢]u-0=[¢]u-0=0 , 


because {1—(@/a)(sin wH/sin BH)}-! in (2.1) becomes zero. 
Next, the dispersion-curves for M“-waves will be given by 


or 


cos (8H/2)=0. 


The former of them, again, was contained in (1.2) and one sees that 


[P]u2)-0= [db] uce<0 


whichever value z may have. 
On the other hand when the latter is satisfied, one sees that 


and 


[P]u@.0=0 , 


[o] uo =[P]u@-0=0 , 


because sin 6H in (2.1) again becomes zero, 


a ee 
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Thus letting rigidity of the plate tend to zero, branches corresponding to cos (@H/2)=0 and 
sin (@H/2)=0 alone can survive, otherwise amplitudes of displacement-potentials become zero. 


The surviving two branches coincide respectively with the branches of the liquid waves of 
the odd order and the even one, as seen in (1.2). 


Birr cS oy nein, PF 


eno fl 1 8 sinh dé “1 ‘ 
[O]u@)9 = i (= ats & Ee ) cosh a( E- oy )eosh a(2— =) ; 


c sin BH 
[P]war—» = -e(s- Flee Ls Sah B/E 1 ) cosh (aH/2) (2.6) 
dedal GE Ye a sinBH 2B/E cos (BH/2) 
x cosh a(E—7 )sin B(z— =) 
in which the characteristic equation is given by 
M®=0 that is e@#/gin ee Sips aed sin Gee +8)=0. (2.7) 
And 
oe an a | B sinh @H \-*..,, Er Gpe if 
@2)_4= — = = eas et 
[bd] u2)~0 ae - Te & sin BH ) sinh a( E 5 )sinh a( 2 5) : 
lees 6 =i 2nw (5- ‘(- B sinh aH ee ) sinh (aH/2) (2.8) 
@WH\U ¢ a@ sinBH 2B/E sin (8H/2) 
x sinh a( B— y )£08 d= ) 
2 2 
in which 


M@=0 that is etH/2 cos (FO —8)—e#l eos os +3)=0. (2.9) 


0 in (2.7) and (2.9) means phase-lag of displacement-potential at PP-reflection on the 
boundaries of the plate. 

Letting the rigidity of the plate tend to zero, the above mentioned 0 will, as illustrated in 
Fig. 3, approach to zero within the very wide region near c=v»y. In this case, (2.7) and (2.9) 
-will be reduced to 


0.4} g » 
fe) fe} 
0.2} a a3. 
EE / 
°°60 Ol 02 03 08 09 10 ¢/v, 


Fig. 3. Phase-lag at PP-reflection on the boundaries. 
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cosh (aH/2) sin (BH/2)=0 for M-waves , (2.10) 


sinh (@H/2) cos (BH/2)=0 for M@-waves. 


Because cosh (aH]/2) cannot be zero, the dispersion-curves for M“-waves must be sin (8H/2) 
=0. Putting this into (2.6), one sees that 


[d]uD-0=[P]u@-0=0 . 
On the other hand, the dispersion-curves for M®-waves will be given by 
sinh (@H/2)=0 or cos (8H/2)=0. 


The former of them makes c equal to vp for any value of EH and one sees that the right 
hand sides of (2.8) become zero at the same time, because U must be equal toc in this case. 
When the latter is satisfied, one sees that 


[d]uv@-0=[P]u@-0=0 , 
because sin 8H becomes zero. 
(iii) vp>vs>c; a=—ia, B=—if. 


A = A way " r 
[d]e@o = one Ge = \a- Le SS) cosh a (S) cosh a iS =) 7 


w@H\U ¢ a sinh BH 2 
Chen ( ae )a- Zou Ay P/E + ie (aH/2) (2.11) 
w@H\U ¢ a sinh BH 28/E / cosh (8H/2) 
x cosh @ (E-¥ sinh B ee a 
Zz 2 
in which the characteristic equation is given by 
M®=0 that is sinh }(a+8)H+A sinh }(a—f)H=0. (2.12) 
And 
; Qnw (1 1 8 sinhaH\-!_.. + Hotes H 
wan Ba hE BBR sn (eH) ani (e 4 
u@)20 ere ee TT sinha|{E 9 sinh a { z arn 
oy ee dno Poul y(a- 6 sinh aH i" ( P/E +1 sinh (@H/2) (2.13) 
CHC, a sinh BH 28/E / sinh (@H/2) 


x sinh a (E- =) cosh B (a5 4 
2 2 
in which the characteristic equation is given by 


M®=0 that is sinh }(a+8)H—A sinh (a—8)H=0. (2.14) 


A in (2.12) and (2.14) means again PP-reflecting coefficient. 


Letting rigidity of the plate tend to zero, A must be confined to —1, as illustrated in Fig 
; ; 


In this case, (2.12) and (2.14) will be reduced to 
cosh (aH/2) sinh (BH/2)=0 for M-.waves , 


sinh (@H/2) cosh (@H/2)=0 for M@-waves . (2.15) 
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0.5, A in Fig. 2, being c/up 
larger than unity, must be —1 
for any value of c/v». 

Turning attention to A for 
general Porsson’s ratios, one sees 
that A will always approach 
again to —1 if c/vpS+1. There- 
fore (2.5) must be satisfied as 
before, for any o if c/vp>1. 

Considering the above two cir- 
cumstances with Fig. 2, it may 
be expected that the nearer o 
approaches to 0.5, the wider be- 
comes the range satisfying (2.5). 

If c/vp=1, (2.5) must be satis- 


: fied again, because A=—1 for 


048 
O45 
0.40 
035 
030 
025 


¥ 
-30+ | u 
| 
| 
| 


a 


(ofofe) 


any value of o. As a=0 in this 

case, sin(8H/2) must always be 

at zero. Thus one has 

EH/2=(mn/2){(vp/vs)—-1}-?, 

m being even for M@-waves. 
(3.1) 

For M@®-waves, on the other 


hand, such a relation as (3.1) 
cannot always be obtained, be- 


cause sin (@H/2) may be zero in 


Fig. 4. PP-reflecting coefficient on the boundaries. (2.5) 


As one sees that 


sinh (@H/2) cosh (a@H/2) cosh (H/2) 
cannot be zero , 


the branches corresponding to sinh (8H/2)=0 
can alone exist, with the result that the 
-right hand sides of (2.11) become zero at the 
same time. Thus one has also 

[¢]w@-0=[$]a1@~0=[b]1@-0=[P] u@-0=0 A 

It has been recognized in this section that 
dispersive RAYLEIGH waves in a solid medium 
should be reduced to superficial waves in a 
liquid at the limit when rigidity is taken as 
zero. In the following section, processes of 
transition will be investigated. 


§3. Numerical calculations of the disper- 
sion-curves for various POISSON’S ratios. 


At the limit when Porsson’s ratio arrives at 


If c/vs=1, (2.10) must be satisfied, because 
6 in Fig. 3 is zero for any value of o. As 
cosh (aH/2) sinh (aH/2) cannot be zero in this 
case, it may be expected that the dispersion- 
curves of M™ and M®)-waves will coincide, 
for EH>1, respectively with sin (@H/2)=0 
and cos(8H/2)=0. At c/us=V/ 2,6 becomes 
always 7/2, therefore the dispersion-curves 
of M® and M®-waves must coincide re- 
spectively with cos (@H/2)=0 and sin (8H/2)=0 
at a point where vp > ¢ > vs. 

As one has seen by now, for any value of 


sina~H=0 and sinfH=0 (372) 


might be the bases with respect to dispersion- 
curves.of dispersive RAYLEIGH waves. 

In Fig. 5 are shown the curves given by 
(3.2), that is 
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Pe ei " M,™ alone, however, must be — 
oc is 5 excluded from the systematic 
pe ad curves above described. Excluding 
™ ae M,@®, all expected matters have 
vii | in practice been also justified as to 
| | M,,)-branches in Fig. 6. 

? ! In Figs. 7 to 10, the dispersion- 
ea | curves for several values of Pots- 
$ son’s ratio are illustrated in the 
| 10 = ane same manner as that in Fig. 6. 

i= \ f 
IE | i ; 
wf H $4. Some noteworthy points. 
1 / i 
_e ; / As M,,-waves could be discussed 
03/- Me similarly to M,™, the latter alone 
E Cees eee, C=Vs will be discussed here. 
% (1) If c/vpyS>1, the dispersion- 
hone ee inl ‘wiaskes | | curves always coincide with any 
Oi 0203 0507 10 +2 3 +#§ 7 10 «2 30 One of the curves given by (2.5), 
whatever value o may have. 
mara: (2) If c/vpS>1, periods of the 
Fig. 5. Chain-lines correspond to @H=Ix, while broken and former and of the latter equation 


dotted lines to BH=mx. T=2r/w. 


aH=Ix and BH=mzx (3.3) 


where / and m represent positive integers, 
indicating orders of these curves. The curves 


given by sin 8H=0 must be shifted to right- 
ward with increase in o, though the curves 
given by sin aH=0 need not. In Pigs 5, for 
an example, the curves for o=0.48 are ex- 
hibited. Deducing from (2.5), odd orders of 
/ and even orders of m will belong to M®, 
On the contrary, even orders of 7 and odd 
orders of m will belong to M@-waves. 

Then the dispersion-curves for o=0.48 were 
calculated, as shown by full lines in Fig. 6 
where the other lines correspond respectively 
to those in Fig. 5. 

All matters expected at the beginning of 
this section have appeared satisfactorily in 
Fig. 6. 

At c/vp=10, M3™ coincides with J=1 while 
Mi, M: and Mi™ coincide respectively 
with m=2, 4 and 6. 

At c/vp=1, all M,@-branches coincide with 
[rv 1%, (6), 1s) hovel Xo) caval, 

At c/vs>1, all Mn™ coincide again with 
the above m’s. 


in (3.3) should coincide with each 
other under the next condition, 


m=(vp/vs)l . (4.1) 


Because even orders of m must correspond 
to odd orders of /, the lowest order of m 
satisfying (4.1) is 2 for /=1, resulting in 
Vp/Ve=2 Or a=0.33. 

If « < 0.33, no curve indicated by m exists 
to the right of the curve indicated by /=1 
with which M,™ always coincides in this 
case. 

If o>0.33, on the other hand, some curves 
indicated by m appear to the right of /=1. 
The nearer o approaches to 0.5, the larger 
becomes the number of m-curves to the right 
of /=1 where the order of m is equal to m/2. 

(3) On dotted and broken lines in Figs. 6 
to 10, as investigated in the second section, 
amplitudes of displacement-potentials are 
zero. Though the dispersion-curves may be 
calculated, no wave can exist there in prac- 
tice. 

(4) It may be expected, on the other hand, 
that on chain lines, amplitudes will be large. 
Indeed amplitudes on full lines differ from 
zero, but the maximum amplitude, being in- 
finitely large, must be attained on chain lines at 


} 
‘ 
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Fig. 8. Dispersion-curves for o=0.33. 


Tv,/H=2/1 . (4.2) 


Finite amplitudes on full lines are trivial in 
comparison with infinitely large one. 

(5) It has been found to be wrong, if 
o>0.33, to say that the lower the order of 
dispersive RAYLEIGH waves is, the larger may 
become its amplitude. 

(6) The order of superficial waves, such as 
Love waves or liquid waves, consisting of 
either S or P-waves alone indicates the num- 
ber of nodes in the layer at the same time. 
However M3 for o=0.48 at c/vpS1, as 
above described, should correspond with the 
liquid wave of /=1 which has only one node 
in the superficial layer. One will see more- 
over that M,™@ for o=0.48 has not four 
nodes. Thus the order of Mn™ and M,,@)- 
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m=2 
t= 


MS 
‘S 


al 02 03° OS O7 10 2 3 
——> Tvp/H 


Fig. 9. Dispersion-curves for o=0.25. 


waves cannot always indicate the number of 
nodes in the layer. 

(7) (4.2) means “half wave-length law”. 
It was a wonder that such a law as (4.2) had 
been satisfied for dispersive RAYLEIGH waves 
in practice, in spite of the existence of P- 
waves as well as of S-waves in the layer. 
But it has been understood in the present 
paper that among many branches of disper- 
sive RAYLEIGH waves a few of them must 


correspond to sinwH=0 and others to sin BH 
=0. Amplitudes of the latter ones cannot be 
so large that the law containing v, alone will 
become important. It must be also noticed 
that the order of dispersive RAYLEIGH waves 
cannot be determined from the “ wave-length 
law ” observed. 
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Fig. 10. Dispersion-curves for =0.00. 
(8) Picking out M1@, M,@ and the curve 


for /=1 alone from Figs. 6 to 10, one will 
obtain Fig. 11. In this figure MM and M,™ 

~ diverge from the curve for /=1 when o ap- 
proaches to 0.5. Thus it should not be ex- 
pected that the two curves, for solid and for 
liquid waves, will coincide with each other 
at the limit of c=0,5. The process of transi- 
tion from solid to liquid cannot be understood 
by Fig; 11. 

According to the present investigation, one 
sees, the process is very complicated. When 
o is near 0.5, the first order of liquid waves 
is constructed approximately of many higher 
orders of M,-waves and the second order 
of liquid waves is constructed approximately 


10 : 1 = 
rt Mm! 
S- 
37 
| 
ea ok 7 
3 ee") 
10 ——— Li See 
} is 0.0 fe epee? 
O7 025 A 
pee 
05-033 
(ee) 
03+ 
O=0A8 a 
lor | i | —— 
To \= 
[ | | M2) 
Site | 
3} | 
oy 
so 2F ey he 
Lor == La ee =F: ] == 
= O00 a, ie 
07 G25 ae Si Yh ; 
OSF Lae a a 
040 
0.3);- 
o=048 
itil | tein 
a2)> 03°05 07, 10 2 3 Salo) 20 
——> T/tl 
Fig. 11. Comparison of dispersion-curves for the 


first order of M-waves to the curve for the 
first order of liquid waves. 


of many higher orders of M,‘-waves and so 


on. 


The present author thinks that conclusions 
reached in this section as to waves in a plate 
may be somewhat common to general disper- 
sive RAYLEIGH waves in superficial layers. 
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